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A  valuable  figure  of  merit  in  evaluating  and  comparing  analytical  techniques 

is   the   limit  of  detection   (LOD),   which   represents   the   minimum  detectable 

concentration  or  amount  of  analyte  in  a  sample.   The  factors  which  influence  the 

value  of  the  LOD  are  theoretically  evaluated  through  the  application  of  estimation 

theory  and  propagation  of  errors  to  the  LOD  concept.  Equations  are  derived  which 

can  be  used  to  estimate  the  magnitude  of  the  random  fluctuations  which  result  from 

the  use  of  sample  statistics  in  the  estimation  of  the  LOD.  The  resulting  confidence 

intervals  can  be  used  to  evaluate  and  compare  the  true  LOD  values  of  analytical 

systems,  and  to  determine  the  most  efficient  procedure  for  quick  estimation  of  LOD. 

Monte  Carlo  simulations  of  typical  analytical  situations  are  used  to  evaluate  the 

effectiveness  of  the  derived  equations. 
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Application  of  conventional  detection  limit  theory  is  not  straightforward  when 
considering  laser  spectroscopic  methods  which  are  capable  of  detecting  single  atoms 
or  molecules  in  the  laser  beam.  Theoretical  considerations  in  the  detection  and 
evaluation  of  these  methods  are  addressed  based  on  a  model  of  a  typical  laser 
spectroscopic  experiment  with  destructive  and  nondestructive  detection  methods. 
Computer  simulations  are  used  to  verify  the  application  of  the  theory  of  single  atom 
detection  (SAD)  to  typical  experimental  situations,  as  well  as  to  discuss  the  scope  of 
SAD  methods  and  the  various  possible  signal  processing  methods  which  can  be  used 
to  continuously  monitor  and  count  the  atoms  or  molecules  which  flow  through  the 
laser  beam. 


vu 


CPIAPTER  1 
INTRODUCTION 


Any  analytical  technique  designed  to  determine  the  concentration  or  amount 
of  analyte  in  a  given  sample  from  the  magnitude  of  the  resulting  signal  can  be 
characterized  by  a  limit  of  detection  (LOD).  The  LOD  is  designed  to  give  an 
mdication  of  the  lower  limit  of  analyte  concentration  or  amount  that  the  given 
analytical  technique  can  distinguish  from  the  background  noise,  which  is  present  even 
when  the  analyte  is  absent.  The  LOD  is  very  often  a  significant  characteristic  of  a 
given  analytical  procedure,  since  calculated  LODs  can  be  used  for  comparison  and/or 
evaluation  of  the  procedure  relative  to  other  analytical  procedures.  For  example,  the 
LOD  may  be  used  to  indicate  the  improvement  of  the  detecting  power  of  a  given 
analytical  protocol;  the  improvement  is  measured  by  comparing  the  LOD  to  the 
values  reported  for  the  procedure  in  the  past.  Alternatively,  a  particular  application 
may  have  certain  sensitivity  requirements;  in  this  case,  the  reported  LOD  value  may 
be  used  as  a  first  basis  of  evaluation  as  to  the  suitability  of  the  method  for  the 
application. 

The  importance  of  the  LOD  to  characterize  a  given  method's  detection 
abilities  signifies  that  great  care  must  be  taken  when  estimating  and  reporting  LODs. 
The  actual  value  for  the  LOD  calculated  for  a  given  analytical  procedure  may  vary 
widely  due  to  either  of  two  factors:  (1)  the  method  used  to  measure  and  calculate  the 


LOD;  and  (2)  the  use  of  different  definitions  to  characterize  a  minimal  detectable 
concentration/ amount  of  analyte.  Despite  the  effort  made  in  the  past  two  decades 
to  define  the  LOD  unambiguously,  and  to  recommend  guidelines  for  its  measurement 
[1,  2],  this  diversity  in  definition  and  measurement  protocol  persists. 

The  random  variation  in  calculated  LOD  is  often  intuitively  sensed  by 
practicing  analytical  chemists,  who  realize  that  LODs  which  are  close  (within  a  factor 
of  2  or  3)  may  not  be  statistically  different.  For  example,  if  the  same  chemist 
determines  the  LOD  of  the  same  analytical  technique  two  consecutive  times,  it  is 
usually  recognized  that  the  same  value  for  the  LOD  will  not  result.  Additionally,  if 
the  LOD  is  measured  on  the  next  day  by  a  different  chemist,  then  it  may  be  even 
more  likely  that  another  value  will  be  obtained.  This  variability  in  calculated  LOD 
is  usually  taken  into  account  only  in  a  general  sense  when  comparing  LODs  from 
different  methods;  i.e.,  two  LODs  which  are  close  are  often  considered  equivalent. 
However,  it  would  be  desirable  to  know  just  how  much  might  be  of  the  variability  in 
calculated  LOD  values  might  be  due  to  random  fluctuations  in  measurement,  and 
how  much  due  to  actual  differences  in  analytical  conditions  (sensitivity,  background 
noise,  etc). 

The  various  different  definitions  used  for  the  LOD  were  recognized  and 
reconciled  in  the  pioneering  works  of  Kaiser  [3]  and  Currie  [4],  and  the  situation  has 
greatly  improved  since  that  time.  However,  there  are  now  emerging  laser 
spectroscopic  methods  which  claim  the  ability  to  detect  very  small  numbers  of  analyte 
atoms  or  molecules  within  the  laser  beam,  all  the  way  down  to  the  level  of  single 
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atoms/molecules.  The  application  of  the  LOD  concept  to  these  newer  methods  is 

not  always  obvious,  but  the  comparison  and  evaluation  of  these  methods  needs 
consistent  and  appropriate  definitions  specifically  designed  to  address  the  difficulties 
which  result  from  the  ability  to  detect  small  numbers  of  atoms. 

The  intent  of  this  dissertation  is  to  address  both  of  the  subjects  discussed 
above.  A  review  of  some  of  the  relevant  concepts  and  past  literature  is  presented 
in  chapters  2  and  3;  these  chapters  form  the  basis  for  the  new  work  presented  in  the 
remainder  of  the  dissertation.  The  variability  of  the  calculated  LODs  will  be 
addressed  from  a  theoretical  standpoint  in  chapter  4;  in  chapter  5,  logical  and  precise 
definitions  will  be  presented  which  are  designed  to  clarify  misunderstandings  which 
may  arise  when  attempting  to  define  detection  limits  and  other  figures  of  merit  for 
laser-based  analytical  methods  capable  of  detecting  single  atoms/molecules.  Later 
chapters  serve  to  confirm,  evaluate,  and  illustrate  the  concepts  introduced  in  chapters 
four  and  five  through  the  use  of  simple  Monte  Carlo  computer  simulations. 


CHAPTER  2 
THEORY  OF  LIMITS  OF  DETECTION 


The  limit  of  detection  (LOD)  is  an  analytical  figure  of  merit  (FOM)  which 
gives  the  minimal  concentration^  of  analyte  in  a  sample  which  can  be  distinguished 
from  the  blank.  The  LOD  of  a  procedure  is  only  one  of  several  possible  FOMs. 
The  purpose  of  FOMs  is  to  objectively  summarize  important  characteristics  of  the 
analytical  procedure.  Other  than  the  LOD,  there  are  a  number  of  important  FOMs, 
such  as  the  linear  dynamic  range,  the  sensitivity  of  the  technique,  the  resolution,  the 
accuracy  and  precision  of  determination,  the  informing  power  of  the  method,  price 
and  time  of  analysis,  interferences,  and  others  which  are  less  commonly  used.  The 
use  and  importance  of  many  of  these  has  been  covered  in  reviews  [5-7]  and  books 
[8-11]  and  will  not  be  covered  here  other  than  to  note  that  the  LOD  is  only  one  of 
a  group  of  FOMs  which  describe  the  total  analytical  procedure,  although  the  LOD 
is  one  of  the  most  important  FOMs. 

The  present  concept  of  the  LOD  has  been  promoted  in  several  landmark 
papers  [3,  4,  12].  These  works  mark  the  begirming  of  applying  a  statistical  approach 
to  the  problem  of  determining  the  minimum  detectable  concentration  of  analyte  for 
a  given  analytical  procedure.  This  chapter  will  summarize  all  the  relevant  principles 


^for  convenience  the  term  "concentration"  will  signify  either  analyte  concentration  or 
amount,  whichever  is  appropriate. 


5 
of  detection  limits  in  both  the  signal  domain  and  in  the  concentration  domain.  There 

will  be  no  attempt  at  a  comprehensive  review  of  past  literature;  these  have  been 

presented  in  various  review  articles  and  textbooks  [8,  10,  11,  13-16]. 

Analvte  Signal  Detection 

Signal  Detection  Limit 

Any  instrumental  analytical  method  is  indirect  in  the  sense  that  the  "result"  of 
a  single  analysis  is  a  mean  instrument  response  in  the  signal  domain  (e.g.,  current, 
charge,  potential  difference)  rather  than  a  direct  reading  of  the  analyte  concentration 
in  the  sample.  In  a  steady-state  measurement,  the  response  of  the  instrument  for  a 
given  sample  is  measured  for  a  certain  time  Tn,,  and  the  average  response  during  this 
time  is  taken  as  a  single  measurement  of  the  sample.  The  relation  of  a  given 
measurement  to  the  concentration  of  analyte  is  found  through  a  calibration  of  the 
response  function  of  the  instrument;  for  the  present  time,  however,  we  will  only  be 
concerned  with  the  instrument  response  in  the  signal  domain. 

In  the  absence  of  analyte  in  the  sample,  there  may  be  a  nonzero  response 
during  T„  due  to  the  blank  only.  This  response  consists  of  a  nonrandom  component 
fj^  and  random  component  characterized  by  the  standard  deviation  a^.  Although  it 
is  certainly  possible  to  compensate  for  the  non-random  contribution  of  the  blank,  the 
random  fluctuation  ct,,  will  still  contribute  to  uncertainty  in  the  measured  signal.  For 
signals  which  are  of  the  same  magnitude  as  a^  there  is  a  need  for  some  criterion  to 
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decide  when  a  given  signal  is  due  to  the  presence  of  analyte  in  the  sample  or  when 
the  measurement  may  be  due  to  spurious  fluctuations  due  to  the  background  signal. 

The  minimal  detectable  signal  in  the  steady  state  case  is  best  illustrated 
through  a  simple  example.  A  typical  situation  is  shown  in  figure  1.  The  response  as 
a  function  of  time  is  shown  for  one  measurement  (where  T^^  =  1000  s  in  this  case). 
Each  individual  value  shown  in  fig.  1(a)  fluctuates  about  the  (unknown)  mean  (m^  with 
variance  a^;  the  distribution  of  the  values  (population  distribution)  is  shown  in  fig. 
1(b).  In  the  particular  case  shown  in  this  figure,  the  population  of  individual 
readings  in  the  blank  is  normally  distributed  with  a  mean  of  50  mV  and  a  standard 
deviation  of  10  mV.  The  average  value  in  fig.  1(a)  during  T^  is  50.259  mV;  this  is 
one  measurement  of  the  blank. 

According  to  the  central  limit  theorem  [17],  the  blank  measurements  will  be 

approximately  normally  distributed,  no  matter  the  form  of  the  original  distribution 

■v 
of  individual  values,  with  a  standard  deviation  given  by 

■ "'  »< 
Ob-  —  [2-1] 

where  ;'  .►^ 

N  =  the  number  of  individual  measurements  during  T^,  ■ .  ^ 

a^  =  the  standard  deviation  of  the  individual  signal  values,  and 
a^  =  the  standard  deviation  of  the  blank  measurement.  '■•/ 

The  above  equation  holds  if  the  fluctuation  of  measurement  values  due  to  the  blank 
is  characterized  by  "white"  noise,  rather  than  long-term  drift  (i.e.,  fUcker  noise).  The 


(I-H) 


.•.•r»*i"i. 


•     'I  l^i; 


> 

8 


09 


O 

O 


a 


t 


■a 


a 
6 

s 


2 

(d 
•o 

o   c  (U 

^    «  J3 

„    On  *' 

II    in  a 

2    O  en    g 

9  « 

.2  13  3  o 

.•s  p  13  ^ 

"H     bo   >     CO 

CJ     t«  — ;     _ 

O     _  t«     C5 

(u   a  c  'Ti 

en    3  o  T3 

S    *  3   iS 

j^   p  C   <u 

a  ^  tS  ■S 

^•^  §  S 

«  .9  3  ,„ 

U ^  &"•- 


o 

10 


10 

CM 


8 


10 
N 


o 

10 


to 


(AID)    IBUSIS 


s 

00 


:  '.-f-     ■%■ 


8 

question  now  becomes,  at  what  point  is  a  sample  measurement  said  to  be  "detected" 

above  the  blank  noise.  In  other  words,  we  seek  a  measurement  value  large  enough 
so  that  the  chance  of  the  value  belonging  to  the  distribution  of  possible  blank  values 
is  negligible. 

The  lowest  value  at  which  the  measurement  is  considered  to  be  due  to  a 
process  other  than  blank  noise  is  the  signal  detection  limit,  X^.  Since  the  standard 
deviation  of  blank  measurements  directly  limits  the  abiUty  to  detect  small  signals,  the 
detection  limit  is  directly  related  to  the  value  a^: 

X^  =  Hi  +  ka^  [2.2] 

where  k  is  known  as  the  confidence  factor.  The  probability  that  a  blank 
measurement  can  give  rise  to  a  measurement  value  greater  than  or  equal  to  the 
detection  limit  is  known  as  the  type  I  error,  or  the  probabiUty  of  false  positive  error: 

PiX.^Xj)   -  a  [2.3] 

where 

a  =  the  probability  of  type  I  error,  and 

Xb  =  one  blank  measurement  (the  mean  during  T„). 
The  value  of  a  is  controlled  by  the  confidence  factor  chosen  in  eqn.  2.2. 

Choosing  the  Confidence  Factor 

According  to  the  central  limit  theorem,  if  T^  is  long  enough  it  is  usually 
reasonable  to  assume  that  the  blank  measurements  follow  a  normal  distribution  N(/ib, 
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a^^).  In  terms  of  the  z-statistic, 


^b-\^b 


z — ^  ~  N{0,1)  [2.4] 

where  the  symbol  ~  means  "distributed  as."  With  this  in  mind  the  confidence  factor, 
k,  can  be  chosen  according  to 

k  -  z,/2  [2.5] 

where  z^/^  is  chosen  from  tables  of  the  z-distribution.  The  factor  k  depends  on  the 
desired  one-sided  confidence  level,  given  by  the  probability  (1-a). 

The  detection  limit  is  defined  in  eqn.  2.2  in  terms  of  the  population 
parameters  n^  and  a^.  However,  these  parameters  are  usually  unknown  and  must  be 
estimated  by  repeated  measurements  on  the  blank.  The  effect  of  substituting  an 
estimate  s^^  for  a^  in  the  z-statistic  is  to  broaden  the  distribution.  In  this  case,  the 
t-statistic  must  be  used,  with  n-1  degrees  of  freedom,  since 

^  ~  V.  [2.6] 

4 


where  n  is  the  number  of  blank  measurements,  each  for  time  T^,  used  to  calculate 
Sb  .  In  addition,  the  effect  of  using  the  estimated  mean  X^  for  /i^  in  the  above 
statistic  can  be  seen: 


Sb  +  Si/n 


tn-2  [2.7] 
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where  the  denominator  reflects  the  increased  variance  in  the  numerator.  The  above 

expression  suggests  that  the  confidence  factor  should  be  chosen  according  to 

^-^n-2.«/2(l^l/«)'^'  [2-8] 

Thus,  the  signal  detection  limit  can  be  found  by  first  choosing  the  desired  confidence 
level  a,  making  repeated  measurements  on  the  blank  signal  (at  least  16-20 
measurements  are  recommended),  calculating  the  sample  mean  and  standard 
deviation,  and  substituting  in  the  above  equation  to  find  the  value  of  k  to  use  in  eqn. 
2.2.  In  practice,  the  (1/n)  term  in  eqn.  2.8  is  usually  ignored  and  the  t-factor  is 
calculated  with  n-1  degrees  of  freedom.  In  essence,  this  is  the  same  as  ignoring  the 
effect  of  using  an  estimate  for  the  parameter  ^i,;  however,  for  typical  values  of  n,  the 
effect  of  ignoring  the  1/n  term  in  eqn.  2.8  is  small. 

The  value  recommended  by  lUPAC  for  the  confidence  factor  is  3  [18];  with 
/Jb  and  a^,  known  this  would  result  in  a  =  0.0014.  Of  course,  even  assuming  a  normal 
distribution  of  blank  measurements,  the  true  value  of  a  would  increase  due  to  the 
imprecise  nature  of  the  estimates  used  for  the  population  parameters.  A  non-normal 
distribution  might  inflate  the  type  I  error  probability  even  further,  although  by 
Chebyshev's  theorem  [19]  this  probability  cannot  exceed  0.11  (for  k=3)  when  the 
blank  noise  is  due  to  random  error  alone.  A  recent  publication  points  out  that  the 
true  value  of  a  cannot  be  accurately  known  due  to  various  factors  such  as  systematic 
error,  long-term  drift,  and  non-normal  distributions  of  blank  measurements  [2].  Since 
strict  interpretation  of  X^  and  LOD  in  terms  of  a  is  usually  difficult,  a  value  of  k=3 
was  recommended  for  consistency. 
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The  detection  limit  for  the  example  given  earlier  (fig.  1)  is  shown  in  fig.  2. 
Figure  2(a)  shows  X^  in  terms  of  the  probability  distribution  of  the  blank 
measurements,  X^,  with  k  =  3.  In  order  to  get  a  better  idea  of  the  magnitude  of  the 
detection  hmit  in  comparison  with  the  background  fluctuation,  figure  2(b)  shows  a 
signal  at  a  level  above  the  detection  limit  (55.102  mV),  compared  with  a  single  blank 
measurement  (50.259  mV).  The  signal  shown  is  barely  distinguishable  from  the 
blank  measurement  by  eye. 

There  are  two  important  points  which  should  be  made  concerning  the  above 
discussion  of  X^: 

(1)  If  a  single  "measurement"  actually  consists  of  two  measurements  -  one  on  the 
sample,  and  one  for  blank  subtraction,  then  the  confidence  factor  should  be 
multiplied  by  a  factor  of  v2.  Simultaneous  blank  subtraction  is  frequently  used  to 
account  for  long-term  drift  in  n^^. 

(2)  The  discussion  above  was  for  steady  state  signals.  The  value  of  the  detection 
limit  will  depend  on  the  value  of  T^  chosen  and  so  this  value  should  always  be  given 
as  part  of  the  experimental  procedure.  Application  of  the  above  theory  to  the  case 
of  transient  signals  is  straightforward  in  the  case  of  peak  detection,  and  o^  is  used  in 
the  above  equations  instead  of  a^^.  The  correct  detection  limit  value  in  the  case  of 
peak  integration  is  somewhat  less  obvious,  and  will  not  be  discussed  here. 
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Minimum  Detectable  Concentration 

'  ,  The  usefulness  of  the  particular  value  of  X^  for  an  analytical  procedure  is 
hmited.  Although  knowledge  of  X^  is  necessary  for  any  analyst  who  wishes  to  detect 
the  presence  of  a  small  signal  as  an  indication  of  whether  an  analyte  is  present,  it  is 
very  difficult  to  use  X^  as  a  comparison  between  different  methods,  or  between 
methods  reported  in  the  literature.  Therefore,  the  value  X^  must  be  transformed 
into  a  useful  measure  of  a  methods  ability  to  detect  the  presence  of  small  amounts 
of  analyte.  Before  this  can  be  done,  the  analytical  procedure  must  be  calibrated;  i.e., 
the  functional  relationship  between  the  measured  signal  and  the  analyte 
concentration  must  be  known. 

Calibration 

The  most  commonly  used  method  of  calibration  is  by  linear  least-squares 
regression.  Textbooks  on  regression  give  detailed  theory  on  various  types  of 
regression  analysis  [20,  21];  this  section  will  only  cover  the  most  basic  type,  simple 
first-order  linear  least-squares  regression. 

In  many  cases  in  analytical  chemistry,  a  linear  relation  between  the  signal  and 
the  analyte  concentration  can  be  assumed  over  the  range  of  interest,  and  the 
following  model  appUes: 

V-y,  -  "o-^i  -^  \^b  roQi 
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where 

A»yi  is  the  true  mean  response  at  X  =  Xj; 

Uq  and  /ib  are  the  true  slope  and  intercept  of  the  caHbration  Une; 

Yj  is  an  observed  (variable)  response  at  X  =  X;;  and 

6;  is  the  true  error  of  an  observed  response  Y;  (the  residual). 
The  first  equation  describes  the  true  mean  value  of  the  signal  at  a  given  value  of 
analyte  concentration;  the  second  equation  describes  the  effect  of  the  random 
variation  in  the  observed  measurement  at  X  =  Xi.  The  parameter  described  by  Oq  is 
also  known  as  the  sensitivity  of  the  analytical  method. 

Estimates  of  the  parameters  a^  and  n^  are  usually  found  by  using  least-squares 
estimators;  equations  for  these  estimators  and  conditions  for  their  validity  are  readily 
available  [20,  21].  Using  these  estimates,  the  model  becomes 

r,  -  a,X,  ^Xi^e,  [2.10] 

where 

a^  =  the  estimated  slope; 

X,,''  =  the  estimated  blank  response  (i.e.,  the  intercept);  and 

Cj  =  the  observed  residual  at  X  =  Xj. 

Of  course,  the  estimators  used  are  subject  to  variation  o^  and  a^  for  the  slope 
and  intercept,  respectively.  The  equations  for  the  least-squares  estimates  and  their 
estimated  standard  error,  Sg  and  s,,,  are  readily  available  in  textbooks,  along  with 
conditions  for  their  vaUdity  [20,  21];  these  values  are  usually  automatically  computed 
in  regression  software  packages.  The  variability  in  the  estimates  for  the  parameters 
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of  the  linear  model  in  eqn,  2.9  means  that  there  will  be  uncertainty  in  a  given 

predicted  response  value  Yj.  It  is  possible  to  construct  a  confidence  interval  within 

which  the  true  value  of  ^yi  will  lie  with  100(l-a)%  reliability.    This  confidence 

interval  is  given  by 


^i   *   ^a/2.n-2^ 


1 i 

N         S„ 


where 


[2.11] 


5„  -  Y:  (^r^) 


[2.12] 


and 


s  =  standard  deviation  of  response  (assumed  constant);  and 

N  =  number  of  points  in  calibration  curve. 

The  interval  above  is  likely  to  contain  the  true  mean  response  at  a  given 
analyte  concentration.  An  interval  which  describes  where  a  future  response  at  X  = 
Xj  is  likely  to  fall  with  (1-a)  probability  is  often  called  the  prediction  interval,  and 
is  given  by: 


•'j    *    ^a/2.n-2^ 


[2.13] 


The  relation  between  the  two  intervals  for  a  typical  calibration  curve  is  shown  in 
figure  3.  The  wider  intervals  are  the  prediction  intervals. 

Weighted  least-squares  regression.  It  should  be  noted  that  the  above  intervals 
were  derived  with  the  assumption  of  constant  variance  in  the  measurement  response 
along  the  calibration  curve.  In  analytical  chemistry,  however,  there  are  a  number  of 
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techniques  with  large  hnear  dynamic  ranges,  where  it  is  Hkely  that  the  noise  on  the 

signal  increases  with  the  concentration.  An  example  was  given  recently  for  atomic 
emission  in  the  inductively  coupled  plasma  [22].  In  addition,  certain  transformations 
of  variables  have  the  affect  of  skewing  the  error  magnitude  even  if  the  assumption 
of  constant  variance  were  valid  to  begin  with  [23].  In  these  cases,  weighted 
least-squares  estimates  must  be  used,  particularly  when  it  is  important  to  obtain 
information  on  the  sizes  of  the  intervals  given  in  the  above  equations. 

One-point  calibration  curves.  If  it  is  only  desired  that  an  estimate  of  the 
sensitivity  ao  near  X^j  be  required,  then  a  single  standard  can  be  used  so  long  as  it 
is  known  that  the  linear  model  applies  up  to  the  standard  concentration,  and  that  the 
standard  is  reUable.  In  this  case,  the  standard  deviation  of  the  sensitivity  estimate 
is  given  by: 


a,  -  -^  [2.14] 


where 

C  =  standard  concentration,  and 

CTy^g)  =  standard  deviation  of  response  at  C. 
Of  course,  multiple  measurements  of  the  standard  are  needed  to  estimate  o^^y  Note 
that  this  fluctuation  in  response  includes  the  uncertainty  in  both  the  blank  and  the 
signal  measurements. 
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Limit  of  Detection 

Once  the  analytical  system  has  been  calibrated,  it  is  possible  to  define  the 
limit  of  detection,  LOD,  as  the  analyte  concentration  which  corresponds  to  the  signal 
detection  limit,  X^: 

LOD  =  ^^sCiil  .  hlb  nis] 

where  all  the  terms  have  been  previously  defined.  The  value  of  LOD  thus  defined 
can  be  used  for  comparisons  between  analytical  procedures.  Thus,  in  addition  to  the 
estimates  for  the  blank  mean  and  standard  deviation  necessary  to  estimate  X,,,  an 
estimate  for  the  sensitivity  must  also  be  used  to  calculate  the  LOD. 

The  confidence  factor,  k,  in  the  definition  of  LOD  is  chosen  as  described  in 
the  section  on  signal  detection.  However,  the  use  of  the  calibration  estimate  of 
sensitivity  introduces  another  source  of  variability  which  may  serve  to  increase  a,  the 
probability  of  type  I  error.  In  the  past,  another  approach  to  estimating  LOD  based 
on  the  calibration  equation  has  been  advocated  in  order  to  compensate  for  this  extra 
uncertainty  [24-27].  In  this  approach,  the  (1  -  a/2)  prediction  interval  for  the 
intercept  (i.e.,  the  response  of  the  blank)  from  the  calibration  curve  is  used.  The 
upper  Umit  of  this  interval  corresponds  to  X^  and  the  corresponding  analyte 
concentration  is  the  LOD.  The  procedure  is  illustrated  in  figure  4.  ^ ''- 

The  prediction  interval  for  the  intercept  can  be  found  by  using  eqn.  2.13  with 
Xi  =  0.  The  procedure  is  equivalent  to  using  a  value  for  the  confidence  factor,  k,  in 
eqn.  2.15  calculated  as  follows: 
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1       J''''' 


^-  Wiv-2    1  -  ^  -  ^i  [2.16] 


■^xr 


The  similarity  between  this  equation  and  eqn.  2,8  can  easily  be  seen;  the  2"**  and  3"^ 
terms  in  the  parenthesis  account  for  the  influence  of  the  calibration  conditions  on  the 
value  of  a.  With  this  procedure  for  calculating  the  LOD,  the  k  term  (and  hence  the 
LOD)  will  depend  on  the  calibration  conditions  such  as  the  number  and  range  of  the 
concentrations  of  standards  used  in  calibration,  and  the  use  of  weighted  or 
unweighted  regression  to  estimate  the  prediction  interval.  Using  the  confidence 
factor  in  eqn.  2.16  compensates  for  possible  change  in  a  due  to  these  calibration 
conditions. 

Limit  of  Guaranteed  Detection 

lUPAC  defines  the  LOD  as  "the  minimum  concentration  or  quantity 
detectable"  and  that  it  is  "derived  from  the  smallest  measure  that  can  be  detected 
with  reasonable  certainty  [i.e.,  the  value  of  XJ"  [1].  This  definition  of  LOD  is 
deceptive  since  it  can  lead  to  the  false  assumption  that  if  the  analyte  is  present  at  or 
above  the  LOD  value  that  it  will  always  be  detected  --  i.e.,  result  in  a  measurement 
above  the  signal  detection  limit.  Conversely,  if  a  given  (unknown)  sample  does  not 
give  a  detectable  signal,  then  it  might  be  falsely  assumed  that  the  analyte  must  be 
present  at  a  concentration  less  than  the  LOD. 

By  the  definition  of  LOD  in  eqn.  2.15,  it  is  apparent  that  the  mean  response 
A»Yi  for  the  analyte  present  at  a  concentration  equal  to  the  LOD  is  the  signal 
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detection  limit,  Xj.  If  the  distribution  of  possible  signal  values  is  symmetrical  about 

the  mean,  this  means  that  in  50%  of  the  measurements  where  an  analyte  is  present 
in  a  sample  at  a  concentration  equal  to  the  LOD,  the  resulting  signal  will  not  be 
detected.  The  probabiHty  that  the  signal  due  to  the  analyte  present  at  a  given 
concentration  does  not  give  rise  to  a  detectable  signal  is  known  as  the  probability  of 
type  n  error,  fi,  or  the  probability  of  a  false  negative.  Thus,  when  the  analyte  is 
present  at  a  concentration  equal  to  the  LOD,  B  =  0.5. 

The  lack  of  a  detectable  signal  does  not  mean  that  the  analyte  concentration 
level  must  be  below  the  LOD  value.  It  would  be  useful  to  know  at  what 
concentration  level  the  analyte  must  be  present  in  order  to  be  detected  with  near 
certainty  (i.e.,  with  very  low  fi).  The  inadequate  nature  of  the  LOD  figure  of  merit 
in  this  regard  has  been  noted  by  several  authors  [3,  4,  14,  16].  It  is  possible  to  define 
a  guaranteed  signal  detection  limit,  Xg,  such  that  for  the  distribution  of  signal 
measurement  values  X^,  with  a  mean  equal  to  X  , 

P{X,<Xj)   -  p  [2.17] 

where  the  value  of  fl  is  chosen  according  to  a  predefined  risk  of  type  n  error.  The 
corresponding  limit  in  the  concentration  domain  is  the  limit  of  guaranteed  detection, 
LOG,  which  is  defined  as  follows: 

LOG  =   ^y  "  ^'  -  ^^  [2.18] 
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where  the  second  part  of  the  equation  can  be  used  to  calculate  LOG  if  the  standard 
deviation  of  the  signal  is  the  same  as  that  of  the  blank;  in  such  a  situation,  a  =  fi. 
Figure  5  shows  the  distributions  of  measurements  for  the  blank  and  analyte 
concentrations  equal  to  the  LOD  and  LOG  values  with  k=3.  As  can  be  seen,  the 
LOG  is  a  useful  FOM  since,  if  a  given  sample  is  not  detected  above  Xj,  the  analyst 
can  confidently  state  to  the  customer  (with  only  B  probability  of  error)  that  the 
analyte  is  not  present  at  a  concentration  at  or  above  LOG. 

Limit  of  Quantitation 

One  final  FOM  should  be  mentioned  which  is  related  to  the  limits  of 
detection  and  guaranteed  detection:  the  limit  of  quantitation  (LOQ),  sometimes 
called  the  limit  of  precision  or  the  Umit  of  determination  [3,  4,  16].  Although  the 
two  limits,  LOD  and  LOG,  are  important  for  the  process  of  analyte  detection,  the 
analyst  is  frequently  most  interested  in  quantitation.  It  is  obvious  that  the  precision 
in  quantitation  is  frequently  degraded  near  the  detection  limit  since  the  signal  and 
the  noise  approach  the  same  magnitudes.  The  quantitation  limit,  X^,  in  the  signal 
domain,  is  defined  as 


X,  -  [1,  +  V«  P.19] 

where 

a    =  true  standard  deviation  on  the  analyte  signal,  and 
1/kq  =  desired  relative  standard  deviation  (RSD). 
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Finding  the  corresponding  value  LOQ  in  the  concentration  domain  is  straightforward. 
The  meaning  of  this  hmit  is  the  LOQ  is  the  lowest  concentration  of  analyte  which 
can  be  determined  at  a  pre-defined  level  of  precision  (RSD).  K  it  is  assumed  that 
a  is  constant  and  10%  RSD  is  required,  then 

LOQ  -  i^^  [2.20] 

The  signal  probability  distribution  (with  10%  RSD)  of  analyte  present  at  the  LOQ 
is  also  shown  in  fig.  5. 

Summary 

The  theory  behind  three  related  and  useful  figures  of  merit,  the  Umit  of 
detection,  the  limit  of  guaranteed  detection,  and  the  limit  of  quantitation,  have  been 
reviewed  in  this  section.  The  purpose  of  the  first  two  FOMs  is  to  give  some 
indication  of  the  analytical  procedure's  ability  to  detect  small  amounts  of  analyte. 
Although  the  LOG  is  a  more  useful  value  in  this  regard,  the  LOD  has  been  far  more 
widely  reported  for  various  analytical  procedures.  From  the  standpoint  of 
comparison  of  techniques'  detection  power,  either  FOM  can  be  used  as  long  as 
consistent  definitions  are  used,  all  the  relevant  experimental  detail  is  given  (the 
measurement  time  and  the  electronic  bandwidth  are  often  ignored)  and  an 
appropriate  experimental  protocol  is  used  to  estimate  values  for  a^  and  ag.  The 
procedure  used  to  obtain  these  last  two  values  should  be  given  lq  any  report  of  LOD 
or  LOG  as  well. 
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CHAPTERS 
THE  DETECTION  OF  INDIVIDUAL  ATOMS  OR  MOLECULES 


Laser  Spectroscopic  Methods  of  Analysis 

Lasers  have  become  a  powerful  and  versatile  tool  in  the  arsenal  of  the 
physical  and  analytical  spectroscopist.  Lasers  possess  a  number  of  unique  and 
valuable  properties  such  as  high  degrees  of  directionality,  intensity,  coherence, 
monochromaticity,  and  polarization;  these  qualities  have  opened  up  a  realm  of 
experiments  previously  impossible  using  only  conventional  light  sources.  In 
particular,  in  the  field  of  ultra-trace  analysis,  the  intensity  and  monochromaticity  of 
the  laser  allow  for  the  unique  blend  of  very  high  sensitivity  and  selectivity,  especially 
in  the  field  of  atomic  spectroscopy. 

The  impact  of  lasers  in  spectroscopy  can  perhaps  be  appreciated  by  reviewing 
briefly  some  of  the  possible  results  of  the  interaction  of  an  electromagnetic  field  with 
an  atom  or  a  molecule  in  the  ground  state,  and  some  of  the  related  methods  which 
use  these  processes  as  basis  for  analysis.  These  processes  are  shown  in  fig,  6. 
Interaction  of  the  ground  state  of  the  analyte  with  radiation  at  a  specific  frequency 
results  in  the  production  of  the  excited-state  species  A*  at  a  rate  that  is  proportional 
to  the  spectral  energy  density  of  the  radiation.  Once  in  the  excited  state,  some  of  the 
processes  which  can  be  detected  for  analysis  are  the  production  of  radiation,  charged  I 
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species,  or  heat  dissipation  into  the  surrounding  medium.  Since  the  number  density 

of  A*  (n^*)  is  directly  proportional  to  the  ground-state  population  present  before 
irradiation,  monitoring  the  events  shown  in  the  figure  can  give  information  related 
to  the  concentration  of  analyte  initially  present  in  the  ground-state  level.  The 
processes  and  some  related  analytical  techniques  shown  in  the  figure  can  be 
described  as  follows: 

1.  Stimulated  absorption  of  incident  radiation:  atomic  or  molecular  absorption 
spectroscopy,  in  which  the  attenuation  of  light  irradiating  the  sample  is 
monitored. 

2.  Stimulated  emission  from  A*:  atomic  and  molecular  stimulated  emission 
spectroscopy.  Methods  based  on  this  process  have  not  been  exploited, 
although  techniques  based  on  this  method  have  been  used  to  monitor  flame 
species  [28,  29]. 

3.  Spontaneous  emission  from  A*:  atomic  and  molecular  fluorescence.  Of 
course,  atomic  emission  spectroscopy  is  also  based  on  this  process,  although 
the  excitation  is  not  provided  by  a  hght  source. 

4.  Collisional  deactivation  of  A':  photothermal  spectroscopy,  where  the 
collisional  heating  of  the  surrounding  medium  is  monitored. 

5.  Radiation-induced  ionization:  photoionization  spectroscopy,  in  which  the 
incident  radiation  produces  the  ion-electron  pairs. 


A**e- 
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Figure  6.  Excitation   of  analyte   by  interaction   with  electromagnetic 

radiation,  and  some  of  the  processes  that  can  occur  as  a 
resuh. 
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6.      Collisionally  induced  ionization:  optogalvanic  spectroscopy,  in  which  the 

enhance  rate  of  production  of  analyte  ion/electron  pairs  with  irradiation  is 
monitored  in  an  external  voltage  field. 

For  a  given  volume  of  irradiated  analyte,  lasers  can  result  in  spectral  energy 
densities  that  are  roughly  4-10  orders  of  magnitude  greater  than  most  conventional 
sources.  In  the  processes  hsted  above,  using  the  laser  as  the  light  source  has  resulted 
in  a  great  enhancement  in  sensitivity  and  selectivity.  Indeed,  analytical  techniques 
based  on  processes  2  and  4-6  are  not  practical  without  using  lasers.  Note  that 
scattering  processes  and  related  techniques  such  as  Raman  spectroscopy,  which  have 
also  become  important  analytical  methods  since  the  use  of  laser  sources,  are  not 
hsted  in  the  figure. 

Single  Atom/Molecule  Detection 

Some  of  the  analytical  techniques  outlined  above  are  so  sensitive  and  selective 
that  it  is  possible  to  detect  analyte-specific  events  when  only  a  few  atoms  or 
molecules  interact  with  the  laser.  Indeed,  laser-based  methods  capable  of  detecting 
single  atoms  or  molecules  were  first  reported  over  a  decade  ago  [30-32].  Interpreting 
these  methods  in  terms  of  the  concepts  of  detection  limits,  as  reviewed  in  chapter  2, 
has  not  proven  to  be  straightforward.  Laser-based  techniques  in  which  the  sensitivity 
is  high  enough  (and  the  noise  is  low  enough)  that  individual  species  can  be  detected 
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shall  be  called  single-atom  detection  (SAD)  methods.^    Strict  requirements  for 

techniques  to  be  termed  SAD  will  be  outlined  in  chapter  5;  all  others  which  possess 
detection  limits  of  only  a  few  atoms  can  be  called  near-SAD  methods. 

Past  SAD  Experiments:  A  Sampling  of  Applications  and  Techniques 

Why  be  concerned  with  the  ability  to  detect  single  atoms?  Certainly  it  would 
seem  that  the  vast  majority  of  analytical  methods,  even  those  based  on  laser-induced 
processes  for  ultra-trace  analysis,  would  fall  far  short  of  requiring  the  capability  of 
SAD.  However,  there  are  several  beneficial  reasons  to  being  aware  of  the  particular 
concepts  which  apply  to  (near-)SAD  techniques;  these  are  now  briefly  outlined. 

First  of  all,  there  are  applications  in  physics  and  chemistry,  including  analytical 
chemistry,  which  do  in  fact  require  a  capability  of  SAD.  Even  if  the  goal  is  not 
quantitation  of  analyte,  the  capability  of  a  laser-based  method  to  detect  single  atoms 
often  must  be  evaluated,  and  many  of  the  concepts  of  SAD  theory  apply. 

Current  ultra-trace  laser  spectroscopic  methods  are  getting  ever  closer  to  the 
SAD  regime  of  analysis.  The  LODs  reported  in  conventional  bulk  analysis  often 
translate  into  very  few  atoms  in  the  laser  during  the  measurement  time.  In  addition, 
for  some  of  the  methods,  the  selectivity  is  so  great  that  it  is  possible  to  reduce  the 
noise  to  almost  nonexistent  levels;  many  of  the  concepts  from  SAD  theory  are 
relevant  in  these  cases. 


The  term  "SAD"  will  be  used  even  though  the  detected  species  may  be  atoms, 
molecules,  ions  or  radicals.  The  term  "atom"  in  reference  to  SAD  techniques  will 
apply  to  all  such  species  as  well. 
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Finally,  the  ultimate  goal  of  an  analytical  method  is  the  detection  of  single 

atoms.  The  ability  to  quantitate  the  amount  of  analyte  in  a  sample  at  the  atomic 
level  certainly  represents  this  ultimate  goal  in  chemical  analysis,  and  many  vmique 
and  interesting  experiments  in  analysis  and  other  fields  would  no  doubt  result  from 
such  an  abiHty. 

Applications  of  SAD 

This  is  by  no  means  meant  to  be  an  exhaustive  list  of  possible  applications  of 
SAD.  There  have  been  reports  in  the  literature  of  possible  applications  of  SAD 
techniques  [33-36];  this  hst  covers  some  of  these  as  well  as  various  others.  The 
applications  Hsted  here  specifically  call  for  the  abiHty  to  detect  the  signal  due  to 
single  atoms. 
Physical  applications 

1.  It  is  possible  to  observe  and  measure  the  transport,  diffusion  or  otherwise,  of 
individual  atoms  through  the  volume  defined  by  the  laser  beam,  as  well  as 
other  statistical  mechanics  applications  in  which  the  fluctuations  in  various 
atomic/molecular  processes  must  be  monitored.  The  observation  of  the 
gas-phase  reactions  of  individual  atoms  or  molecules  [37]  might  be  included 
in  this  type  of  application. 

2.  The  (mean)  Ufetimes  of  various  excited  states  of  individual  atoms  or  • 
molecules  can  be  reported  and  studied  in  various  environments.     This 
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information  is  usually  only  available  as  an  ensemble  average  over  many 

atoms/molecules. 

3.  The  spectroscopic  features  of  free  atoms  or  ions  [37,  38]  can  be  studied;  the 
spectroscopy  of  molecules  in  various  sites  in  a  solid  matrix  can  also  be 
investigated  [39-42]. 

4.  Sorting  of  individual  species  (atoms,  ions,  molecules)  requires  the  ability  to 
"tag"  and  detect  them  with  the  laser.  This  is  an  example  of  Maxwell's  sorting 
demon  [35,  43]. 

5.  The  detection  of  rare  events,  such  as  solar  neutrinos,  by  their  effect  on 
individual  atoms  [33,  34,  44]. 

6.  The  observation  of  the  orientation  of  individual  species  in  clusters.  SAD  has 
been  applied  to  clusters  of  ions  in  a  laser-cooled  ion  trap  [45]. 

Analytical  applications 

1.  Detection  of  very  low  bulk  concentration  of  analyte,  where  the  limiting  "noise" 
on  the  signal  is  due  to  the  statistical  appearance  of  individual  atoms  in  the 
laser. 

2.  Detection  of  very  rare  isotopes  and  application  in  various  related  fields  such 
as  geochemistry  and  environmental  analysis  [43, 46]  and  cosmochemistry  [35]. 

3.  Application  of  SAD  methods  in  surface  analysis  is  necessary  for  even 
relatively  "high"  bulk  and  surface  concentrations  since  very  few  analyte  atoms 
will  be  analyzed  at  any  time  by  combinations  of  sputtering  methods  with 
laser-based  detection  [47]. 
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4,      Detection  of  anaiyte  in  difficult  matrices.  Simple  dilution  is  often  an  effective 

method  of  correction  for  matrix  effects,  but  is  limited  by  the  sensitivity  of  the 
analytical  technique.  With  the  development  of  methods  capable  of  SAD  in 
large  samples,  the  sample  can  be  diluted  to  an  arbitrarily  high  degree. 

Techniques  Capable  of  SAD 

Techniques  which  have  achieved  (near-)SAD  in  the  past  can  be  broadly 
categorized  as  using  either  destructive  or  nondestructive  methods  of  detection.  In 
the  former  class,  the  atom  is  consumed  during  the  detection  process,  producing  at 
most  one  single  "count";  in  the  second  category,  each  atom  can  produce  multiple 
detectable  events  during  its  interaction  with  the  laser.  In  this  work,  most  attention 
will  be  focused  on  two  laser-based  methods  which  have  had  the  most  success  at 
detecting  single  atoms:  resonance  ionization  spectroscopy  (RIS)  and  related  methods, 
and  laser-induced  fluorescence  (LIF).  These  two  methods  can  represent  the  two 
classes  of  SAD  methods:  RIS  is  almost  always  destructive,  while  with  LIF  it  is 
frequently  possible  for  each  atom  to  emit  many  photons  during  its  interaction  with 
the  laser.  ~     * 

Resonance  ionization  spectroscopy 

Resonance  ionization  spectroscopy  involves  the  laser-assisted  production  and 
subsequent  detection  of  anaiyte  ions.  Production  of  ion-electron  pairs  using  lasers 
can  approach  100%  efficiencies.  Since  the  manipulation  of  charged  particles  is 
relatively  easy  and  detection  of  these  particles  can  also  be  highly  efficient,  it  is  not 
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surprising  that  there  have  been  a  number  of  (near-)SAD  reports  using  RIS  [31,  32, 

43,  44,  48-54]  as  well  as  a  number  of  review  articles  which  cover  various  theoretical 
and  practical  aspects  of  the  technique  [35,  36,  55,  56]. 

There  are  two  basic  ways  in  which  the  ion  has  been  produced  through 
laser-induced  processes;  these  are  shown  in  figure  7.  The  method  in  fig.  7(a),  direct 
photoionization  by  the  laser,  has  been  developed  at  Oak  Ridge  National  Laboratories 
as  outlined  by  Hurst  and  Payne  [35];  the  detection  of  the  ion  in  a  buffer  gas  has 
usually  been  achieved  by  a  proportional  counter.  The  process  shown  in  fig.  7(b) 
involves  excitation  to  a  long-lived  Rydberg  level  in  a  vacuum  with  subsequent  field 
ionization,  and  detection  by  a  secondary  electron  multipUer  (with  or  without  a  mass 
spectrometer).  The  advantages  of  excitation  to  a  Rydberg  level  is  that  less  powerful 
lasers  are  necessary  for  saturation,  which  may  result  in  a  more  general  technique  and 
less  laser-induced  background.  The  development  of  this  method  by  Letokhov  and 
Bekov  has  been  outlined  in  a  recent  book  [36]. 
Laser-induced  fluorescence 

Techniques  of  involving  the  detection  of  spontaneous  emission  of  laser-excited 
analyte  atoms  have  been  characterized  by  extremely  high  selectivity  and  sensitivity. 
This  high  sensitivity  has  been  accompanied  by  a  number  of  reports  of  claiming  SAD 
[30,  37,  38,  45,  57-69].  There  are  a  large  number  of  possible  combinations  of 
excitation/detection  schemes  using  LIF;  figure  8  shows  three  of  them  which  represent 
the  types  of  schemes  which  can  be  used. 
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Figure  7.  Two  possible  ionization  methods  for  RIS. 

(a)  photoionization  from  an  intermediate  level. 

(b)  field  ionization  from  a  Rydberg  level. 
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The  technique  of  LIF  can  be  either  destructive  or  nondestructive,  depending 

on  the  particular  analyte,  the  environment  and  the  time-scale  of  the  interaction  with 
the  laser.  The  presence  of  a  metastable  level  can  act  as  an  effective  "trap"  for  atoms 
during  the  measurement;  for  example,  if  a  scheme  similar  to  the  one  shown  in  fig. 
8(b)  is  used  when  the  middle  level  has  a  very  long  lifetime  and  the  analysis  takes 
place  in  a  vacuum,  then  LIF  will  be  a  destructive  technique.  Similarly, 
photodestruction  of  molecules  by  the  laser  may  occur  and  effectively  limit  the 
emitted  photons  to  a  small  number.  On  the  other  hand,  when  there  is  good 
coUisional  and  optical  coupling  between  all  the  levels  involved  in  the 
excitation/detection  scheme,  then  "cycling"  of  the  atom  back  to  the  ground-state  is 
possible,  where  the  atom  can  then  further  interact  with  the  laser. 
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CHAPTER  4 
THE  VARIABIUTY  OF  ESTIMATED  LIMITS  OF  DETECTION 


The  value  of  the  LOD  calculated  for  a  given  analytical  procedure  can  vary 
due  to  a  number  of  factors.  In  this  chapter,  these  factors  will  be  investigated,  and 
the  variation  in  calculated  LOD  due  to  random  fluctuations  will  be  theoretically 
evaluated.  Before  proceeding  further,  however,  a  brief  introduction  to  estimation 
theory  is  necessary  to  fully  appreciate  the  concepts  involved  in  calculating  an  LOD 
for  a  given  analytical  procedure. 

Estimation  Theory 

Variables  are  defined  in  terms  of  the  population  parameters  of  their 
probabihty  distribution,  such  as  the  mean,  n^,  and  variance,  a/,  of  the  variable  x. 
However,  the  true  values  of  these  parameters  are  very  rarely  known  exactly  and  must 
be  estimated  from  a  sample  of  the  population.  Functions  of  the  sample  which 
provide  such  estimates  are  known  as  sample  statistics;  common  sample  statistics 
include  the  sample  mean  x  and  variance  s^^^.  The  value  of  a  single  sample  statistic 
is  also  known  as  a  point  estimate  of  the  corresponding  population  parameter. 

The  value  of  a  point  estimate  depends  upon  the  particular  sample  chosen; 
thus,  point  estimators  are  themselves  variables,  with  properties  dependent  on  the  size 
of  the  sample  chosen  and  the  estimation  function  used.  The  theory  of  the  behavior 
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of  sample  statistics  as  variable  estimators  of  the  population  parameters  is  known  as 

estimation  theory.  Besides  the  two  already  given,  some  other  common  estimators  - 
include  least-square  estimators  (seen  in  chapter  2)  for  the  calibration  parameters  of 
slope,  tto,  and  intercept,  Mf  Desirable  characteristics  of  sample  statistics  include  the  ^ 
property  of  being  unbiased  estimators,  when  the  mean  of  the  point  estimates  is  the 
equal  to  the  population  parameter,  and  the  property  of  efficiency,  which  is  indicated 
by  a  low  variance  of  the  estimate.  The  standard  deviation  of  the  statistic  is  also 
known  as  the  standard  error.  ,,  '  », 

If  information  is  known  about  the  standard  deviation  of  the  estimate,  then  it 
is  possible  to  construct  an  interval  about  the  point  estimate  within  which  the  true 
value  of  the  population  parameter  is  likely  to  lie.  Such  an  interval  is  known  as  the 
confidence  interval.  The  probability  that  the  parameter  Ues  outside  the  interval,  a, 
is  inversely  proportional  to  the  size  of  the  interval.  Recall  that  confidence  intervals 
were  briefly  mentioned  earlier  in  relation  to  estimating  the  true  analytical  response 
from  a  calibration  curve. 

The  Limit  of  Detection  as  a  Population  Parameter 

The  LOD  was  defined  in  equation  2.15  according  to 

LOD  =  ^  [4.1] 
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where 

k  =  the  confidence  factor, 

a^  =  the  standard  deviation  in  blank  measurement,  and 

oeq  =  the  analytical  sensitivity. 
In  hght  of  the  previous  discussion  on  estimation,  it  can  be  easily  understood  that 
since  the  LOD  is  defined  in  terms  of  population  parameters,  the  LOD  is  itself  a 
parameter  of  the  analytical  system,  which  can  be  estimated  according  to 

LOD  -  ^  [4.2] 

where 

LOD  =  calculated  estimate  of  true  LOD  value. 

Thus,  regardless  of  how  the  confidence  factor  k  is  chosen  in  the  above 

A, 

equations,  the  value  of  LOD  has  a  given  variance  associated  with  it.  Since  the  true 
value  of  the  LOD  is  frequently  of  great  importance  as  a  performance  characteristic 
of  a  given  analytical  method,  it  is  of  interest  to  estimate  this  variance,  ct^qd- 

A 

Estimating  the  variabiUty  of  a  given  LOD  would  serve  several  purposes:  (1) 
confidence  intervals  within  which  the  true  LOD  value  would  he  could  be  calculated; 

A, 

(2)  hypothesis  tests  with  different  LOD  values  can  be  made  in  order  to  determine  if 
there  is  a  significant  difference;  and  (3)  the  influence  of  various  experimental 
procedures  on  CTlod  can  be  assessed.  In  many  research  laboratories,  for  example, 
calibration  curves  may  only  be  used  to  provide  an  estimate  of  the  sensitivity,  Uq,  to 
use  in  calculating  LOD.  Time  saving  methods  of  providing  this  sensitivity  estimate 


>-7<r^5 


40 

which  do  not  adversely  affect  the  quality  of  the  LOD  estimate  (ie  increase  CTlod) 

would  be  welcome. 

Variability  of  LOD 

To  derive  an  equation  for  aLOD»  ^  propagation  of  errors  approach  can  be  used 
with  eqn.  4.2: 


'LOD 
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[4.3] 


To  solve  this  equation,  it  is  necessary  to  estimate  the  variance  of  s^.  The  variance 
of  this  estimate  can  be  related  to  the  variance  of  the  estimator  for  a^,  again  by 
propagation  of  errors: 


The  estimated  variance,  s\,  for  n  measurements  is  distributed  as  follows: 


[4.4] 
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where  the  x^  distribution  has  a  mean  of  (n-1)  and  a  variance  of  2(n-l)  [70].  From 
this,  we  can  deduce 
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Now  we  can  estimate  the  variance  of  the  background  standard  deviation  as 


si 


C2 

'"       2(n-l) 


[4.7] 


Substituting  this  into  eqn.  4.3  gives  the  following  estimate  for  aLoo- 
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where 


s^  =  estimate  of  the  standard  error  of  the  slope  (sensitivity),  and 
Slod  =  estimate  for  the  standard  error  of  LOD. 
Rearrangement  allows  for  two  useful  forms  of  the  above  equation: 
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where 

CVjj  =  the  coefficient  of  variation  of  variable  (x  =  ajy.^,  and 

n  =  number  of  blank  measurements  to  calculate  s^. 
Both  of  the  above  equations  are  valid  no  matter  how  the  confidence  coefficient  is 
chosen  in  calculating  the  LOD  by  eqn.  4.2 
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Confidence  Limits  and  Comparing  Values  of  LOP 

The  validity  of  the  propagation  of  errors  approach  used  to  derive  an  equation 
depends  on  the  following  two  conditions:  (1)  the  errors  in  the  slope  and  background 
variances  must  be  independent;  and  (2)  the  term  CV3  must  be  small  (below  about 
0.10).  Effects  of  violation  of  the  first  condition  (i.e.,  including  blank  values  in  the 
calibration  curve)  are  probably  only  slight  when  a  reasonable  number  of  points  are 
included  in  the  calibration  set.  The  second  restriction  arises  due  to  the  non-linear 
relationship  between  the  LOD  and  the  slope;  for  high  coefficients  of  variation  of  the 
latter,  the  propagation  of  errors  approach  breaks  down  [71];  effects  of  high  CV, 
values  will  be  investigated  in  chapter  7. 

Keeping  the  above  two  restrictions  in  mind,  the  utility  of  eqns.  4.9  and  4.10 
are  as  follows.  Equation  4.10  gives  an  indication  of  the  expected  fluctuation  in  LOD 
estimate  due  only  to  the  variable  nature  of  point  estimates;  this  value  can  be  used 
to  construct  confidence  intervals  within  which  the  true  LOD  is  Ukely  to  he. 
Examination  of  eqn.  4.9  shows  how  the  variabihty  of  LOD  can  be  divided  between 
the  uncertainty  in  the  estimates  of  s^  and  s^.  For  a  given  set  of  conditions  used  to 
calculate  LOD,  eqn.  4.9  can  be  used  to  determine  the  relative  contribution  of  the  two 

A 

error  sources.   The  method  used  to  determine  the  confidence  interval  about  LOD 
depends  on  the  relative  magnitude  of  the  two  terms.    If  the  first  term  dominates, 
then  most  of  the  variation  in  the  estimate  is  due  to  the  uncertainty  in  estimating  a^. 
This  situation  is  probably  the  most  common  in  analytical  chemistry.  In  these  cases, 

A, 

the  X"  distribution  can  be  used  to  construct  confidence  intervals  on  the  LOD  value 
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in  the  same  manner  as  for  the  sample  standard  deviation  [72];  Kaiser  has 

demonstrated  similar  calculations  [3].  For  example,  for  a  two-sided  95%  interval  and 
20  measurements  of  the  blank, 

0.16(LOD)    i  LOD  i  1.46  (LdD)  [4.11] 

Confidence  intervals  for  different  numbers  of  blank  measurements  can  easily  be 
constructed  by  using  the  appropriate  degrees  of  freedom.  For  comparison  of  LOD 
values  obtained  under  similar  conditions,  F-tests  can  be  used  [72]. 

As  the  second  term  in  eqn.  4.9  becomes  more  important,  the  distribution  of 

A, 

LOD  will  come  to  approximate  a  normal  distribution;  in  these  cases,  standard 
t-tables  can  be  used  to  construct  confidence  intervals  of  the  type  LOD  ±  kSLoo.  and 
to  compare  LOD  estimates  with  the  t-test.  In  either  case,  eqn.  4.10  is  valid  and  gives 

A. 

a  good  idea  of  the  variability  of  a  given  value  of  LOD  due  to  random  fluctuation  of 
the  sample  statistics. 

Summary 

Inspection  of  eqns.  4.1  and  4.2  leads  to  the  conclusion  that  there  are  two 
sources  of  variation  in  calculated  values  of  LOD:  the  first  source  is  due  to  actual 
changes  in  population  parameters  a^  and  a^  of  the  analytical  technique  -  a  shift  in 
alignment,  slightly  increased  background  noise,  the  presence  of  interferences, 
improved  analytical  methodology  -  and  the  second  source  is  due  to  the  variability 
of  the  point  estimates  used  in  eqn.  4.2.  The  source  and  magnitude  of  the  latter 
fluctuation  in  LOD  can  be  seen  by  appUcation  of  eqns.  4.9  and  4.10.  One  benefit  of 
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viewing  the  calculated  value  of  the  LOD  as  an  estimate  of  the  true  value  of  the  LOD 
is  that  these  two  sources  of  fluctuation  of  LOD  can  be  separated;  if  the  LOD 
fluctuates  by  an  amount  much  greater  than  the  calculated  value  of  Slqd  indicates, 
then  the  probable  source  of  the  change  in  the  observed  LOD  is  an  actual  change  in 
the  analytical  parameters  of  the  system.  The  performance  and  use  of  eqns.  4.9  and 
4.10  will  be  further  investigated  in  chapter  7. 


CHAPTERS 
THEORY  OF  THE  DETECTION  AND  COUNTING  OF  ATOMS 


Introduction 

Some  of  the  practical  aspects  of  detecting  individual  atoms  or  molecules  with 
lasers  were  outlined  in  chapter  3.  As  stated  in  that  chapter,  certain  laser-based 
methods  have  achieved  high  enough  sensitivity  so  that  single  atom  detection  may  be 
possible.  In  the  Uterature  there  are  reports  of  methods  which  were  able  to  detect 
single  atoms  with  very  high  S/N  ratio  [30,  37,  44].  However,  most  of  the  reports  of 
SAD  have  had  a  comparatively  low  S/N  ratio  for  each  atom;  in  these  cases,  there  is 
a  question  of  whether  or  not  single  atoms  can  actually  be  detected  above  the 
background  (if  any  background  is  present).  Evaluation  of  these  methods  often 
involves  the  following  questions: 

1.  Can  individual  atoms  can  be  detected? 

2.  Under  what  conditions  is  single  atom  detection  possible? 

3.  If  it  is  possible  to  detect  single  atoms,  is  it  also  possible  to  count  the  numbers 
of  atoms  passing  through  the  laser  beam? 

4.  What  are  the  characteristics  of  these  (possible  SAD)  methods  when  utilized 
in  a  more  conventional  sense,  ie  to  measure  bulk  concentration  of  analyte? 
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There  axe  problems  in  applying  the  theory  of  detection  limits  summarized  in 

chapter  2  in  attempting  to  answer  these  questions.  In  comparison  with  practical 
aspects,  theoretical  considerations  for  (possible)  SAD  methods  have  received  very 
little  attention  in  the  literature.  There  have  been  a  handful  of  papers  which  have 
evaluated  laser  spectroscopic  techniques  with  respect  to  the  potential  to  detect  atoms 
[73-78];  these  include  some  theoretical  discussion  of  SAD  methods.  Recently,  several 
papers  have  attempted  to  deal  with  the  problem  of  verifying  that  single  molecules 
were  being  detected  by  LIF  as  they  flow  through  the  laser  [65-68].  However,  the 
work  of  Alkemade  remains  the  only  in-depth,  systematic  general  theoretical 
treatment  of  SAD  to  date;  this  work  is  presented  in  two  classic  papers  [79,  80]  and 
has  recently  been  reviewed  and  extended  [81]. 

The  purpose  of  this  chapter  is  to  answer  the  questions  posed  earUer  which 
may  arise  for  methods  which  produce  a  signal  for  individual  atoms  which  might  be 
detectable  above  the  background  noise.  The  concepts  from  chapter  2  are  applied  in 
a  logical  manner  to  SAD  methods,  and  strict  definitions  and  important  figures  of 
merit  for  (near-)SAD  methods  are  presented.  Several  other  important  factors  in 
evaluating  possible  SAD  methods  will  also  be  discussed.  The  treatment  of  SAD  in 
this  chapter  owes  much  to  Alkemade's  original  treatment  of  the  subject;  many  of  the 
terms  used  are  identical,  although  the  meanings  may  have  been  modified.  The  intent 
of  this  treatment  is  to  generalize  Alkemade's  pioneering  work  and  to  present  a 
general  theory  of  SAD  in  a  form  useful  in  the  development  of  methods  which  may 
achieve  the  goal  of  detecting  and  counting  atoms.   Most  of  the  concepts  presented 


47 

are  verified  and  further  illustrated  in  chapter  8  through  the  use  of  computer  models 

of  SAD  experiments. 

Definition  of  an  SAD  Method 

The  definition  of  an  SAD  method  is  a  generalization  of  Alkemade's  four 
criteria  for  true  SAD  [80]  and  is  the  following: 

A  method  is  an  SAD  method  if  each  and  every  atom  which  interacts 
with  the  laser  can  be  detected  above  the  background  noise. 

The  above  definition  will  be  re-stated  in  a  more  rigorous  form  later  in  this 
chapter;  nevertheless,  the  general  concept  of  an  SAD  method  can  be  appreciated. 
Implicit  in  the  above  statement  is  that  we  are  concerned  only  with  the  atoms  which 
actually  interact  with  the  laser.  Equally  important  is  to  notice  the  difference  between 
a  method  in  which  some  (but  not  all!)  individual  atoms  can  be  detected  above  the 
background  and  an  SAD  method:  it  is  not  enough  to  simply  have  a  certain  likelihood 
of  detecting  single  atoms,  but  every  atom  which  is  probed  by  the  laser  must  be 
detected. 
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General  Model  of  SAD  Methods 

The  Poisson  Process 

Many  of  the  processes  involved  in  SAD  methods  v^dll  be  assumed  to  be  related 
to  the  Poisson  distribution.  These  include  the  number  of  atoms  probed  by  the  laser, 
the  detection  probability  of  an  atom,  and  the  number  of  detected  events  per  atom. 
To  fully  understand  the  nature  and  limitations  of  SAD  model  presented  in  this 
chapter,  a  firm  grasp  of  the  properties  of  Poisson  variables  is  necessary  [82,  83], 

Experiments  which  measure  the  (variable)  number  of  discrete  occurrences  of 
an  event  in  a  certain  length  of  time,  or  in  a  given  area  of  space,  frequently  deal  with 
variables  possessing  a  Poisson  distribution.  For  a  Poisson  variable,  the  probability 
of  X  number  of  events  occurring  during  a  given  fixed  interval  time  or  space,  t,  is 
given  by: 

PiX)   .  m±e:l  [5.1] 


where 

0  =  flux  of  events  per  unit  time/space,  and 

A  variable  which  is  truly  a  Poisson  variable  (characteristic  of  a  Poisson 
process)  possesses  the  following  quahties: 

1.       The  probability  of  an  event  occurring  vidthin  the  interval  of  time  or  space  is 
small  compared  to  the  probability  that  it  can  occur  elsewhere. 
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2.       The  probability  of  an  individual  event  occurring  within  the  length  of  time  or 

space  is  independent  of  all  other  events  which  have  occurred,  either  during 

that  length  of  time/space,  or  outside  it.  / 

The  Poisson  distribution  is  closely  related  to  several  other  distributions, 

including  the  uniform  random  distribution,  and  the  Gamma  and  exponential 

distributions.  This  relationship  can  be  understood  by  considering  a  Poisson  variable 

in  time.    When  a  given  event  can  occur  randomly  in  time  (i.e.,  the  probability 

distribution  of  the  time  of  occurrence  is  a  uniform  distribution),  then  the  number  of 

events  during  a  given  time  interval  follows  a  Poisson  distribution.  This  is  the  basis 

of  a  typical  Poisson  process.  For  such  a  process,  the  probability  distribution  of  the 

time  interval,  t,  between  events  is  given  by  an  exponential  distribution 

Pit)   -  4)  e-*'  [5-2] 

where  the  meaning  of  0  is  the  same  as  in  eqn.  5.1.  The  mean  time  interval  between 
events  is  given  by  0'\  The  exponential  distribution  is  a  simplified  form  of  the 
Gamma  distribution.  The  variable  in  the  Gamma  distribution  is  the  amount  of  time, 
t^,  for  a  specified  number,  v,  of  events  to  occur: 

Pit)   -      »^^"^"*^' 

;  [5.3] 
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Typical  SAD  Experiment 

The  general  form  of  a  laser-based  SAD  method  can  be  illustrated  with  the  aid 
of  figure  9.  In  the  figure,  analyte  atoms  flow  past  a  region  of  interaction  with  a  laser 
beam;  atoms  which  interact  with  the  laser  produce  a  number  of  detectable  events. 
The  method  of  detection  can  be  either  destructive  (e.g.,  ion  detection  in  an  atomic 
beam)  or  nondestructive  (e.g.,  fluorescence  detection  of  molecules  in  a  flowing 
stream).  The  number  of  events  detected  during  a  measurement  time  T„  are  counted. 
During  this  time,  a  certain  volume  of  sample  containing  analyte,  V^,  flows  past  the 
laser  beam;  for  the  analyte  atoms  within  this  volume  there  is  a  probability  of  entering 
a  region  in  which  it  is  possible  to  interact  with  the  laser  beam  and  emit  detectable 
events.  This  volume  is  the  probe  volume,  Vp,  and  is  defined  by  the  region  of 
intersection  of  the  flowing  stream  with  the  laser  beam  which  can  be  viewed  by  the 
detector.  A  certain  number  Np  of  atoms  which  enter  Vp  interact  with  the  laser  beam 
during  T^;  a  given  atom  interacts  with  the  laser  beam  for  time  t^.  The  atom's 
interaction  time  is  a  function  of  various  factors  such  as  the  magnitude  of  T„,  the  type 
of  laser  used  (probed,  continuous  or  modulated),  and  the  atom's  residence  time,  t^, 
within  Vp.  The  value  of  t^  is  usually  a  variable,  with  mean  t^  depending  on  such 
factors  as  velocity,  diffusion,  and  size  and  shape  of  Vp.  The  atom's  corresponding 
interaction  time  t;  may  also  be  a  variable,  with  mean  Tj,  depending  on  the  relative 
magnitude  of  T^  and  t^,  and  the  conditions  of  the  experiment. 
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Detection  efficiency:  general  definition 

For  the  general  layout  of  a  typical  SAD  experiment  as  given  above,  it  is 
convenient  to  define  the  detection  efficiency  of  atoms  which  enter  the  probe  volume. 
The  detection  efficiency,  e^,  is  defined  as  the  probability  that  any  given  atom,  during 
its  interaction  with  the  laser,  produces  a  signal  that  can  be  distinguished  from  the 
background  which  arises  during  the  measurement  time.  Alkemade  gives  a  similar 
definition  for  the  detection  efficiency  [80],  but  there  are  two  important  differences: 
(1)  Alkemade  was  only  concerned  with  the  case  where  there  was  no  background 
noise;  hence,  e^  was  simply  the  probability  that  an  atom  which  interacted  with  the 
laser  produced  at  least  one  detectable  event;  (2)  Alkemade's  term  only  appUed 
during  a  single  probing  time  (e.g.,  one  single  laser  pulse  in  a  pulsed  experiment) 
rather  than  during  an  atom's  entire  interaction  time  with  the  laser.  One 
characteristic  which  the  two  definitions  have  in  common  is  that  they  are  only 
concerned  with  atoms  which  actually  interact  with  the  laser  (i.e.,  atoms  for  which 
ti  >  0).  .      . 

Signal  production 

This  section  will  introduce  terms  relating  to  the  signal  detected  during  a  single 
measurement;  assumptions  which  will  be  made  regarding  the  probabiUty  distribution 
of  these  signals  in  the  SAD  model  will  be  discussed  in  the  next  section.  During  the 
measurement  time  T„j,  there  may  be  a  certain  (variable)  number  of  background 
counts,  lb,  with  a  mean  given  by  '  '        ,'■■">        -'^ .  *-'■  *  '*^'« 
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\^t  -<i>6i;  ^•'!  ^    -.^^  H  [5.4] 
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where  '  "  '  ,  . 

0^  =  mean  flux  of  background  noise  (count/s)         , 

/ij,  =  mean  noise  during  T„  (counts). 
The  total  signal,  \,  recorded  during  T„  is  due  to  the  contribution  of  noise  and 
analyte  signal.  For  nondestructive  detection,  this  can  be  given  as 

^t-h^E  is  [5.5] 

where 

ij  =  number  of  detected  events  due  to  each  individual  atom  which  flows 
through  V  .  This  variable  has  a  mean  given  by 

^  -  *,  r,  [5.6] 

where 

05  =  mean  flux  of  signal  from  a  given  atom  (count  s'^  atom'^). 
The  mean  total  signal  for  a  destructive  technique  is  given  by 


/,  -  Ht  +  BjN 


[5.7] 


where  N    is  the  number  of  atoms  which  were  probed  during  T„  and  e^  is  the 
detection  efficiency,  the  probability  that  a  given  atom  will  be  detected. 

The  physical  meaning  of  the  term  03  depends  on  the  method  used.  If  a 
nondestructive  method  such  as  cyclic  LIF  is  used,  the  term  refers  to  the  actual  rate 
of  detected  events  from  single  atoms.  For  destructive  methods,  however,  the 
meaning  is  obviously  different  since  there  can  be  no  more  than  one  detected  event 
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per  atom.  In  this  case,  the  term  is  actually  the  reciprocal  of  the  mean  detection  time 
of  the  atom  within  the  laser.  In  other  words,  for  a  given  atom  which  mteracts  with 
the  laser  within  Vp  the  mean  time  until  a  detected  event  is  produced  is  0/ . 
Theoretical  expressions  for  0^  for  a  number  of  cases  for  LIF  and  RIS  can  be  found 
in  the  literature  [78,  81,  84,  85]. 

Figure  9  has  only  shown  one  type  of  possible  SAD  experiment,  specifically  for 
a  case  where  atoms  continuously  flow  through  Vp.  Figure  10  depicts  several 
alternative  examples  of  SAD  methods.  As  can  be  seen,  various  possible  relative 
magnitudes  of  T„  and  t,  are  possible:  e.g.,  when  a  continuous-flow  atomizer  is  used 
with  a  pulsed  dye  laser,  if  the  detected  events  are  counted  for  each  pulse,  then  it  is 
frequently  true  that  t^  >  >  T„,  and  the  atoms  are  "frozen"  during  the  measurement. 
This  situation  was  termed  the  stationary  case  by  Alkemade  [80].  The  nonstationary 
case  occurs  when  the  measurement  time  T„  is  of  the  same  magnitude  or  larger  than 
the  typical  value  of  t^  For  example,  in  a  heated  cell,  the  atoms  or  molecules  may 
be  free  to  diffuse  in  and  out  of  Vp  during  T„;  a  case  such  as  shown  in  fig.  9  with  a 
continuous  wave  (CW)  laser  would  be  nonstationary. 

Basic  Assumptions  for  the  SAD  Model 

Number  of  probed  atoms.  Np 

In  a  solution  or  sample  which  contains  analyte  atoms,  it  is  often  reasonable 
to  assume  that  the  analyte  atoms  are  randomly  distributed  throughout  the  sample. 
Such  being  the  case,  it  can  be  assumed  that  the  appearance  of  analyte  atoms  in  Vp 
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Figure  10.  Examples  of  possible  SAD  methods.  Typical  examples  of 
nonstationary  methods  are  methods  (a)  and  (b);  methods  such 
as  (c)  and  (d)  are  often  under  stationary  conditions  during  T, 
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is  a  Poisson  process,  and  the  number  N^,  is  a  Poisson  variable  governed  by  eqn.  5.1 

with  a  mean  which  is  dependent  on  Vp,  T„,  and  the  analyte  concentration.    This 
assumption  will  abnost  always  be  valid  in  the  absence  of  severe  clustering  effects  or 
very  small  probe  volumes  (ie  when  analyte  atoms  cannot  be  treated  as  infinitely 
small);  these  cases  violate  the  requirements  of  a  Poisson  process. 
Number  of  detected  events.  I^ 

The  distribution  of  the  background  counts,  I^,  will  be  assumed  to  follow  a 
Poisson  distribution.  In  other  words,  the  limiting  noise  will  be  background  shot 
noise.  Even  though  this  assumption  is  made  for  convenience  (and  applies  in  many 
counting  situations),  it  is  quite  possible  for  the  background  noise  SAD  experiment 
to  be  flicker  noise  [69];  in  such  cases,  the  concepts  in  this  and  other  chapters  still 
apply,  but  with  some  slight  modification. 

The  detection  process  in  SAD  methods  will  also  be  assumed  to  be  a  Poisson 
process.  For  nondestructive  methods,  such  as  cyclic  LIF,  this  means  that  the  number 
of  counts  per  atom,  i„  detected  during  T^,  is  a  Poisson  variable  in  time  such  as  was 
described  in  an  earlier  section;  all  of  the  detected  photons  are  randomly  distributed 
during  the  atom's  interaction  time  with  the  laser.  Since  lb  is  also  assumed  to  follow 
a  Poisson  distribution  with  a  mean  of  Mb  =  <PiJm>  then  when  a  single  atom  passes 
through  Vp  during  T„,  I^  has  a  Poisson  probability  distribution  with  a  mean  of  ^t  + 
0sti- 

For  a  destructive  detection  method,  such  as  RIS,  the  detection  process  is  also 
a  Poisson  process,  although  in  a  more  subtle  way.  A  large  number  of  atoms  N  all 
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interacting  with  the  laser  beam  at  one  time  would  result  in  a  mean  flux  0^  of 

detected  ions;  the  detection  times  of  these  ions  are  assumed  to  be  uniform  random 
variables,  and  the  time  intervals  between  detection  follow  an  exponential  distribution 
such  as  in  eqn.  5.2,  with  mean  time  between  detection  of  0/^.  K  this  assumption  is 
true,  then  the  probabiHty  that  at  least  one  single  ion  will  be  detected  by  time  t;  can 
be  found  by  integrating  the  exponential  distribution  from  t  =  0  to  t  =  ti: 

PiO^t^ti)   -  1  -«-♦'''  [5.8] 

Since  this  is  assumed  to  be  a  Poisson  process,  the  production  and  detection  of  ions 
are  independent  events;  thus,  the  above  equation  applies  even  if  there  is  only  one 
atom  irradiated  by  the  laser.  However,  the  meaning  of  the  <p^  term  has  changed, 
since  a  "flux"  of  ions  is  of  course  not  possible  with  only  one  possible  ion.  As 
explained  earlier,  <p^  is  now  considered  to  be  the  reciprocal  mean  time  to  detection 
of  the  ion  produced  from  the  single  atom.  As  with  any  Poisson  process,  the 
distribution  of  detection  times  is  an  exponential  distribution.  Note  that  eqn.  5.8 
applies  also  for  nondestructive  detection,  which  is  also  assumed  to  be  a  Poisson 
process;  in  this  case,  the  variable  "detection  time"  is  the  time  before  one  event  due 
to  a  single  atom  is  detected. 

The  mean  signal  due  to  N  atoms  irradiated  during  T^  was  given  in  equation 
5.7  as  e<jN  .  We  now  consider  the  distribution  of  the  signal  produced:  the  number 
of  detected  ions,  Nj,  when  a  fixed  number  of  atoms,  N  ,  are  probed  by  the  laser  is 
given  by  the  binomial  distribution: 
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PiN,)   -[^';](e^)^'(l-e^)V'  [5.9] 


However,  as  explained  above,  Np  is  assumed  to  be  a  Poisson  variable;  hence,  the 
probability  distribution  of  Nj  will  also  follow  a  Poisson  distribution,  with  mean  e^N  . 
Obviously,  when  only  one  atom  crosses  through  Vp  during  T„,  the  signal  due  to  the 
analyte  is  either  zero  (not  detected)  or  one  (detected). 
Variability  of  I^ 

As  described  in  the  previous  section,  the  detection  process  for  the  signal 
produced  by  atoms  in  the  laser  beam  is  considered  to  be  a  Poisson  process.  In  this 
section,  the  variabihty  in  the  total  signal,  represented  by  the  value  of  the  total 
variance,  a^^,  will  be  investigated  for  cyclic  LIF. 

Equation  5.5  gives  the  total  signal  due  to  Np  atoms  interacting  with  the  laser 
beam  during  a  measurement  time  T„.  Since  the  background  is  assumed  to  possess 
a  Poisson  probability  distribution,  and  the  number  of  detected  photons  from  each 
individual  atom  also  follows  a  Poisson  distribution  with  a  mean  given  by  eqn.  5.6,  it 
would  seem  that  I^  should  also  be  a  Poisson  distribution  with  mean  and  variance 
given  by 

T,-a]-ii,^  (j)/,.  [5.10] 

However,  eqn.  5.10  is  only  valid  when  both  0,  and  t;  are  constant  for  every 
atom  which  can  interact  with  the  laser.  Although  this  can  be  true,  depending  on  the 
experimental  conditions,  it  can  easily  be  the  case  that  both  0^  and  tj  are  variables, 
with  mean  i^  and  Tj,  respectively.  The  values  of  <p^  and  tj  for  a  particular  atom  may 
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depend  on  the  path  of  that  atom  through  Vp.   For  example,  0^  depends  upon  the 

optical  collection  efficiency,  and  this  may  not  be  constant  over  the  entire  probe 
volume;  another  source  of  variation  in  0^  with  path  is  if  the  transition  of  the  atom 
is  not  saturated,  and  the  laser  intensity  is  not  constant  throughout  V  .  The 
interaction  time,  tj,  of  an  atom  will  not  be  constant  if  diffusion  effects  play  a 
significant  role  during  T„,  or  if  the  shape  of  Vp  is  such  that  (for  example)  an  atom 
travelUng  down  the  center  of  Vp  will  have  a  larger  interaction  time  than  one  which 
skirts  the  edge. 

Thus,  although  the  number  of  detected  photons  from  a  given  atom  which 
interacts  with  the  laser  will  follow  a  Poisson  distribution  with  a  mean  given  by  eqn. 
5.6,  the  overall  distribution  of  photoelectrons  due  to  a  single  atom,  I^  will  have  a 
mean  given  by 


[5.11] 


and  Is  will  not  follow  a  Poisson  distribution.  The  mean  of  the  total  signal  can  thus 
be  written  for  the  general  LIF  case: 


It-^t,T^-N^<b,x, 


[5.12] 


The  effect  of  having  either  03  or  tj  variable  is  to  increase  the  variance  of  I^. 
The  variance  of  I^  (for  a  fixed  value  of  Np)  can  be  partitioned  between  the  variance 
due  to  the  Poisson  detection  process  in  both  the  background  and  signal  counts  ("shot 
noise")  and  the  "extra"  variance  due  to  any  variabihty  in  0^  and  X{. 


U... 
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a]  -T,^   {N^)^oHl,)  [5.13] 

where  the  second  term  in  the  equation  is  due  to  the  "extra"  variance.  Of  course, 
when  05  and  tj  can  both  be  assumed  to  be  reasonably  constant,  then  this  term 
approaches  zero  and  \  will  be  approximately  Poisson. 

Note  that  this  section  only  discussed  the  case  of  LIF.  Of  course,  similar 
considerations  are  involved  in  any  technique  based  on  destructive  detection.  The 
case  of  variable  interaction  time  and  LIF  detection  will  be  addressed  in  chapter  8. 

Signal  Detection  Limit  for  the  SAD  Model 

With  the  SAD  model  and  definition  as  given  above,  the  application  of 
detection  limit  theory  is  as  follows.  For  a  given  measurement  time  T„,  the 
distribution  of  \  is  Poisson  with  mean  and  variance  Mb  •  From  eqn.  2.3  ,  the  signal 
detection  limit  X^  is  set  according  to  a  pre-defined  tolerance  (denoted  by  a,  the 
probability  for  type  I  error)  for  false  positives: 

P{I,  ^Xj)   -  a  [5.14] 

A  value  for  a  must  be  chosen  before  the  experiment;  the  value  of  X^  is  set  so  that 
the  observed  probability  of  one  or  more  false  counts  during  T^  (due  to  background 
noise)  is  at  or  below  this  level.  Note  that,  since  I^  is  a  discrete  variable,  the  value 
of  a  will  not  be  uniformly  decreased  by  increasing  the  value  of  X^.  K I^  is  a  Poisson 
variable  then  estimating  the  mean  Pt  during  T„  and  application  of  eqn.  5.1  will  allow 
Xj  to  be  set  through  the  use  of  tables  of  the  Poisson  probability  distribution  [86]. 
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K  it  is  found  that,  for  given  values  of  T„  and  a, 

P(/i>0)   ^  a  [5.15] 

then  there  is  essentially  no  background  noise  during  T^;  this  is  called  the 
intrinsic-noise  limit,  since  the  only  noise  on  the  analyte  signal  is  due  to  variance  of 
the  signal  itself.  At  the  intrinsic  hmit,  a  single  count  (or  more)  indicates  the 
presence  of  analyte;  the  value  of  X^  is  one  count. 

Detection  Efficiency  of  a  near-SAP  Method 

The  general  definition  for  the  detection  efficiency,  e^,  has  been  given 
previously  as  that  probability  of  a  given  atom  will  result  in  a  signal  detectable  above 
the  background  noise.  Now  we  can  state  more  clearly  that  the  detection  efficiency 
is  defined  such  that,  when  a  single  atom  interacts  with  the  laser  during  T^, 

e^  -  P(/,  i  X^)  [5.16] 

where  X^  is  chosen  according  to  a  pre-defined  tolerance  for  false  positives.  At  the 
intrinsic  limit,  the  detection  efficiency  is  limited  by  the  noise  inherent  on  the  signal 
itself.  Since  we  have  assumed  a  Poisson  detection  process  for  both  destructive  and 
nondestructive  detection,  we  can  see  from  eqn.  5.8  that,  for  ^^  =  0,  the  probabiUty 
that  an  atom  entering  Vp  will  be  detected  is  given  by: 

e^  -  1  -  exp(-(J)/,.)  [5.17] 
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This  equation  is  equivalent  to  the  one  used  by  Alkemade  [80],  who  was  only 

concerned  with  the  intrinsic-limited  case.   When  noise  is  indeed  present,  however, 
eqn.  5.16  must  be  used. 

Requirements  for  SAD 
General  Requirement 

The  general  definition  of  an  SAD  method,  given  earlier,  is  that  the  method 
detects  each  and  every  atom  which  interacts  with  the  laser  with  near-certainty.  This 
requirement  is  now  given  in  a  more  succinct  form:  an  SAD  method  is  a  method  in 
which 

e^  ^  1  -  P  [5.18] 

where  the  detection  efficiency  is  defined  in  eqn.  5.16  and  fi  is  the  probability  of  type 
n  error  (false  negative)  as  described  in  chapter  2.  The  highest  allowable  value  of  B 
must  be  decided  prior  to  the  evaluation  of  the  (possible)  SAD  method. 

The  application  of  the  above  general  requirement  for  SAD  is  different  in  the 
cases  of  RIS  and  cycUc  LIF.  The  difference  between  the  two  as  SAD  methods  can 
best  be  understood  by  studying  figure  11  for  fixed  Np  and  equivalent  detection 
efficiencies.  The  intrinsic-limited  case  is  shown  in  the  figure,  and  e^  can  be 
calculated  by  using  eqn.  5.17.  It  is  assumed  that  0^  and  tj  are  constant  for  all  N 
atoms  in  both  cases  (such  a  situation  is  reasonable  for  an  atomic  beam  experiment 
with  a  pulsed  dye  laser  and  a  small  Vp). 
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Achieving  SAD  with  RIS  (Destructive  Detection't 

With  a  destructive  method  such  as  RIS,  a  single  atom  can  only  give  rise  to  one 
coimt  at  the  most.  The  detection  efficiency  in  this  case  is  simply  the  binomial 
probabiUty  of  "success"  (eqn.  5.9)  --  for  RIS,  the  probability  that  the  given  atom  will 
be  ionized  during  its  interaction  time  tj.  This  probability  is  given  in  eqn.  5.8  for  a 
Poisson  process.  However,  the  only  way  in  which  the  requirement  for  SAD  will  be 
met  as  set  forth  in  eqn.  5.18  is  if  X^  =  1;  i.e.,  the  intrinsic-limited  case.  Thus,  for 
true  SAD  using  RIS  (or  any  other  destructive  technique)  the  following  two  conditions 
must  hold  (from  eqns.  5.15,  5.17  and  5.18): 

P(/t>0)    ^  a 

Note  that  the  second  condition  assumes  constant  0^  and  tj.  The  effect  of  variable 
values  of  these  parameters  on  the  overall  detection  efficiency  must  be  taken  into 
account  if  necessary. 

Achieving  SAD  with  LIF  (Nondestructive  Detection^) 

Guaranteed  detection  limit  (X^V  Recall  that  the  concept  of  a  guaranteed 
signal  detection  Hmit,  Xg,  was  introduced  in  chapter  2.  The  value  of  Xg  is  helpful  in 
clarifying  the  requirements  for  a  nondestructive  technique  to  be  a  true  SAD  method. 
For  a  given  value  of  X^,  the  value  Xg  is  defined  so  that  the  probability  of  a  variable 
in  a  distribution  with  mean  Xg  being  less  than  X^  is  negligible  (less  than  a  desired 


"g 


65 

probability  B  of  type  11  error).  K  both  the  background  and  the  signal  are  described 

by  Poisson  probability  distributions,  it  is  easy  enough  to  assign  values  to  X^  and  X 
for  any  value  of  Pt  by  using  tables  of  Poisson  values  [86].  Table  1  shows  these  values 
for  a  number  of  cases,  and  different  values  of  a  and  B.  The  procedure  in 
determining  these  signal  detection  limits  is  as  follows:  from  the  value  of  A»b .  X^  is 
chosen  so  that  F(\  >  X^)  «  a.  From  this  value  of  X^,  a  Poisson  distribution  is  found 
such  that  P(X  >  X^,)  «  1-B.  The  mean  of  this  distribution  is  Xg. 

If  it  is  assumed  that  both  <p^  and  t;  are  constant,  then  if  Xg  is  found  from 
Poisson  tables  as  described,  the  requirement  for  SAD  by  LIF  is 

h  ^  X^-<i>,T„  [5.19] 

Figure  12  shows  a  situation  with  n^,  =  1  count  and  SAD  is  possible  by  LIF  detection 
(a  =  B  =  0.0014).  When  /i^  =  1  count  during  T„,  X,,  =  6  counts  and  Xg  =  16 
counts;  thus,  by  eqn.  15,  SAD  is  possible  with  i^  >  15  counts/atom.  Note  from  table 
1  that  even  in  the  intrinsic-Hmited  case  (^^  =  0),  a  value  of  i^  =  6.6  counts/atom  is 
needed  for  SAD  by  LIF  (at  the  99.86%  confidence  level).  •'       '  '       <  ■• '^ 

'  .'  '■■■''  l* 

Summary:  Requirements  for  SAD  .  .  ,, 

The  basic  requirement  for  SAD  by  any  laser-based  method  is  given  by  eqn. 
5.18.  The  practical  consequences  of  this  requirement  for  both  destructive  and 
nondestructive  cases  have  also  been  discussed  in  this  section.  True  SAD  is  possible 
with  destructive  detection  only  at  the  intrinsic-limit;  however,  SAD  is  possible  by 
nondestructive  means  even  in  the  presence  of  noise  if  the  sensitivity  is  high  enough. 


Table  1 
The  Two  Limits  for  an  SAD  Experiment 
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Mean  blank  level  (p^,) 

Detection  limit  (X^) 
a  <0.0014 
(a  <  0.05) 

Guaranteed  limit  (X  ) 
(3  «  0.0014 
0  «  0.05) 

0.00 

1 
(1) 

6.6 
(3.3)                  1 

0.05 

2 
(1) 

8.9 
(3) 

0.25 

4 
(2) 

12.6 
(4.7) 

1.00 

6 
(4) 

16 
(7.8) 

5.00 

14 
(10) 

28 
(16) 

10.00 

22 
(16) 

39 
(23) 

100.00 

132  [130]^ 
117  [117f 

[169]^ 

[136]^                 1 

The  limits  are  given  for  the  signal  domain,  with  all  the  signals  given  as  counts. 

^The  values  in  the  square  brackets  were  found  by  assuming  a  Gaussian  distribution 
with  the  appropriate  k  values.  At  these  higher  signal  levels,  the  Poisson  distribution 
can  be  approximated  by  a  Gaussian  with  \i  =  cr. 
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It  should  be  emphasized  that  it  is  possible  for  a  given  technique  to  detect  individual 

atoms  and  still  not  be  a  true  SAD  method;  for  example,  with  Mb  =  0  the  requirement 
for  SAD  (at  the  99.86%  confidence  level)  with  LIF  is  6.6  counts/atom  (assuming  0^, 
and  tj  constant);  however,  if  i^  =  1  count/atom,  then  individual  atoms  would  still  be 
detected  quite  often  (e^  =  0.632).  This  is  an  example  of  a  near-SAD  method,  where 
detection  of  an  atom  in  the  laser  beam  during  T^  is  possible  (and  perhaps  likely)  but 
not  certain. 

Detection  Efficiencv  as  a  FOM.  Notice  that  even  though  the  detection  limit 
theory  of  chapter  2  was  applied  to  the  SAD  case  in  the  signal  domain,  it  is  somewhat 
difficult  to  speak  of  near-SAD  and  SAD  methods  in  terms  of  LOD.^  Intuitively,  it 
is  not  possible  to  have  LOD  <  1  atom,  or  as  a  non-integral  number  of  atoms.  This 
would  seem  to  indicate  that  LOD  (and  LOG)  is  not  an  ideal  FOM  for  the  evaluation 
and  comparison  of  (near-)  SAD  methods,  as  it  is  for  more  conventional  methods. 
Another  possible  FOM  is  to  compare  S/N  for  one  (or  more)  atoms,  where  S/N  is 
the  ratio  of  the  mean  signal  due  to  one  atom  to  the  background  noise.  A  more 
informative  FOM  than  S/N  for  near-SAD  methods  is  the  detection  efficiency,  e^. 
This  parameter  contains  information  about  the  magnitudes  of  both  the  mean  signal 
and  the  background  noise,  as  well  as  the  noise  on  the  signal  due  to  single  atoms. 
Improvements  in  near-SAD  methods  should  be  evaluated  by  improvements  in  the 
value  of  e^  rather  than  an  increase  in  S/N  (which  may  be  at  a  cost  to  e^  if  there  is 


Vith  LOD  in  terms  of  numbers  of  atoms  in  the  laser  beam,  not  bulk  concentration 
(this  aspect  will  be  discussed  later). 
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an  increased  variance  in  the  signal  due  to  one  atom).  Once  SAD  has  been  achieved 

(e^  «  1),  then  improvements  in  the  SAD  method  would  best  be  described  in  terms 
of  increases  in  S/N  due  to  single  atoms. 

Precision  of  Counting  Atoms 

Thus  far  in  this  chapter  we  have  treated  the  likelihood  of  detecting  single 
atoms.  However,  the  question  remains  as  to  whether  it  is  possible  to  precisely  count 
the  number  of  atoms  which  pass  through  the  laser  if  Np  >  1  atom.  Recall  that  in 
chapter  2,  a  FOM  called  the  limit  of  quantitation  (LOQ)  was  introduced  which 
specified  an  analyte  concentration  above  which  it  was  possible  to  determine  an 
unknown  sample  with  a  predefined  degree  of  precision  (usually  with  RSD  =  0.10, 
following  the  suggestion  of  Kaiser  [3]).  For  an  SAD  method,  it  is  possible  to  define 
a  similar  FOM  in  terms  of  atoms  which  pass  through  the  laser  beam  during  T„;  in 
this  section,  a  related  FOM  shall  be  investigated:  the  minimum  sensitivity,  (i^  \, 
necessary  to  count  atoms  with  a  pre-defined  precision  at  all  levels  of  N  . 

Precision  of  Signal 

In  the  measurement  of  a  fixed  Np  in  the  SAD  model,^  the  signal  distributions 
in  the  LIF  and  RIS  experiments  are  given  by  the  Poisson  and  binomial  distributions. 


It  IS  assumed  that  (p^  and  tj  are  constant  throughout  the  section  on  precision  of 
counting  atoms. 


m 

and  by  substituting  the  appropriate  values  for  a  from  these  distributions,  the 
following  equations  are  obtained  for  the  intrinsic-limited  case: 


RIS:  RSD 


I  1-^d 


\\^A 


p . 


UF:  RSD  - 


[5.20(a)] 


[5.20(b)] 


From  eqn.  5.20(a),  it  can  be  calculated  that  for  a  precision  of  10%  or  better  with  RIS 
it  is  necessary  that  e^  >  0.99.  Since  the  RSD  improves  as  Np  increases  for  RIS,  an 
RIS  method  with  this  detection  efficiency  or  better  is  capable  of  counting  atoms  for 
any  value  of  N  . 

The  situation  for  LIF  is  different,  however,  since  many  events  can  be  detected 
from  a  single  atom.  Substituting  Np  =  1  atom  in  eqn,  5.20(b)  results  in  a 
requirement  of  100  photoelectrons/atom  for  RSD  =  0.10  with  Np  =  1.  Thus,  it 
would  seem  that  for  precise  counting  of  atoms  with  LIF,  a  sensitivity  of  at  least  100 
photoelectrons/atom  during  tj  is  required  at  the  intrinsic  limit.  This  is  a  far  more 
stringent  requirement  than  the  6.6  photoelectrons/atom  which  are  necessary  for  SAD 
(at  the  99.86%  level). 

There  is  a  problem,  however,  when  the  RSD  is  calculated  using  the  above 
formula  that  stems  from  the  difference  between  the  precision  of  signal  measurement, 
RSD^,  and  the  precision  of  counting  atoms,  RSD^.,  when  using  a  nondestructive  SAD 
method.    This  problem  does  not  arise  with  RIS  since,  when  SAD  is  possible,  it  is 


■^;^.' 
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essentially  true  that  every  atom  gives  rise  to  a  single  count;  thus,  the  signal  exactly 

follows  the  number  of  atoms. 

Consider  the  situation  for  LIF  with  100  photoelectrons/atom,  shown  in  figure 
13.  This  figure  illustrates  the  difference  between  the  signal  precision  and  the 
counting  precision.  Although,  by  eqn.  5.20(b),  this  situation  corresponds  to 
RSD^  =  0.1  for  Np  =  1  atom,  it  is  obvious  that  there  is  very  little  possibiHty  of 
incorrectly  counting  atoms  when  Np  =  1  since  there  is  almost  no  overlap  between  the 
distributions.  Obviously,  the  value  for  RSD„  is  not  a  reflection  of  the  counting 
precision  of  the  LIF  method. 

Counting  Precision 

This  section  is  concerned  with  nondestructive  detection  only,  since  the 
requirements  for  precise  counting  by  RIS  are  essentially  the  same  as  the 
requirements  for  SAD.  For  a  nondestructive  technique,  the  value  of  N  is  estimated 
according  to  the  following  equation: 


^s 


where  the  RND  function  rounds  the  expression  in  the  parentheses  to  the  nearest 
whole  number,  and  the  integer  N„  is  the  number  of  measured  atoms  (i.e.,  the 
estimate  for  Np).   We  can  define  the  counting  precision,  RSD^,  as 

RSD,  =  -^-^  [5.22] 
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where  o{NJ  is  the  standard  deviation  in  the  number  of  measured  atoms  with  fixed 

The  difference  between  RSD^  and  RSD^  can  be  seen  in  figure  14,  which 
shows  the  signal  probability  distribution  with  Mb  =  0  counts,  i^  =  20  counts/atom  and 
N  =  5  atoms.  The  top  axis  displays  the  values  of  N^  which  would  be  calculated  at 
a  given  signal  value  using  eqn.  5.21.  The  dashed  lines  in  the  figure  show  the  portions 
of  the  probability  distribution  which  would  result  in  values  of  4,  5,  or  6  atoms  for  N„. 
Notice  that  N„  is  an  integer;  thus,  a  signal  of  75  photoelectrons,  for  example,  would 
result  in  a  value  of  N^  =  4  atoms  (and  not  3.75  atoms!)  as  the  best  estimate  for  the 
true  value  of  Np. 
Calculation  of  RSD^ 

In  order  to  investigate  the  counting  precision  of  an  SAD  method,  and 
calculate  a  value  of  (i^  )c  for  a  given  background,  the  theoretical  value  of  RSD^  for 
various  fixed  N  must  be  determined.  The  theoretical  value  of  RSD^  is  not  as  easily 
calculated  as  for  RSD„  in  eqn.  5.20.  To  do  so,  an  expression  for  a(N„)  in  eqn.  5.22 
must  be  found  for  given  conditions  of  Np,  sensitivity  and  noise.  From  the  definition 
of  the  variance  of  a  discrete  variable  [87],  it  is  known  that 

oHNJ   -J:^N„-N^)'  P{NJ  [5.23] 

where  P(N^)  is  the  probability  distribution  of  N„.  This  probability  distribution  can 
be  written  as 
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Pif^m)   -  E^(A)  t5.24] 


where  the  summation  limits,  X,  and  X„  are  found  as  follows. 


For  N„  =  0: 


X,  -  0 
X„  -  X^-1 


For  N„  =  1: 


X  -  X 


For  N„  >  1: 


Xt  -  INT[{N„-0.5)<i,^i  +  jiJ+1 
X„  -  /iVr[(iV„^-0.5)<l)/i  +  \x^] 

In  all  cases  the  INT  function  represents  the  integral  part  of  the  expression  in  the 
parentheses.  ....  i    r 

The  above  expressions  are  tedious  to  solve  manually;  a  computer  program  can 
be  written  to  evaluate  these  expressions  for  the  SAD  model  for  various  values  of  Np, 
/lb  and  0s  to  determine  (i^  )c  if  the  distribution  of  It  is  knovra.  The  results,  and 
comparison  of  the  theoretical  value  with  the  observed  value  of  RSD^  from  computer 
simulations,  are  presented  in  chapter  8.  For  the  intrinsic  limit  the  above  equations 
give  a  value  of  (i^  )c  =  35  counts/atom  for  RSD^  <  0.1  for  all  values  of  Np. 

In  any  practical  situation,  the  value  of  RSD^  would  be  almost  impossible  to 
determine  since  N  would  not  be  fixed  but  would  be  a  variable  that  would  change  for 
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different  measurements  (according  to  a  Poisson  distribution  for  the  model  presented 

in  this  chapter).  Nevertheless,  the  purpose  behind  this  section  is  to  illustrate  that  for 
a  nondestructive  technique  such  as  LIF,  the  requirement  of  SAD  is  not  sufficient  to 
count  atoms  during  T^  to  an  arbitrary  precision  (unlike  the  case  with  a  destructive 
SAD  method).  There  is  an  additional  increase  in  sensitivity  required  before  this  is 
possible  by  LIF. 

Scope  of  an  SAD  Method 

Effect  of  the  Measurement  Time  on  an  SAD  Method 

The  requirements  for  an  SAD  method  have  been  presented;  as  an  illustration 
of  what  an  "SAD  method"  means  in  a  more  conventional  analytical  situation,  let  us 
assume  that  somehow  the  flow  of  atoms  can  be  directed  so  that  every  atom  in  a 
given  sample  interacts  with  the  laser.  If  the  technique  is  capable  of  SAD,  does  this 
mean  that  the  method  possesses  infinite  detection  capabiUties?  In  other  words,  can 
any  concentration  of  analyte  in  the  sample  can  be  detected? 

The  general  SAD  model  presented  the  requirements  for  an  SAD  method  for 
a  given  measurement  time,  T„.  The  meaning  and  limitation  of  this  requirement 
should  be  very  clear:  a  method  which  is  truly  capable  of  SAD  is  capable  of  detecting 
above  the  background  noise  every  atom  which  crosses  the  laser  during  T^.  Thus,  in 
partial  answer  to  the  above  questions,  it  is  not  possible  to  simply  state  that  a  given 
technique  is  an  "SAD  method"  without  stating  the  conditions  under  which  this  is 
possible  -  most  particularly  the  measurement  time,  T^.  A  simple  numerical  example 
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with  SAD  by  pulsed  LIF  will  illustrate  this  point  and  serve  to  answer  the  above 

questions. 

Numerical  example:  LIF  with  pulsed  lasers 

Let  us  imagine  the  interaction  of  a  beam  with  a  pulsed  dye  laser,  where 
stationary  conditions  apply.^  During  a  single  laser  pulse  it  is  found  that  ^b  =  0.25 
counts;  i.e.,  there  is  one  "noise"  count  every  four  laser  pulses  on  the  average.  During 
a  single  laser  pulse  (T„  is  the  pulse  duration),  we  can  say  for  the  SAD  model 
presented  in  this  chapter  that  X^  =  3  (a  =  0.00216)  and  Xg  =  10.8  (fi  =  0.00143). 
Thus,  from  eqn.  5.19,  we  see  that  if  i^  =  10.55  counts/atom,  then  SAD  is  possible 
during  a  single  laser  pulse. 

Thus,  this  is  a  true  SAD  technique.  Now  suppose  that  a  given  sample  is 
analyzed  and  will  take  10,000  laser  shots  to  completely  flow  past  the  laser.  Again, 
assume  that  every  atom  in  the  sample  interacts  with  the  laser  for  time  tj.  Is  it 
possible  to  reliably  detect  a  single  atom  in  the  sample  with  the  technique  just 
described? 

The  answer  is  no,  because  the  scope  of  the  above  SAD  method  is  only  a  single 
laser  pulse.  When  T„  is  changed  to  10,000  laser  shots,  the  requirements  for  SAD 
will  change.  If  X^  =  3  counts  is  used  as  a  criterion  to  distinguish  the  presence  of 
analyte  atoms  from  the  background  noise,  then  from  the  value  of  a  for  one  laser 


^each  atom  interacts  with  the  laser  for  one  pulse  only  and  tj  is  fixed.  It  is  assumed 
that  0s  is  constant  as  well.  It  will  be  assumed  that  the  number  of  detected  events  (ie 
photoelectrons  above  a  discriminator  level)  can  be  unambiguously  counted.  In 
reality,  this  may  not  be  so  easy  with  typical  pulsed  dye  laser  experiments. 
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pulse,  we  see  that  during  T„  =  10,000  shots  there  will  be  216  false  positives  on  the 
average.  Obviously  it  would  be  impossible  to  detect  a  single  atom  with  this  value  for 

With  the  new  value  of  T„,  we  must  choose  X^  high  enough  so  that  a  is  at  the 
desired  value.  When  /it  =  0-25  counts/pulse  and  the  background  follows  a  Poisson 
distribution  it  can  be  calculated  that 

P{Il,i7)    -  9.734   X  10"^ 

SO  that  for  10,000  shots  with  X^  =  7  counts,  a  =  0.000973,  an  acceptable  level.  This 
value  of  Xj  gives  Xg  =  18  (13  =  0.001043);  the  requirement  to  detect  a  single  atom  in 
the  sample  is  that  i^  =  17.25  counts/atom,  instead  of  10.55  counts/atom. 

This  simple  illustration  shows  that  a  given  SAD  method  has  a  certain 
"scope"  --  i.e.,  a  certain  value  of  T^  over  which  the  technique  can  detect  a  single 
atom.  Beyond  this  measurement  time,  the  method  can  no  longer  detect  single  atoms 
with  the  required  values  of  a  and  fi.  If  the  sensitivity  of  the  method  is  very  high, 
however,  the  scope  of  the  SAD  method  (measurement  time  over  which  SAD  is 
possible)  may  be  so  long  as  to  be  practically  infinite.  In  other  words,  X^  can  be  set 
so  high  that  it  is  extremely  unlikely  that  I^  will  ever  exceed  X^  during  a  single  laser 
pulse,  no  matter  how  many  pulses  are  counted.  If  SAD  is  still  possible  with  such  an 
Xj  value,  then  the  method  is  truly  capable  of  detecting  a  single  analyte  atom  in  any 
(reasonable)  size  sample. 

Counting  precision.  RSD^.  The  value  of  RSD^,  and  the  requirement  for 
precise  counting  of  Np  in  each  laser  pulse  depends  only  on  the  value  Mb  (and  on  the 
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variability  of  I,;  however,  we  assume  0^  and  tj  are  constant).   Thus,  increasing  the 

length  of  T^  does  not  affect  the  value  of  minimum  sensitivity  necessary  to  precisely 

count  the  number  of  atoms  present  in  each  laser  pulse. 

Continuous  Monitoring  of  Atoms 

One  of  the  most  promising  methods  of  achieving  SAD  is  by  using  continuous 
lasers  with  LIF  detection  (CW-LJF).  This  method  has  been  used  in  the  past  to 
detect  single  molecules  [62,  65-68]  and  atoms  [59,  60,  63]  as  they  flow  through  the 
laser.  The  evaluation  of  near-SAD  methods  based  on  CW-LIF  is  a  task  which  must 
be  carefully  approached.  This  section  will  discuss  some  possible  signal-processing 
methods  which  apply  the  general  SAD  theory  discussed  thus  far.  In  chapter  8,  some 
of  the  methods  and  ideas  discussed  in  this  section  will  be  demonstrated. 

Signal  Processing  Methods 

Simple  integration  over  the  measurement  time.  T^ 

Most  analyses  based  on  atomic  or  molecular  fluorescence  of  bulk  analyte  in 
a  sample  solution,  in  the  simplest  case,  will  integrate  the  signal  for  the  measiurement 
time  and  the  sum  (or  its  normalized  analog,  the  average)  will  be  the  measurement 
value.  This  situation  was  depicted  in  fig.  1  and  fig.  2.  Such  an  approach  is  of  course 
possible  with  a  method  based  on  CW-LIF.  Application  of  the  SAD  theory  when 
using  this  "simple  integration"  signal  processing  method  is  straightforward:  the 
detection  efficiency  is  measured  for  a  given  value  of  X^  (chosen  for  the  pre-defined 
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tolerance  to  false  positive  detection)  and  eqn.  5.18  is  used  to  determine  if  SAD  is 

possible  during  T^. 

The  simple  integration  method  is  very  inefficient  when  T^  >  >  Tj  since  the 
mean  signal  due  to  a  single  analyte  atom  will  never  exceed  i^r^  while  the  mean  blank 
value  will  increase  linearly  with  T„.  In  addition,  there  is  no  temporal  information 
on  the  exact  time  the  atom  is  in  the  laser  beam;  it  is  only  known  that  the  atom 
entered  Vp  sometime  during  T^.  Nevertheless,  there  may  be  certain  situations  in 
which  the  simplicity  of  the  method  has  its  advantages.  In  the  analysis  of  discrete 
samples  which  flow  quickly  through  the  laser  beam  (e.g.,in  LIF  of  analyte  atoms 
atomized  in  a  furnace,  or  in  flow  injection  analysis),  then  SAD  may  be  possible  when 
T„  is  chosen  so  that  the  entire  sample  is  analyzed. 

Far  more  efficient  methods  for  the  analysis  of  continuously  flowing  analyte 
solution  through  Vp  are  based  on  the  use  of  a  time  "window,"  of  length  t^  apphed 
repeatedly  during  T„.  Three  such  methods  will  be  presented  here.  Their  common 
feature  is  that  the  only  the  number  of  photoelectron  counts  within  the  window  t^  are 
considered  at  any  one  time,  and  the  duration  of  t^  is  chosen  to  be  of  the  magnitude 
of  Tj.  Thus,  the  S/N  ratio  due  to  single  analyte  atoms  is  increased. 
Sequential  Application  of  Time  Window  (t^   <  t  ) 

It  may  be  that  the  sensitivity  of  the  CW-LIF  method  is  so  high  that  it  is  still 
possible  to  detect  the  presence  of  a  single  atom  even  if  t^  «  t^/ 10.  The  S/N  due  to 
a  single  atom  would  decrease  relative  to  a  situation  in  which  t^  =  t,;  however,  if 
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SAD  is  still  possible  by  eqn.  5.18,  then  this  is  the  best  method  to  use.  Such  situations 
have  only  rarely  been  reported  in  the  Uterature  [30,  37]. 

The  application  of  this  method  is  simply  the  use  of  sequential  integrations  for 
time  t„  for  the  entire  measurement  time  T„.  In  choosing  the  value  of  X^  for  t„  it 
should  be  remembered  that  X^  should  be  high  enough  so  that  the  total  number  of 
false  positives  from  T^/t^  windows  will  be  at  the  desired  tolerance  level.  The 
disadvantage  of  this  method  is  that  higher  sensitivity  is  required  to  achieve  SAD  than 
by  either  of  the  next  two  methods  which  use  a  time  window.  However,  if  this 
sensitivity  is  available,  then  the  method  is  very  simple  to  use  and  can  count  atoms  in 
real  time  as  they  traverse  the  laser  beam.  For  CW-LIF  which  is  on  the  borderUne 
of  becoming  a  true  SAD  method  and  does  not  have  a  high  enough  sensitivity  for  this 
method,  one  of  the  following  two  methods  can  be  used. 
Photon  burst  method 

This  method  has  been  reported  in  the  Uterature  as  a  method  of  recording 
spectra  of  single  atoms  in  an  atomic  beam  [61,  63].  Figure  15  demonstrates  how  the 
"photon  burst"  method  may  be  applied  to  CW-LIF.  A  single  photoelectron  count 
triggers  open  the  gate,  which  is  pre-set  to  count  the  photoelectrons  for  the  gate  time 
t^  The  number  of  photoelectrons  (including  the  trigger  pulse)  constitutes  the  "signal" 
during  t^  The  gate  duration  is  set  so  that  t^  «  t^.  The  photon  burst  method  can  be 
implemented  during  real  time  so  that  bursts  above  a  pre-set  value  (X^)  can  signal  the 
presence  of  an  analyte  atom  in  the  laser  beam.  Alternately,  the  number  of  counts 
from  all  the  photon  biursts  can  be  stored  in  a  computer  and  analyzed  later;  only  the 
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bursts  with  a  sum  equal  to  or  greater  than  X^  can  be  considered  due  to  analyte 

atoms  with  1-a  probability  of  false  positive.  All  the  normal  requirements  given 
previously  for  an  SAD  method  apply.  The  advantage  of  the  photon  burst  method 
over  simple  integration  is  that  only  the  noise  during  t„  can  contribute  to  false 
positives,  and  consequently  a  lower  value  of  X,,  can  be  chosen  and  the  sensitivity 
necessary  for  SAD  is  not  as  high. 

False  positives  during  T^.  As  explained  previously,  as  T^^  increases,  it  is 
frequently  necessary  to  increase  X^  to  keep  the  probability  of  a  false  positive  during 
T„  down  to  an  acceptable  value.  For  the  "simple  integration"  and  "sequential 
window"  methods  above,  it  is  relatively  easy  to  determine  the  value  of  a  for  given 
values  of  X^  and  T^.  The  theoretically  calculated  value  of  a  (assuming  a  Poisson 
distribution  of  I^)  for  the  "photon  burst"  method  is  slightly  more  complicated.  In  the 
presence  of  background  noise  only,  the  sum,  S,  of  the  number  of  photoelectrons  in 
a  given  photon  "burst"  follows  a  Poisson  distribution  based  on  the  mean  background 
level  such  that, 

PiS)   -  PUt-S-l) 

e-'^iVi,)'-^  [5.25] 

(5-1)  ! 

since  S  >  1  count.  Thus,  for  a  certain  background  level,  the  probabiHty  of  a  given 
burst  giving  a  false  positive  is 


PiSiXj)   -  f;  Pil^)  [5.26] 
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so  that  during  the  measurement  time, 

a  -355; />(/,)  [5.27] 

where 

b  =  average  number  of  bursts  during  T„. 
For  a  given  value  of  ^b  and  T^,  the  average  number  of  bursts  can  be  calculated 
according  to 

The  dependence  of  a  on  T„  can  readily  be  seen  from  eqns.  5.27  and  5.28.  However, 
the  denominator  in  eqn.  5.28  ensures  that  a  increases  very  slowly  (compared  to  the 
simple  integration  method)  for  a  given  X^  as  T„  increases. 

The  above  calculations  of  the  theoretical  value  of  a  can  aid  to  choose  a 
correct  value  for  X^,  in  cases  where  the  shot  noise  of  the  background  is  limiting.  In 
most  practical  situations,  however,  this  would  need  to  be  confirmed  by  performing 
many  blank  measurements  and  observing  the  effect  of  different  X^  values  on  the 
number  of  false  positives. 
Sliding  sum  method 

The  "shding  sum"  signal  processing  method  is  an  intuitively  obvious  technique. 
After  the  raw  data  is  collected,  which  consists  of  the  counts  as  a  function  of  time 
during  T„,  a  data  transformation  is  applied  in  which  the  value  sum  of  the  number 
of  counts  over  a  width  t„  is  assigned  to  the  middle  of  the  time  window.  The  window 
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is  then  moved  one  step  and  the  process  is  repeated.  For  continuous  monitoring  of 

atoms  with  CW-LIF,  the  step  size  should  be  a  fraction  of  the  residence  time  of  the 
atom  within  Vpi  e.g.,  At «  Tj/20.  Smaller  step  sizes  are  of  course  better,  but  only  up 
to  the  point  where  the  signal  processing  step  becomes  too  long.  The  optimum  value 
for  the  duration  of  the  moving  sum  is  t^  «  t,  (and  no  longer  than  the  largest  possible 
residence  time). 

Peak  Detection.  The  sliding  sum  peak  maximum  will  occur  when  t^  and  t^ 
exactly  coincide.  Thus,  it  would  seem  desirable  to  use  the  sliding  sum  peak 
maximum  value  as  a  signal  of  the  presence  or  absence  of  the  analyte  atom.  The 
distribution  of  peak  maximum  will  follow  the  distribution  of  I,;  thus,  application  of 
the  SAD  theory  from  the  previous  section  is  straightforward.  When  the  maximum 
exceeds  X^  (chosen  based  on  the  values  of  a,  ^b  and  T„),  then  the  sliding  sum  peak 
is  presumed  to  be  due  to  the  presence  of  an  atom  in  V^.  In  a  similar  manner,  peak 
maxima  can  be  used  to  evaluate  the  detection  efficiency  and  determine  if  the  method 
is  truly  SAD. 

A  problem  with  using  the  peak  value  as  an  indicator  of  I,  will  occur  when  the 
analyte  concentration  is  high  enough  so  that  there  are  problems  with  peak  overlap; 
ie,  there  is  more  than  one  atom  in  Vp  at  a  given  time.  Sliding  sum  peak  detection 
for  counting  Np  during  T„  is  only  practical  when  there  is  a  very  small  probabiUty  that 
such  an  overlap  will  occur.  '  ''      •?       ,  *-v  • 

False  positives  with  sliding  sum  peak  detection.  Calculation  of  the  theoretical 
value  of  a  when  using  sliding  sum  peak  detection  is  compUcated;  a  good  first 
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approximation  can  be  calculated  with  the  use  of  the  Gamma  distribution  (if  a  Poisson 

distribution  of  the  background  can  be  assumed)  if  X^  is  large  enough  that  F(\>X^) 
is  negUgible.  The  details  will  not  be  given  here,  but  the  most  practical  method  of 
detennining  a  for  given  values  of  X^,  0^,  T„,  and  t^  is  from  the  distribution  of  peak 
maxima  from  repeated  blank  measurements.  Theoretically,  the  values  of  a  from  this 
method  and  the  "photon  burst"  method  are  very  similar  for  the  same  conditions.  A 
comparison  of  the  number  of  false  positives  for  these  two  methods  will  be  shown  in 
chapter  8. 

Peak  area  detection.  The  problem  with  overlapping  peaks  when  using  peak 
heights  of  the  sHding  sums  was  mentioned  above.  One  method  of  alleviating  this 
problem  is  to  simply  use  the  integrated  peak  area  of  sliding  sum  peaks  as  the  "signal"; 
when  the  residence  times  of  two  atoms  overlap,  the  resulting  peak  will  be  longer  and 
the  area  will  be  the  sum  of  the  contributions  of  the  signal  due  to  both  atoms. 

A  single  photoelectron  count  in  the  raw  data  array  results  in  a  contribution 
of  t„/At  counts  to  the  resulting  sliding  sum  peak  area  as  the  window  is  moved  passed, 
where  At  is  the  step  size.  Thus,  a  cluster  of  I^  counts  (due  to  analyte  atoms  and 
background  noise)  will  result  in  a  peak  of  area  (t^/At)I,.  In  general,  the  distribution 
of  peak  areas,  I^,  will  have  the  following  characteristics 


^^^    '  [5.29] 
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where 

CTg  =  standard  deviation  of  peak  area. 

The  advantage  of  using  peak  area  detection  instead  of  peak  heights  is  that, 
with  no  loss  of  information  or  detection  efficiency,  higher  concentrations  of  analyte 
can  be  analyzed.  The  total  peak  area  from  overlapping  atoms  is  a  measured  of  the 
number  of  atoms  which  contributed  to  the  peak.  The  theory  from  the  counting 
precision  can  be  directly  appUed  in  this  situation  in  order  to  count  the  number  of 
atoms  which  contribute  to  a  given  sliding  sum  peak;  i.e.,  it  may  be  possible  to 
precisely  count  the  number  of  atoms  contributing  to  a  given  peak  (in  addition  to  the 
total  number,  Np,  which  pass  through  Vp  during  T„)  if  the  sensitivity  is  high  enough. 

Overall  Efficiency  of  Detection  "      ,  C     •'*^  > 

Much  of  this  chapter  has  been  concerned  with  detection  efficiency;  indeed,  an 
SAD  method  is  defined  essentially  as  a  method  which  has  almost  unity  detection 
efficiency.  However,  it  is  important  to  recall  the  e^,  applies  only  to  atoms  which 
interact  with  the  laser  beam.  A  high  detection  efficiency  does  not  necessary  result 
in  a  technique  with  a  corresponding  low  LOD  in  terms  of  bulk  concentration  of 
analyte  in  the  sample,  since  it  is  quite  possible  that  most  of  the  analyte  atoms  in  the 
sample  never  interact  with  the  laser  beam.  A  better  parameter  to  indicate  the 
detection  power  of  a  laser-based  method  is  the  overall  efficiency  of  detection  [80,  81], 
e^,  which  is  given  by  a  product  of  efficiency  terms: 


e,  -  e„e^e,e^  [5.30] 

where 

e,  =  the  efficiency  of  atomization, 

Cp  =  the  spatial  probing  efficiency, 

e^  =  the  temporal  probing  efficiency,  and 

ej  =  the  detection  efficiency  of  atoms  which  interact  with  the  laser. 

The  atomization  efficiency  describes  the  probability  that  an  analyte  atom  in 
the  sample  will  be  converted  into  free  atoms  in  an  energy  level  suitable  for 
interaction  with  the  laser.  In  the  case  of  molecules,  this  term  is  the  probabiUty  that 
the  analyte  molecules  in  the  sample  will  be  prepared  (perhaps  by  a  chemical  reaction 
to  tag  the  molecules  with  a  fluorophore)  in  a  state  suitable  to  produce  an 
analyte-specffic  signal.  The  spatial  probing  efficiency,  e^,  describes  the  probability 
that  the  free  analyte  atoms  will  pass  through  Vp,  the  portion  of  the  analyte  "volume" 
probed  by  the  laser.  The  temporal  probing  efficiency,  e„  is  the  probability  that  an 
atom  which  passes  through  Vp  will  interact  with  the  laser.  This  fraction  of  analyte 
atoms  is  determined  by  the  duty  cycle  of  the  laser  and  the  length  of  T^.  Finally,  the 
detection  efficiency  is  of  course  a  familiar  term  by  now;  this  term  in  the  above 
equation  takes  into  account  the  probability  of  detecting  a  signal  above  the 
background  noise  level. 

Figure  16  depicts  the  various  stages  of  a  typical  laser-based  analytical 
measurement.  Note  that  the  detection  efficiency,  e^,  as  defined  earlier,  only  comes 
into  play  at  the  last  stage.   We  have  discussed  "SAD"  methods  in  this  chapter,  but 
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fig.  16  indicates  the  difiference  between  SAD  in  Vp  and  true  single  atom  detection 

within  the  sample.  For  such  a  feat  to  be  possible,  e^  must  be  nearly  unity.  Figure 
16  also  shows  clearly  that  increasing  e,  or  e^  at  the  expense  of  any  other  term  in  eqn. 
5.30  will  not  necessarily  result  in  a  more  sensitive  analytical  technique,  even  if  SAD 
is  achieved  thereby.  For  example,  focusing  the  laser  beam  may  result  in  a  higher 
value  for  c^  but  may  decrease  e^  and  result  in  a  lower  value  for  e^. 

Conventional  LOP  for  SAD  Methods 

The  scope  of  a  given  SAD  method  can  be  used  to  determine  the  value  for 
LOD  in  terms  of  bulk  analyte  concentration  if  the  overall  efficiency  of  detection  is 
known.  Obviously,  for  the  maximum  value  for  T^  over  which  SAD  is  possible  (ie  the 
scope  of  the  method),  it  is  possible  to  detect  every  atom  which  appears  within  Vp  and 
interacts  with  the  laser.  If  we  have  knowledge  of  the  amount  of  sample  consumed 
during  T„,  and  have  reasonable  estimates  of  the  efficiency  terms  in  eqn.  5.30  then  it 
seems  that  we  should  be  able  to  calculate  the  value  of  the  LOD  in  terms  of  bulk 
concentration  of  analyte  in  the  sample,  rather  than  in  terms  of  atoms  which  interact 
with  the  laser. 

If  the  product  e  e^e^  is  known  for  an  SAD  method  and  is  much  less  than 
unity,  then  the  value  for  the  LOD  of  the  technique  is 

LOD  -   (e^e.e,)-^  [5.31] 
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in  terms  of  atoms  of  analyte  per  volume  of  sample  consumed  during  T„.  Note  that 

when  analyzing  large  samples  which  continuously  flow  through  Vp  (as  opposed  to 
discrete  sample  amount),  the  value  of  LOD  in  terms  of  bulk  concentration  of  analyte 
per  sample  analyzed  will  be  improved  through  the  use  of  larger  values  of  T^,  even 
though  the  method  may  no  longer  be  considered  SAD.  The  reason  for  this 
improvement,  of  course,  is  that  the  amount  of  sample  analyzed  increases  linearly  with 
Tm  while  (in  the  shot  noise  limit)  the  noise  of  the  background  has  a  square  root 
dependence  on  T^^. 


CHAPTER  6 
EXPERIMENTAL 


Monte  Carlo  computer  simulations  of  simple  analytical  experiments  were  used 
as  a  means  for  demonstration  and  verification  of  the  theoretical  work  set  forth  in 
chapters  four  and  five.  These  simulations,  particularly  those  based  on  the  SAD 
model,  were  also  helpful  in  the  formulation  of  the  theory  which  has  been  presented. 
The  intention  of  experiments  based  on  these  computer  simulations  is  not  to  perform 
an  exhaustive  study  of  all  aspects  of  the  models  discussed,  but  merely  to  prove  their 
vaUdity  and  provide  some  insight  into  their  usefulness.  This  chapter  will  outline  the 
general  form  of  some  of  the  programs  which  were  used  in  these  simulations;  more 
details  will  be  given  when  appropriate. 

General 

All  Monte  Carlo  simulations  were  carried  out  on  IBM-compatible 
microcomputers  which  were  equipped  with  either  a  25  MHz  80386  CPU  and  a  80387 
math  co-processor  chip,  or  a  20  MHz  80286  CPU  with  a  80287  co-processor.  All 
programs  were  written  and  compiled  in  Microsoft  QuickBASIC  (version  4.5, 
Microsoft  Corp.,  Redmond,  WA).  Algorithms  to  generate  random  numbers 
according  to  normal  distributions  and  Poisson  distributions  were  written  according 
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to  guidelines  presented  by  Knuth  [88];  these  routines  are  also  presented  in  the 

Appendix,  and  are  based  on  QuickBASIC's  pseudo-random  number  generator. 

Simulations  to  Investigate  the  Variance  of  the  LOP 

Chapter  4  presented  the  concept  of  the  calculated  LOD  of  an  analytical 
procedure  as  a  variable  estimator  of  the  true,  unknown  LOD  of  the  technique.  Two 
different  types  of  models  were  used  to  investigate  the  properties  of  the  LOD 
estimator  given  by  eqn.  4.2,  in  light  of  eqns.  4.9  and  4.10:  (1)  a  generic  situation  with 
arbitrary  standards  and  various  background  noise  levels;  and  (2)  a  model  based  upon 
conditions  found  in  the  trace  analysis  of  metals  by  electrothermal  atomization  in  a 
commercial  graphite  furnace  and  laser-induced  fluorescence  detection  of  the  analyte 
atoms  (ETA-LIF).  The  ETA-LIF  model  conditions  were  based  on  the  recent 
analysis  of  thallium  at  the  sub-femtogram  level  [89].  The  variation  of  the  LOD 
estimate  was  determined  under  different  experimental  conditions;  specifics  will  be 
given  when  appropriate. 

Experimental  Determination  of  Shot  Noise 

In  order  to  simulate  conditions  of  a  typical  ETA-LIF  experiment,  it  was 
necessary  to  have  a  knowledge  of  the  shot  noise  in  a  photomultiplier  for  a  given 
anodic  current  output.    Although  in  most  cases  the  shot  noise  can  be  calculated  1 

theoretically  using  well-known  formulas  [90],  in  the  case  of  boxcar  detection  the  task 
is  considerably  more  difficult  since  the  "effective"  electronic  bandwidth  is  unknown. 
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An  experimental  estimate  of  the  relationship  between  output  signal  and  shot  noise 
was  found  in  the  following  manner. 

An  LED  was  used  as  a  light  source  due  to  its  very  stable  output;  thus,  any 
flicker  component  in  the  source  was  minimized.  A  monochromator  was  tuned  to  the 
maximum  of  the  LED,  and  various  calibrated  neutral  density  filters  were  placed 
between  the  LED  source  and  the  entrance  slit  of  the  monochromator.  A  terminator 
of  1.2  kn  was  used  to  convert  the  anodic  current  of  the  photomuhiplier  (model  R928, 
Hamamatsu  Corporation,  Bridgewater,  NJ)  to  a  voltage;  the  signal  was  amplified  20X 
by  a  voltage  amplifier  (Evans  model  4163,  Evans  Associates,  Berkeley,  CA)  and  input 
into  a  boxcar  averager  (SRS  model  SR250,  Stanford  Research  Systems,  Sunnyvale 
CA).  The  boxcar  was  internally  triggered  at  50  Hz,  and  a  gate  of  300  ns  was  used. 
Various  values  for  the  averaging  was  used:  last  sample,  3  and  10  averages.  More 
averages  were  not  attempted  because  an  overly  large  time  constant  would  be 
inappropriate  for  the  transient  pulses  in  ETA-LIF.  The  conditions  chosen  above 
closely  mimic  those  which  are  found  in  a  typical  ETA-LIF  experiment. 

The  shot  noise  is  proportional  to  the  square  root  of  the  mean  output  anodic 
current;  the  proportionality  constant  was  determined  by  linear  least-squares  fit  to  a 
plot  of  the  variance  vs.  the  anodic  current.  The  fitted  line  relating  shot  noise  and 
anodic  current  was  used  to  find  the  shot  noise  at  any  signal  level  in  the  simulations 
of  ETA-LIF  conditions. 
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Experimental  Variation  in  Estimated  LOD  "^ 

Using  the  model  conditions  for  the  background,  sensitivity  and  signal  noise, 
a  computer  program  was  written  which  would  calculate  s^  and  slope  for  1000 
calibration  curves,  and  estimate  the  LOD  according  to  eqn.  4.2  for  each  one.  The 

A. 

observed  value  of  Slod  calculated  from  the  1000  values  of  LOD  can  be  compared  to 
the  theoretical  value  predicted  by  eqn.  4.10.  In  addition,  the  performance  of 
confidence  intervals  about  the  estimated  LOD  can  be  evaluated  by  how  often  the 
true  LOD  value  falls  within  the  interval.  The  following  steps  give  a  brief  outline  of 
the  algorithm  used: 

1.  The  model  parameters  were  defined:  true  values  for  the  sensitivity, 
background  mean  and  standard  deviation,  and  calibration  standards.  Thus, 
the  "true"  value  for  the  LOD  is  known. 

2.  The  signal  level  for  each  calibration  standard  was  "measured".  For  both 
generic  and  ETA-LIF  models,  the  noise  on  the  signal  was  assumed  to  be 
normally  distributed. 

3.  The  slope  and  expected  standard  error  in  the  slope  were  calculated  with 
either  unweighted  or  weighted  least-squares  regression,  using  conventional 
formulas  [20]. 

4.  For  each  of  the  1000  "experiments",  a  number  of  separate  "measurements" 
were  taken  of  the  blank  using  the  normal  random  number  generator. 


% 

Estimates  of  a^  were  calculated  from  the  results;  however,  the  blanks  were 
not  included  in  the  calibration  curve. 

5.  For  each  experiment,  the  LOD  and  Jlod  values  were  estimated  according  to 
eqns.  4.2  and  4.10,  and  a  95%  confidence  interval  was  formed  assuming  either 
a  normal  distribution  of  estimated  LOD  values  or  with  the  use  the  x^  tables, 
as  in  eqn.  4.11.  If  the  "true"  LOD  was  not  within  the  interval,  then  this  fact 
was  noted. 

6.  The  estimated  variance,  Slod>  ^^om  the  1000  estimates  of  the  true  LOD  values 
was  calculated. 

As  noted  earlier,  two  models  were  used  for  the  experimental  simulations  of 
calibration  curves  and  blank  measurements:  a  generic  model  to  examine  the  behavior 
of  eqns.  4.9  and  4.10  under  various  conditions,   and  a  model  based  on  the 
determination  of  thallium  by  ETA-LIF.  These  models  will  now  be  described. 
Simulations  using  a  generic  model 

The  general  validity  of  eqn.  4.10  was  investigated  using  an  arbitrary  model. 
The  effect  of  increasing  slope  error  and  decreasing  numbers  of  blank  measurements 
(i.e.,  increasing  either  term  in  eqn.  4.9  which  contribute  to  cTlod)  was  investigated  to 
determine  the  range  within  which  the  expressions  may  be  used.  The  model 
conditions  are  as  follows:  ■,.,.. 

Standards  (x-axis):  220,  240,  260,  280,  300  (arbitrary  units). 

Sensitivity:  1. 
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All  measurements  exhibited  homogeneous  variance,  a^.  The  value  of  a  ;. 
determines  the  LOD  and  the  value  of  CV^  in  eqn.  4.9.  Two  different  trends  were 
followed  to  determine  their  effect  on  the  calculated  and  observed  value  of  Slod^  ^^st 
the  slope  error  was  gradually  increased  by  increasing  the  value  of  a  (and  the  true 
LOD)  while  the  number  of  blank  measurements  was  held  constant  at  twenty.  Next, 
the  number  of  blank  readings  used  to  estimate  \  was  varied  while  the  value  of  a  was 
held  constant. 
Simulations  using  an  ETA-LIF  model 

A  model  based  on  the  determination  of  thallium  by  ETA-LIF  was  used  to 
demonstrate  eqn.  4.10  in  a  shghtly  more  realistic  setting.  The  sensitivity  was  2.5 
mV/fg  (peak  detection).  The  limiting  noise  was  either  the  photomultipher  dark 
current  (in  which  case  the  true  LOD  would  be  0.12  fg)  or  laser-induced 
scatter/fluorescence  (LOD  =  1.367  fg).  Wei  et  al.  have  reported  that  in  the  second 
case,  the  noise  is  white  in  nature  [91];  it  is  assumed  here  to  be  due  to  photomultipher 
shot  noise,  and  the  magnitude  is  calculated  as  explained  in  the  experimental  section. 
The  noise  on  the  "signal"  was  due  to  two  sources:  (1)  photomultipher  shot  noise,  as 
just  described,  and  (2)  injection  imprecision  of  the  sample.  The  first  contribution  was 
important  at  lower  concentrations.  The  second  contribution  was  assumed  to  have  a 
value  of  5%  RSD.  The  program  would  calculate  the  signal  for  a  given  standard  as 
follows:  first  the  true  injected  amount  of  sample  would  be  calculated  by  assuming  a 
normal  distribution  (5%  RSD)  centered  on  the  true  standard  value.  The  mean  signal 
for  this  injected  amount  was  calculated  from  the  true  (model)  sensitivity,  and  the 
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actual  "measured"  signal  was  normally  distributed  around  the  mean  with  a  variance 

due  to  photomultipher  shot  noise. 

Weighted  least  squares  fitting  was  used  for  all  program  runs  in  which  a 
calibration  curve  was  used  to  determine  the  sensitivity.  It  is  important  in  practice  to 
determine  whether  weighted  or  unweighted  least  squares  is  necessary.  The  second 
term  in  eqn.  4.9  wiU  not  be  a  good  estimation  of  the  slope  coefficient  of  variation  if 
unweighted  least  squares  is  incorrectly  used  to  fit  heteroscedastic  data. 

There  were  three  different  types  of  simulations  run  with  the  ETA-LIF  model. 
In  the  first  case,  ten  standards  were  each  "measured"  once.  The  standards  covered 
5  orders  of  magnitude  from  50  fg  to  1  ng  injections  into  the  furnace.  These  values 
are  realistic  for  the  analysis  of  thallium,  which  is  present  at  very  low  levels  in  the 
environment  --  contamination  of  the  standards  is  not  a  problem.  Weighted  least 
squares  regression  was  used  with  the  weights  equal  to  the  inverse  of  the  true  variance 
of  the  signals  at  each  concentration.  In  the  second  run,  only  3  standards  covering  2 
orders  of  magnitude  (1  pg  -  100  pg)  were  used,  but  each  was  measured  four  times 
to  provide  estimates  for  the  weights  which  should  be  used  in  weighted  least  squares 
regression.  This  is  a  more  realistic  situation  than  the  case  in  the  first  run,  where  the 
true  weights  were  known.  In  the  final  run,  calibration  curves  were  not  used  but 
instead  a  single  calibration  standard  (1  pg)  was  measured  four  times  to  provide  an 
estimate  for  the  sensitivity  to  use  in  estimating  the  LOD  by  eqn.  4.2.  The  CV^  term 
in  eqn.  4.9  was  estimated  by  accounting  for  the  error  in  both  signal  and  background 
averages,  and  then  using  eqn.  2.14  to  estimate  a.. 
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■"  -^ 
V  Simulations  of  SAD  bv  LIF 

Chapter  5  presented  a  general  model  of  an  SAD  technique  by  laser 
spectroscopic  methods,  and  outlined  various  associated  concepts  of  SAD  methods 
based  on  this  model.  In  Alkemade's  original  work  on  SAD  theory  [79,  80],  the 
possibihty  of  detecting  single  atoms  in  the  presence  of  noise  was  ignored,  and  an 
SAD  method  was  considered  possible  only  in  the  absence  of  noise  (i.e.,  at  the 
intrinsic-noise  lunit).  In  this  work,  a  more  general  model  broadens  the  definition  of 
a  true  SAD  method  to  include  the  presence  of  noise;  however,  SAD  is  only  possible 
under  these  conditions  by  a  nondestructive  technique.  The  computer  simulations 
which  will  be  presented  emphasize  this  fact  by  modeling  an  SAD  method  based  on 
cyclic  LDF  of  atoms  flowing  through  the  laser  beam.  Many  of  the  aspects  of  SAD  in 
the  presence  of  noise  presented  in  chapter  5  will  be  illustrated  using  these 
simulations. 

Models  for  the  Interaction  Time  with  the  Laser 

The  simulations  for  SAD  by  cyclic  LIF  were  based  on  an  experiment  in  which 
the  analyte  species  flow  through  the  probe  volume  and  interact  with  the  laser,  which 
is  a  CW  laser.  The  simulations  can  be  classified  into  two  types,  according  to  the 
method  used  to  calculate  the  residence  time,  t^  of  an  atom  within  Vp;  these  models 
will  be  called  the  Simple  and  the  Cylindrical  Probe  models.  In  all  cases,  it  was 
assumed  that  0,  was  constant  throughout  Vp.  All  other  assumptions  of  the  general 
SAD  model  presented  in  chapter  5  apply. 
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Simple  model 

In  simulations  based  on  the  simple  model,  it  was  assumed  that  the  value  of 
tj,  and  hence  l^,  is  constant  for  all  atoms  which  enter  Vp.    In  this  situation,  the 

It*' 

variable  I,  is  a  Poisson  variable,  and  the  "extra"  variance  term  in  eqn.  5.13  is  zero. 
Such  a  situation  can  occur  with  pulsed  laser  LIF  experiments  or  in  situations  where 
the  shape  of  Vp  is  such  that  the  analyte  atoms  have  (ahnost)  constant  t^  values  for 
any  path  through  Vp  -  e.g.,  if  a  cylindrical  lens  was  used  to  focus  the  laser  in  a  line 
across  the  analyte  stream,  or  the  analyte  flow  was  constricted  so  that  it  all  flowed 
through  the  center  of  an  unfocused  or  expanded  laser  beam  (e.g.,sheath  flow).  The 
use  of  spatial  masking  to  shape  Vp  might  also  decrease  the  variation  in  t,  to  a  very 
low  value. 
Cylindrical  probe  volume  model 

The  purpose  of  this  model  is  not  necessarily  to  provide  a  more  realistic 
simulation  of  typical  experimental  conditions  but  to  demonstrate  the  effect  of  a 
non-zero  "extra"  variance  in  eqn.  5.13.  In  this  model,  the  extra  variance  is  due  to  a 
variable  interaction  time  with  the  laser  beam.  The  model  is  based  on  an  experiment 
such  as  was  shown  in  figure  9  in  chapter  5,  in  which  the  probe  volume  as  defined  by 
the  laser  beam  and  the  collection  optics  (with  a  slit  as  a  limiting  apermre)  is  a 
cylinder. 

In  this  model,  it  is  assumed  that  tj  is  a  variable  that  depends  on  the  path  of 
an  individual  atom  through  Vp.  The  situation  used  in  the  simulations,  along  with 
labelled  axes,  is  shown  in  figure  17.    Effects  due  to  diffusion,  flow  profiles  and 
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velocity  distribution  are  neglected;  all  atoms  enter  Vp  with  the  same  velocity  in  the 
z-direction,  v^.  Both  v^  and  Vy  components  of  the  velocity  are  assumed  to  be  zero. 
The  atoms  are  spread  uniformly  throughout  the  analyte  volume,  V„  which  flows  past 
V  during  T„.  Since  the  number  of  atoms  within  this  volume  will  vary  according  to 
a  Poisson  distribution,  the  number  of  atoms  which  interact  with  the  laser,  Np,  would 
normally  also  follow  a  Poisson  distribution  in  an  experiment  following  this  model. 
The  Cartesian  coordinate  system  shown  in  fig.  17  has  its  origin  in  the  center  of  Vpi 
if  the  position  (x,y)  of  the  atom  within  V,  is  such  that  x  <  Xp  and  y  <  rp  then  this 
atom  interacts  with  the  laser.  For  an  atom  which  enters  Vp,  the  value  of  tj  is 
dependent  only  on  the  position  along  the  y-axis  (since  v^^  =  Vy  =  0). 

For  this  model,  it  is  easy  to  see  that  the  probing  efficiency  is  given  by      '  . 


,  -  ^  [6.1] 


4X, 


The  interaction  time  of  an  atom  which  enter  Vp  can  also  easily  be  calculated  for  the 
Cylindrical  V  model  according  to 


u  -  ^{^^l-y^  [6.2] 


Vz 


Note  that  t;  depends  only  on  the  (uniformly  random)  position  of  the  atom  along  the 
y-axis.  It  is  assumed  that  y  is  uniformly  distributed  within  Vp^;  under  these 
conditions,  the  mean  value  of  t  can  be  calculated  to  be  ^        . -<    •  - . 


2  „  2\H 


Hhe  condition  necessary  for  this  is  that  rp  <  {r^  -Xp  )  . 
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In  light  of  the  Cylindrical  V^  model,  the  Simple  model  with  fixed  tj  as 
described  above  can  be  viewed  as  a  "rectangular"  Vp  model  with  v^  =  Vy  =  0  and 
constant  v^,  as  above.  Comparisons  between  the  two  models  were  made  by  choosing 
and  r  (or  vj  in  eqn.  6.3  such  that  Tj  for  the  Cylindrical  Vp  model  is  equal  to  the 
fixed  value  of  t;  assumed  in  the  Simple  model  --  in  this  manner  the  mean  signal  per 
atom,  \  ,  remained  the  same  for  both  models. 

General  Outline  of  Simulations 

For  either  of  the  models  described,  the  general  format  of  a  typical  simulation 
used  to  generate  the  raw  data  (detected  photons  as  a  function  of  time)  was  as 
follows: 

1.  The  model  parameters,  T„,  (p^,,  <p^,  X^,  were  defined.  The  signal  detection 
limit,  Xj,  was  set  by  the  values  of  <p^„  T„  and  the  desired  a;  different  values 
of  Xj  were  often  chosen  to  investigate  the  effects  on  the  number  of  false 
positives. 

2.  Values  for  Np  and  t;  were  chosen  or  calculated.  Although  Np  would  normally 
follow  a  Poisson  distribution  in  an  SAD  experiment,  it  was  usually  fixed  at  a 
given  value  for  the  computer  simulations.  The  value  of  tj  was  assigned 
according  to  the  model  used;  if  a  Cylindrical  Vp  model  was  used,  then  it  was 
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calculated  according  to  eqn.  6.2  with  the  position  on  the  y-axis  a  uniform 

random  variable  between  -t^  to  t^. 

3.  The  contribution  of  the  background  noise  was  added  to  the  data  matrix.  The 
number  of  counts  due  to  the  noise  varied  according  to  a  Poisson  distribution 
with  a  mean  of  (p^T^.  If  this  mean  was  high  enough  (above  50),  a  normal 
distribution  was  used  to  approximate  the  Poisson  distribution.  The  detection 
times  for  each  of  these  photoelectrons  was  then  generated  according  to  a 
uniform  random  distribution  between  0  and  T^. 

4.  The  entry  time  of  each  Np  atom  into  Vp  was  randomly  chosen  from  a  uniform 
distribution.  The  signal  due  to  each  of  these  atoms  was  a  Poisson  variable 
with  mean  (p^t^,  where  l,  was  calculated  for  each  atom  according  to  the  model 
used.  The  signal  counts  were  uniformly  distributed  during  the  atom's 
residence  time  within  Vp. 

Once  the  raw  data  file  was  generated  according  to  this  procedure,  any  signal 
processing  method  could  be  chosen  with  an  appropriate  value  of  t^  One  of  the 
advantages  of  using  a  computer  simulation  such  as  the  one  above  to  investigate  SAD 
theory  is  that  the  number  Np  could  be  fixed,  something  which  may  be  impossible  to 
achieve  in  an  SAD  experiment.  The  number  of  "positives"  -  signals  at  or  above 
Xj  --  could  be  compared  to  the  true  number  of  Np.  This  ability  was  necessary  to 
investigate  the  behavior  of  the  counting  precision,  RSD^,  at  various  values  of  0^,  <p^, 
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and  Np.  Also,  the  performance  of  the  signal  processing,  the  sliding  sum  and  photon 

burst  methods,  could  be  more  easily  evaluated. 

In  many  instances,  simpler  versions  of  the  above  simulation  were  used.  For 
example,  when  investigating  RSD^  the  entire  raw  data  file  was  not  necessary;  only  the 
analyte  signal  was  of  interest.  In  such  cases,  the  value  of  I,  with  the  analyte  in  the 
probed  region,  was  a  Poisson  variable  with  mean  /i^,  +  N  05ti,  and  X^  was  chosen 
according  to  the  value  of  Mb-  This  is  similar  to  the  situation  of  the  "simple 
integration"  signal  processing  presented  in  chapter  5,  with  the  value  of  N^  fixed  at  a 
given  level. 


CHAPTER? 

INVESTIGATING  THE  VARIABILITY  OF  THE 

LIMIT  OF  DETECTION 


Introduction 

In  chapter  4,  the  variability  of  the  calculated  LOD  of  an  analytical  procedure 
was  explained  to  be  the  result  of  two  factors:  true  changes  in  the  population 
parameters  of  the  procedure  (such  as  sensitivity  and  background  noise)  and  the 
imprecision  of  the  sample  statistics  used  to  estimate  the  true  (unknown)  LOD  of  the 
system.  In  this  chapter  computer  simulations  are  used  to  investigate  the  magnitude 
and  nature  of  this  second  source  of  fluctuation  in  estimated  LOD,  and  particularly 
to  illustrate  the  limitations  and  uses  of  eqns.  4.9  and  4.10  derived  in  chapter  4. 

The  first  set  of  simulations  investigate  the  validity  of  the  derived  equations 
with  various  magnitude  of  contribution  from  the  two  terms  in  eqn.  4.9.  These 
experiments  use  the  generic  model  presented  in  chapter  6.  The  application  and 
utility  of  these  equations  in  a  typical  analysis  of  ultra-trace  levels  of  thallium  by 
ETA-LTF  is  demonstrated  by  the  second  set  of  simulations. 
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Effect  of  Increasing  Values  of  Slope  Error 

The  effect  of  increasing  values  of  slope  error  on  observed  and  calculated  Slod 
in  simulations  run  using  the  generic  model  is  shown  in  table  2.  The  program  takes 
30  blank  "measurements"  for  each  simulation.  The  value  of  the  parameters  a^,  and 
CVj  are  shown  in  the  first  two  columns.   The  entries  in  the  third  column  are  the 

A. 

average  value  and  standard  deviation  of  LOD  calculated  for  1000  caUbration  curves. 
The  true  LOD  value  and  calculated  ctlod  (using  eqn.  4.10  with  the  population 
parameters)  are  shown  in  column  four.  The  confidence  intervals  shown  in  the  table 
are  calculated  from  tables  of  the  %  distribution  or  with  the  formula  LOD  ±  2Slod- 
The  percentage  of  runs  in  which  the  true  LOD  value  lies  below  the  lower  limit  are 
given  in  the  colunm  labelled  "lower";  the  percentage  above  the  upper  limit  are  given 
in  the  next  column.  The  expected  values  in  both  columns  is  about  2.5%. 

The  agreement  between  the  predicted  and  actual  behavior  of  the  LOD 
estimated  by  eqn.  4.10  is  quite  close  while  the  slope  error  remains  relatively  small. 
Above  a  value  of  CV^  of  roughly  10%,  the  x^  confidence  intervals  begin  to  show  a 
noticeable  deviation  from  the  expected  values;  at  this  point,  however,  the  observed 
and  calculated  values  of  Slqd  still  largely  agree.  The  final  two  runs,  however,  show 
a  marked  departure  of  the  observed  Slqd  froni  the  value  predicted  by  eqn.  4.10;  this 
is  due  to  the  approximate  nature  of  the  propagation  of  errors  approach  since  the 
LOD  estimate  is  non-linearly  related  to  the  slope  value.  Note  that  at  higher  values 
of  CVa,  the  2Slod  intervals  seem  to  perform  better  than  the  x^  intervals,  since  the 
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Generic  case:  effect  of  increasing  a^  (homogenous  variance) 


^b 

CVa 

LOD 
estimate^ 

LOD 

calc'd'' 

X^  intervals'^ 

2Slod  intervals'* 

lower 

upper 

lower 

upper 

1 

1.6% 

2.98 
0.39 

3.00 
0.40 

2.5% 

1.5% 

0.6% 

5.4% 

2 

3.2% 

5.97 
0.77 

6.00 
0.81 

2.7% 

1.7% 

0.1% 

4.2% 

3 

4.7% 

8.91 
1.20 

9.00 
1.13 

3.5% 

1.3% 

0.3% 

3.2% 

5 

7.9% 

14.92 
2.38 

15.00 
2.30 

5.9% 

3.6% 

0.8% 

6.0% 

7.5 

11.9% 

22.72 
3.97 

22.50 
3.98 

8.7% 

4.8% 

0.8% 

5.8% 

10 

15.8% 

30.63 
6.31 

30.00 
6.17 

12.5% 

6.0% 

0.7% 

5.4% 

12.5 

19.8% 

38.94 
11.31 

37.5 
8.9 

16.9% 

10.0% 

0.7% 

8.1% 

15.0 

23.7% 

49.74 
47.40 

45.0 
12.2 

21.3% 

13.1% 

1.3% 

8.9% 

^calculated  with  eqn.  4.2  for  each  loop.   First  number  in  each  entry  is  the  average 
of  1000  loops;  second  is  the  observed  standard  deviation. 

''theoretical  values  calculated  with  eqns.  4.1  and  4.10. 

'^e  X   confidence  intervals  are  calculated  for  each  loop  from  the  %   tables  with  29 
degrees  of  freedom  and  a  =  2.5%  on  either  side. 

^95%  intervals  calculated  using  eqn.  4.10  and  eqn.  4.2  for  each  loop. 
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CV  term  in  eqn,  4.9  begins  to  become  more  important.  In  the  next  section  there 

will  be  more  discussion  on  the  two  types  of  intervals. 

Effect  of  Increasing  Number  of  Blank  Measurements 

For  this  next  set  of  simulations  using  the  generic  model,  the  value  of  CV^  was 
kept  at  about  3%  while  N,  the  number  of  blank  measurements,  was  increased.  The 
results  are  shown  in  table  3;  the  calculated  Slod  in  the  third  column  is  the  average 
of  the  1000  estimates  for  each  value  of  N  in  the  first  column.  It  appears  that  for  all 
values  of  N,  the  variabihty  of  the  LOD  estimate  is  well  predicted  by  eqn.  4.10.  One 
aspect  of  table  3  is  immediately  obvious:  the  LOD  estimates  shown  in  column  two 
are  all  negatively  biased  compared  to  the  true  LOD  values  listed  in  the  third  column. 
Upon  examination  of  eqn.  4.2,  however,  the  bias  in  the  LOD  estimate  is  perfectly 
logical,  since  the  sample  standard  deviation  is  a  biased  estimator  of  the  population 
standard  deviation  (although  the  sample  variance  is  unbiased)  [92].  Thus,  the 
calculated  LOD  is  a  biased  estimator  of  the  true  LOD  of  the  analytical  procedure. 
The  bias  is  worse  with  fewer  degrees  of  freedom;  this  is  one  reason  for  choosing  a 
relatively  large  number  of  blank  measurements  to  calculate  s^.  However,  the  bias 
was  small  relative  to  the  value  of  Slod  ^^  any  value  of  N  used;  hence,  the  decrease 
in  Slod  with  increasing  N  is  still  the  biggest  motivation  for  a  large  number  of  blank 
measurements  when  estimating  the  LOD. 

The  confidence  intervals  obtained  here  with  the  x^  distribution  are  superior 
to  those  obtained  by  assuming  a  normal  distribution  of  the  LOD  estimate;  however, 


Table  3 

The  effect  of  increasing  numbers  of  blank  measurements 
on  LOD  estimates  (CV^  kept  at  3%) 
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N 

LOD 

estimate^ 

LOD 

calc'd^ 

X    intervals'^ 

2Slod  intervals'^ 

Lower 

Upper 

Lower 

Upper 

3 

5.37 
2.74 

6.00 
2.69 

1.7% 

2.2% 

0.0% 

21.5% 

5 

5.60 
2.02 

6.00 
1.99 

2.4% 

2.2% 

0.0% 

16.5% 

10 

5.82 
1.29 

6.00 
1.38 

1.9% 

L6% 

0.0% 

7.5% 

15 

5.88 
1.11 

6.00 
1.13 

2.3% 

1.4% 

0.2% 

6.6% 

20 

5.96 
0.95 

6.00 
0.98 

2.5% 

2.3% 

0.2% 

5.2% 

30 

5.96 
0.76 

6.00 
0.80 

2.9% 

1.5% 

0.5% 

4.4% 

^from  eqn.  4.2. 

''from  eqns.  4.1  and  4.10. 


'^from  X    tables  with  N-1  degrees  of  freedom, 
'rrom  eqns.  4.2  and  4.10  for  each  loop. 


Ill 

the  latter  intervals  become  more  acceptable  as  N  increases  to  values  of  15  or  better. 
The  performance  of  the  confidence  intervals  based  upon  the  normal  distribution  was 
at  first  puzzling;  even  when  N  or  CV^  were  increased  to  the  point  where  the 
distribution  of  LOD  was  well  characterized  by  a  normal  distribution  (e.g.,  N  >  30) 
the  performance  of  the  confidence  intervals  did  not  improve.  In  all  cases,  there  was 
a  larger  fraction  of  cases  where  the  true  LOD  value  was  larger  than  the  upper 
confidence  limit;  at  the  same  time,  in  very  few  cases  was  the  true  LOD  value  less 
than  the  lower  confidence  Umit. 

The  distribution  of  the  1000  values  of  LOD  with  30  blank  measurements  are 
shown  in  figure  18,  along  with  a  Gaussian  distribution  centered  around  the  true  LOD 
value  and  with  a  standard  deviation  predicted  by  eqn.  4.10.  As  can  be  seen,  the 
deviation  fi-om  a  normal  distribution  is  only  slight.  Indeed,  values  of  the  x^ 
distribution  above  30  degrees  of  freedom  are  commonly  calculated  by  assuming  a 
normal  distribution.  And  yet,  there  is  asymmetry  in  all  simulations  when  2Slod 
confidence  intervals  are  used  --  the  true  LOD  tends  to  be  above  the  upper 
confidence  limit  more  than  expected.  This  behavior  is  all  the  more  puzzhng  since 
the  slight  asymmetry  shown  in  the  figure  is  skewed  to  the  right,  and  would  tend  to 
suggest  that  the  opposite  trend  would  prevail,  i.e.,  the  true  LOD  would  be  less  than 
the  lower  confidence  limits  of  the  extra  LOD  values  on  the  upper  wing  of  the 
distribution. 

The  behavior  can  be  understood  by  realizing  that,  as  shown  in  eqn.  4.10,  the 
estimated  value  Slqd  (and  hence  the  size  of  the  confidence  interval)  depends  on  the 


112 


140 


LOD  estimate 


Figure  18.  Histogram  of  calculated  LOD  values  from  1000  calibration 

curves.  The  symbols  represent  observed  frequency  of  LOD 
estimates  with  30  blank  measurements  and  CV^  =  3%.  The 
curve  is  a  normal  distribution  with  ji  =  LOD  and  a  value  of 
^LOD  calculated  by  eqn.  5.10. 
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magnitude  of  the  estimated  LOD.     Remember  that  confidence  intervals  are 

constructed  around  the  calculated  value  of  LOD;  the  true  value  of  LOD  is  unknown. 

A. 

Values  of  LOD  on  the  lower  tail  of  the  distribution  in  fig.  18  have  smaller  confidence 
intervals  than  the  values  on  the  upper  tail;  thus,  a  disproportionate  number  of  values 
on  the  lower  side  of  the  distribution  of  LOD  will  have  confidence  intervals  which 
do  not  include  the  true  mean  of  the  distribution  (ie  the  true  LOD).  Likewise,  values 
on  the  right  wing  of  the  distribution  will  have  larger  confidence  intervals  and  a  larger 
number  of  these  values  will  include  the  true  LOD  value  within  the  interval. 

Using  confidence  intervals  calculated  from  y^  tables  gives  the  best  confidence 
intervals  in  the  cases  shown  in  table  3;  however,  using  LOD  ±  3Slod  confidence 
intervals  is  ahnost  certain  (>99%  from  past  experience)  to  contain  the  true  value 
LOD,  and  may  be  more  appropriate  when  the  CV^  term  is  large,  as  can  be  seen  in 
table  2.  In  addition,  in  some  cases,  very  large  values  for  N  might  be  possible 
depending  on  the  source  of  the  limiting  background  noise;  for  example,  if  the 
background  noise  is  not  sample  related  but  is  due  to,  for  example,  laser  scatter,  s^ 
can  be  measured  with  very  large  degrees  of  freedom.  In  this  case,  as  long  as  CV,  is 
low  enough,  the  confidence  interval  can  easily  be  calculated  by  using  SSlqd  limits 
from  eqn.  4.10. 

Application  to  HTA-TTF 

There  are  two  sources  of  noise  on  a  given  calibration  point  in  any  ETA 
experiment:  noise  due  to  variation  of  the  signal  (y-value)  and  noise  due  to  imprecise 
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injection  (variation  in  x-value).  It  is  normally  stated  in  any  textbook  on  regression 

that  least-squares  estimation  of  the  regression  parameters  assumes  that  there  is  no 
error  in  the  x-axis.  In  many  cases  in  analytical  chemistry,  it  is  assumed  that  errors 
in  the  independent  variable  are  negligible  compared  to  the  range  of  this  variable; 
hence,  any  random  errors  in  the  calibration  standards  are  ignored.  However,  for 
some  of  the  cases  described  here  with  ETA-LIF  this  assumption  is  not  valid. 

Berkson  [93]  has  noted  a  special  class  of  situations  when  errors  in  both 
variables  can  still  be  treated  as  a  simple  least-squares  regression.     When  the 
x-variable  can  be  controlled,  then  any  random  error  in  the  independent  variable  can 
be  attributed  to  an  additional  random  error  component  in  the  dependent  variable. 
Berkson's  case  has  been  discussed  by  Mandel  [71]  and  Neter  et  al.  [94]  and  should 
apply  in  this  case  even  though  the  error  is  not  constant  but  linearly  proportional  to 
the  magnitude  of  the  independent  variable.  Table  4  shows  the  magnitudes  of  the  two 
error  components  ~  due  to  signal  shot  noise  and  injection  imprecision  ~  and  the 
calculated  and  observed  total  noise  in  the  signal  at  each  calibration  standard.  Since 
the  observed  noise  closely  follows  the  calculated  values,  it  appears  that  Berkson's 
model  does  indeed  apply  to  ETA-LIF  (when  the  error  in  the  injected  amount  is 
randomly  distributed  about  the  expected  value);  additionally,  it  is  obvious  that  shot 
noise  in  the  signal  is  quickly  overshadowed  by  the  injection  imprecision  for  the 
parameters  used  in  this  model. 

Table  5  shows  the  results  of  the  ETA-LIF  simulation  run  under  the  three 
conditions  described  in  the  experimental  section.    The  calculated  value  for  S. 
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Table  4 

Observed  and  calculated  noise  on  signal  in  ETA-LIF  model; 
run  1(a)  (10  standards,  dark  current  limited) 
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Cone 
(Pg) 

Observed 

signal.noise 

(V) 

Expected 

signal.noise 

(V) 

Shot  noise 
(V) 

Cone  flicker 

(injection) 

(V) 

0.05 

0.125 
9.32  X  10"^ 

0.125 
9.24  x  10'^ 

6.81  X  10"^ 

6.25  X  10-^ 

0.1 

0.249 
1.525  X  10-2 

0.250 
1.575  x  10-2 

9.59  X  10"^ 

1.25  X  10-2 

0.5 

1.248 
6.633  X  10-2 

1.250 
6.603  X  10-2 

2.13  X  10-2 

6.25  X  10-2 

1 

2.509 
0.134 

2.500 
0.129 

3.01  X  10-2 

0.125 

5 

12.499 
0.623 

12.500 
0.629 

6.71  X  10-2 

0.625 

10 

24.957 
1.251 

25.000 
1.254 

9.50  X  10-2 

1.250 

50 

125.10 
6.16 

125.00 
6.25 

0.212 

6.25 

100 

250.40 
12.67 

250.00 
12.50 

0.300 

12.50 

500 

1.252  X  10^ 
62.474 

1.250  X  10^ 
62.504 

0.671 

62.50 

1000 

2.510  X  10^ 
125.00 

2.500  X  10^ 
125.00 

0.949 

125.00 

■'  'W 
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Table  5 
Estimating  LOD  for  Thallium  in  ETA-LIF 
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Case^ 

LOD 

estimate'^ 

LOD 

calc'd'^ 

X^  intervals'* 

2Slod  intervals^ 

lower 

upper 

lower 

upper 

1(a) 

0.119 
0.019 

0.120 
0.019 

2.4% 

2.0% 

0.1% 

5.6% 

1(b) 

1.343 
0.216 

1.367 
0.219 

2.4% 

2.0% 

0.4% 

5.9% 

2(a) 

0.119 
0.019 

0.120 
0.020 

3.4% 

1.6% 

0.2% 

4.9% 

2(b) 

1.351 
0.218 

1.367 
0.221 

2.8% 

2.2% 

0.4% 

4.7% 

3(a) 

0.118 
0.019 

0.120 
0.020 

2.4% 

2.5% 

0.2% 

4.9% 

3(b) 

1.352 
0.223 

1.367 
0.224 

3.6% 

1.5% 

0.4% 

5.5% 

^Case  1  uses  weighted  linear  regression  with  10  standards  and  the  true  weights 
known;  case  2  uses  only  three  standards  and  measures  each  4  times  to  estimate  the 
weights;  case  3  only  measures  one  standard  four  times  to  estimate  the  sensitivity. 
In  all  cases,  (a)  is  dark  current  limited,  and  (b)  is  laser-induced  shot  noise  limited 
(ie,  due  to  scatter  or  fluorescence).    See  the  experimental  section  for  more  details. 

''from  eqn.  4.2. 

'^from  eqns.  4.1  and  4.10.        .-  '  --       -       ,-., 

'  \  \  ' 

from  X    tables  with  19  degrees  of  freedom, 
^from  eqns.  4.2  and  4.10  for  each  loop.  * 


■  \ 
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shown  in  the  third  column  is  the  average  of  the  individual  estimates  of  each  of  the 
1000  calibration  curves.  It  is  easily  seen  that  eqn.  4.10  predicts  the  observed 
variation  in  the  LOD  estimate  quite  well;  by  examining  the  relative  magnitudes  of 
the  terms  in  eqn.  4.9,  in  all  cases,  the  first  term  dominates,  and  the  variation  in  LOD 
is  caused  for  the  most  part  by  the  imprecision  in  the  %  estimate  (N  =  20  in  all 
simulations  here).  This  fact  is  taken  advantage  of  in  the  final  case  (3a  and  3b):  only 
one  standard  is  measured  to  estimate  the  sensitivity  of  the  method  with  no  loss  in 
efficiency  of  the  LOD  estimate.  Thus,  when  changing  various  parameters  of  the 
analytical  system  (e.g.,  laser  excitation  lines  used,  monochromator  slit  width,  etc.)  to 
improve  the  LOD,  only  one  standard  need  be  measured  to  calculate  the  new  LOD 
with  no  loss  in  precision  of  the  estimate,  as  long  as  it  is  known  that  the  calibration 
function  is  linear  up  to  the  standard  used. 

All  of  the  above  calculations  of  LOD  have  used  the  lUPAC  recommendation 
for  the  confidence  factor  (i.e.,  k  =  3).  It  is  instructive  to  compare  the  true  LOD  of 
the  ETA-LIF  model  using  k=3  with  the  "calibration  band"  approach  of  calculating 
the  detection  limit,  depicted  in  fig.  4  in  chapter  2,  and  using  a  confidence  factor  given 
by  eqn.  2.16.  For  run  2(a)  (only  three  standards  in  the  calibration  set),  calculation 
of  the  confidence  factor  according  to  eqn.  2.16  gives  the  true  value  of  LOD  at  0.16 
fg,  as  opposed  to  0.12  fg.  However,  the  value  of  Slqd  calculated  was  0.02,  and  the 
upper  tolerance  limit  (2a)  is  0.16  fg  --  i.e.,  95%  of  the  values  of  calculated  LOD 
estimates  using  k  =  3  will  be  in  the  interval  0.08-0.16  fg.  Clearly,  ignoring  the 
variation  in  LOD  (due  in  this  case  to  the  variation  in  s^)  is  at  least  as  serious  as 
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ignoring  the  caHbration  error.  Actually,  using  only  three  standards  has  magnified  the 

confidence  factor  using  the  "calibration  band"  approach;  with  a  larger  number  of 

standards  in  a  calibration  curve  with  ETA-LIF,  the  value  of  aLoo  will  easily 

overshadow  the  calibration  error.    It  should  be  emphasized  that  even  when  using 

different  values  of  confidence  coefficients  eqns.  4.9  and  4.10  are  still  valid  indicators 

of  the  variability  of  the  calculated  LOD  due  to  the  imprecision  in  the  sample 

estimators. 

Conclusions 

In  this  chapter,  the  effect  of  two  parameters  (the  estimated  sensitivity  and 
blank  standard  deviation)  on  the  quality  of  the  calculated  LOD  as  an  estimator  of 
the  true  (unknown)  LOD  has  been  investigated.  When  optimizing  or  comparing 
analytical  techniques  by  using  the  value  of  LOD,  it  is  necessary  to  be  aware  of  these 
factors;  in  addition,  when  quick  estimates  of  the  LOD  are  needed  during 
optimization,  judicious  and  careful  use  of  eqns.  4.9  and  4.10  allow  an  evaluation  of 
the  quality  of  the  estimate. 

It  cannot  be  overemphasized  that  the  fluctuations  in  the  LOD  estimate  as 
presented  in  this  chapter,  and  chapter  4,  are  due  solely  to  the  imprecise  nature  of  the 
sample  statistics  Sb  and  ao  for  the  background  noise  and  analytical  sensitivity, 
respectively.  Long-term  drift,  day-to-day  or  even  hourly  fluctuations  in  the 
population  parameters  themselves  could  certainly  account  for  a  much  larger  variation 
in  calculated  LOD  values  for  a  given  analytical  procedure.  However,  there  is  clearly 
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value  in  knowing  the  estimated  magnitude  of  Slqd  as  discussed  in  this  chapter  since 

the  two  sources  of  variation  can  be  separated,  in  addition  to  assessing  the  effect  of 

the  measurement  protocol  used  to  estimate  the  LOD  on  Slqd- 


CHAPXFR  8 

INVESTIGATING  THE  PROCESS  OF  SINGLE  ATOM  DETECTION 

WITH  LASER-INDUCED  FLUORESCENCE 


Introduction 

The  purpose  of  this  chapter  is  to  illustrate  the  concepts  of  SAD  presented  in 
chapter  5  for  the  case  of  detection  by  cyclic  LIF.  Detection  by  cyclic  LIF  is  an 
example  of  a  nondestructive  detection  method;  the  advantage  of  a  technique  based 
on  nondestructive  detection  is  that  SAD  is  possible  even  in  the  presence  of  noise. 
The  general  outline  of  the  computer  simulations  used  in  this  chapter  were  presented 
in  chapter  6;  the  full  program  used  to  generate  the  raw  data,  as  well  as  the  signal 
processing  methods  used  and  the  algorithm  for  generating  random  numbers 
according  to  a  Poisson  distribution,  is  presented  in  the  Appendix. 

In  all  simulations  in  this  chapter,  unless  noted,  T^  =  8000  units  and  t^  =  20 
units.^  All  times  (and  corresponding  mean  flux  values)  given  in  this  chapter  are  in 
terms  of  arbitrary  time  units.  To  provide  a  basis  of  comparison  for  possible  SAD 
methods  based  on  flowing  atoms/molecules  through  Vp,  it  has  been  reported  that 
T,  =  10-100  MS  for  atomic  beam  experiments  [63,  76]  and  t,  =  1-10  ms  for  SAD 
experiments  with  flowing  solutions  [66,  68]. 


^or  T,  =  20  units  for  the  Cyhndrical  Vp  model 
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The  first  sections  of  this  chapter  will  illustrate  a  few  of  the  more  basic 
concepts  of  SAD  from  chapter  5.  The  basic  idea  of  e^  will  be  presented  for  the 
intrinsic-limited  case,  and  then  SAD  in  the  presence  of  noise  will  be  illustrated.  In 
addition,  the  precision  of  counting  atoms  will  be  calculated  theoretically  and 
compared  to  observed  values  from  the  computer  simulations.  Finally,  the 
consequences  of  a  variable  t^  will  be  demonstrated  through  the  use  of  computer 
simulations  based  on  the  Cylindrical  Probe  model.  For  the  demonstration  of  these 
concepts,  all  considerations  regarding  the  entry  times  of  atoms  into  Vp,  and  the 
detection  times  of  individual  photoelectrons,  were  ignored.  This  situation  would  be 
found,  for  example,  if  the  simple  integration  method  was  used,  in  which  the  total 
number  of  counts  during  T„  were  taken  as  the  signal.  All  values  of  X^,  0„  and  (p^, 
quoted  in  these  sections  will  be  based  on  simple  integration  as  the  signal  processing 
method,  with  T„  =  8000  and  t^  =  20  units.  Another  possible  scenario  where 
detection  times  and  atoms'  entry  times  into  Vp  would  be  unimportant  is  when  a 
pulsed  dye  laser  (pulse  lengths  10-100  ns)  where  used  (since  tj  <  <  tJ;  in  such  cases, 
t;  would  be  the  pulse  duration  and  the  number  of  photoelectrons  from  each  laser 
pulse  would  be  measured. 

The  application  of  SAD  theory  to  the  task  of  monitoring  atoms/molecules 
flowing  through  Vp  by  CW-LIF  will  be  demonstrated  through  the  use  of  the  full 
computer  simulations  to  generate  raw  data  files  (which  consist  of  the  detection  times 
of  photoelectrons,  which  are  due  to  either  atoms  passing  through  Vp  or  due  to 
background  noise).  The  photon  burst  and  sliding  sum  signal  processing  methods  will 
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be  used  on  the  data,  and  compared.   In  addition,  the  effect  of  T„  -  the  "scope"  of 

the  SAD  method  -  will  also  be  illustrated  through  these  simulations. 
Detection  Efficiency  at  the  Intrinsic  Noise  Limit 

If  the  value  of  0b  is  so  low  that  Mb  ^  o;  during  the  desired  T„  then  there  is  no 
reason  not  to  simply  integrate  the  signal  over  T„.  For  T„  =  8000  units  this  would 
correspond  to  0b  ^  l-^S  x  lO"''  counts  (unit  time)''  for  a  <  0.0014.  At  the  intrinsic 
noise  limit,  the  value  of  X^  is  1  single  count  during  T^.  Table  6  presents  the 
observed  detection  efficiency,  e^,  for  various  values  of  0^;  the  Simple  Model  was 
assumed,  with  constant  0^,  and  t,  =  20  units.  The  value  of  Np  was  fixed  at  1  atom 
for  8000  measurements;  when  the  signal  was  at  least  one  count  then  the  atom  is 
"detected"  during  T^. 

This  simulation  represents  a  test  of  the  Poisson  number  generator  as  much  as 
an  illustration  of  SAD  at  the  intrinsic  limit.  Table  6  presents  the  expected  and 
observed  values  of  e^,  i^  ,  and  a,.  As  expected  for  a  Poisson  detection  process,  the 
observed  mean  and  variance  of  the  signal  due  to  the  single  atoms  are  almost  equal 
for  all  values  of  0,.  SAD  is  achieved  at  the  95%  confidence  level  if  i,.  =  3 
counts/atom,  and  at  the  99.86%  confidence  level  if  i^  =  6.6  counts/atom,  as 
predicted  from  table  1  and  eqn.  5.19.  It  is  important  to  realize  that  the  detection  of 
individual  atoms  is  possible  for  all  these  conditions;  the  value  of  e^  as  a  measure  of 
the  effectiveness  of  a  near-SAD  technique  is  clearly  illustrated. 


Table  6 
Detection  Efficiency  at  the  Intrinsic  Limit 
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Sensitivity 

Observed 

Detection  efficiency 

True  mean 
(count/atom) 

<f>s' 

average 
(count) 

variance 
(count)^ 

expected 



observed 

0.1 

0.005 

0.0853 

0.0852 

9.52% 

8.16% 

0.5 

0.025 

0.5124 

0.5149 

39.35% 

39.99% 

1.0 

0.050 

1.0090 

1.0105 

63.21% 

63.33% 

3.0 

0.150 

3.0093 

2.9628 

95.02% 

95.28% 

6.6 

0.330 

6.6526 

6.4640 

99.86% 

99.84% 

9.0 

0.450 

8.9946 

8.8835 

99.99% 

99.99% 

The  value  of  X^  for  the  above  table  is  1  count;  an  atom  is  "detected"  if  a  signal  of 
at  least  one  count  is  detected  when  Np  =  1  atom.  The  value  of  <t>^  assumes  t-  =  20 
time  units.  The  values  of  the  expected  detection  efficiency  are  taken  from  Poisson 
distribution  tables. 


^The  signal  flux  has  units  of  count  (atom)'^  (time  unit)"^ 


'.    .-v'i. 
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SAD  in  the  Presence  of  Noise 

Once  the  level  of  noise  during  T^^  becomes  significant,  then  the  value  of  X^ 
must  increase  to  keep  the  probabiUty  of  false  positive  at  the  desired  level.  The 
inefficiency  of  the  simple  integration  signal  processing  method  in  the  presence  of 
noise  has  already  been  mentioned;  two  alternative  methods  will  be  illustrated  later 
in  this  chapter.  However,  if  0^  is  still  fairly  low,  and  discrete  samples  are  analyzed 
(low  values  of  T„),  then  the  simplicity  of  the  method  may  be  desirable.  The  process 
of  detecting  single  atoms  in  the  presence  of  noise  with  simple  integration  is 
illustrated  in  table  7  for  various  values  of  <p^  and  X^,  assuming  the  Simple  Model  with 
integration  over  T^  =  8000  units  and  with  t^  =  20  units.  All  values  of  X^,  Xg,  a  and 
B  shown  in  the  table  are  taken  from  tables  of  the  Poisson  distribution  [86]. 
Detection  at  two  different  significance  levels  are  illustrated,  represented  by  the  two 
different  pairs  of  values  for  X^  and  Xg  for  each  noise  level.  In  the  first  row  of  value 
of  Mb,  a  and  &  are  chosen  to  be  approximately  5%,  while  these  probabilities  are  kept 
around  0.14%  for  the  second  row. 

Choosing  the  signal  detection  limits.  Xj  and  X^:  The  signal  detection  Umit, 
Xj,  for  a  given  noise  level  is  chosen  so  as  to  keep  the  probabiUty  of  false  positive 
during  T„  (denoted  by  a)  below  a  pre-defined  level.  However,  due  to  the  discrete 
nature  of  X^,  it  is  not  possible  to  choose  a  at  any  level,  since  a  takes  on  discrete 
values  for  different  X^j,  The  exact  values  of  a  for  the  chosen  Xj  are  shown  in  the 
fourth  column  in  the  table.  The  nature  of  choosing  X^  at  low  levels  of  noise  is 
clearly  indicated:  it  is  rarely  possible  to  choose  a  to  be  exactly  the  value  desired  (5% 


Table  7 

Detection  and  guaranteed  detection  of  atoms 
in  the  presence  of  noise 
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^b 

h 

Xd 

a 

^g 

^ 

Observed  e^              | 

<^sti  =  Xd-/ib 

^sti=Xg-pb 

1 

0.00125 

4 

1.90% 

7.8 

4.85% 

57.19% 

95.58% 

6 

0.06% 

16 

0.14% 

54.78% 

99.88% 

2.5 

0.00313 

6 

4.20% 

10.5 

5.04% 

55.55% 

94.94% 

9 

0.11% 

21 

0.11% 

54.93% 

99.93% 

5 

0.00625 

10 

3.18% 

16 

4.33% 

54.68% 

95.50% 

14 

0.07% 

28 

0.13% 

52.79% 

99.89% 

10 

0.01250 

16 

4.87% 

23 

5.20% 

54.10% 

94.29% 

22 

0.07% 

39 

0.12% 

52.30% 

99.95% 

All  units  are  in  counts,  except  for  0^,  which  is  in  counts  (time  unit)" 
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or  0.14%);  instead,  upper  tolerance  must  be  defined,  and  a  chosen  to  be  below  these 

levels.  Naturally,  similar  considerations  hold  for  Xg  and  fl;  the  value  of  Xg  (which 

is  defined  in  terms  of  X^  by  the  procedure  outlined  in  chapter  5)  is  chosen  so  that 

fi  is  as  close  to  the  desired  value  as  possible. 

When  the  sensitivity  of  the  technique  was  such  that  the  signal  distribution  was 

centered  at  X^  (S/N  «  2-3),  the  observed  detection  efficiency  values  at  the  various 

levels  of  noise  were  those  given  in  column  six  in  the  table.  Due  to  the  asymmetry 

of  the  Poisson  distribution,  the  expected  value  of  e^  is  slightly  higher  than  50%.  It 

is  obvious  that  true  SAD,  as  defined  in  chapter  5,  is  not  achieved  at  these 

sensitivities.   On  the  other  hand,  when  the  sensitivity  is  such  that  the  mean  of  the 

distribution  of  Ij  is  exactly  Xg,  then  the  expected  value  of  e^  is  1-fi  (required  for 

SAD);  comparison  of  the  values  of  fi  and  the  observed  e^  show  that  this  was  indeed 

the  case.  ' '       !!'      • 


i     Counting  Precision  ,  ^  ^, 

Counting  at  the  Intrinsic  Limit  "  '       -     ^    J.":,^  '\ 

The  difference  between  the  signal  measurement  precision,  RSD„,  and  the 
counting  precision,  RSD^.,  when  a  given  number  of  atoms  pass  through  Vp  during  a 
measurement,  was  illustrated  in  figures  13  and  14  in  chapter  5.  Calculation  of  the 
theoretical  value  of  RSD^,  for  fixed  values  of  Np,  i^  ,  and  fi^,  is  represented  by  eqns. 
5.21-5.24.  The  theoretical  values  are  computed  in  the  case  of  Poisson  distribution 
of  Ij  with  a  simple  computer  program,  since  the  manual  computation  can  be  rather 
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tedious.   Particularly,  the  minimum  value  of  the  sensitivity  necessary  to  achieve  a 

certain  pre-defined  counting  precision,  is  of  interest  as  a  comparison  to  the  sensitivity 
necessary  for  SAD. 

A  computer  simulation  can  be  used  to  compare  the  observed  value  of  RSD^ 
with  the  theoretically  calculated  values,  since  in  a  computer  simulation  the  number 
of  atoms  probed  can  be  held  constant  for  any  number  of  measurements.  In  the 
simulations,  the  measured  number  of  atoms,  N„,  are  computed  for  each 
measurement  according  to  eqn.  5.21;  the  standard  deviation  of  N^  is  then  used  to 
calculate  the  observed  value  of  RSD^  and  comparisons  to  the  theoretical  value  can 
be  made.  Table  8  shows  the  results  of  8000  measurements  at  the  intrinsic  limit  for 
various  values  of  i^  and  Np,  along  with  comparisons  to  the  theoretical  values  for  both 
RSD,  andRSD^. 

A  number  of  interesting  observations  can  be  made  from  table  8.  First  of  all, 
in  all  cases  the  observed  and  calculated  values  of  RSD^  agreed  very  well.  The 
column  with  a  sensitivity  of  6.6  counts/atom  represents  the  requirement  of  SAD  at 
the  0.14%  confidence  level;  clearly,  RSD,  is  higher  than  0.1  (the  desired  precision) 
for  Np  <  16  atoms.  This  observation  illustrates  that  fact  that  a  sensitivity  which  is 
necessary  to  achieve  SAD  may  not  be  sufficient  to  count  any  number  of  atoms  with 
precision;  at  lower  analyte  concentrations  the  counting  precision  is  quite  high. 

Another  interesting  observation  involves  the  difference  between  RSD^j  and 
RSDj.  As  predicted  by  eqn.  5.20,  the  value  of  RSD^  decreases  as  Np  increases. 
However,  the  values  of  RSD^  at  low  numbers  of  probed  atoms  is  markedly  different 
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from  the  expected  RSD„,  and  increases  to  a  maximum  value  before  once  again 

decreasing  with  increasing  Np.   From  both  the  calculated  values  of  RSD^  and  the 

simulations,  the  value  of  (i,  ),  necessary  for  RSD,  <  0.1  is  approximately  35 

counts/atom  at  the  intrinsic  limit.    This  value  is  roughly  5-6  times  the  necessary 

sensitivity  for  true  SAD. 

Finally,  the  comparison  of  RSD„  and  RSD,  at  100  counts/atom  is  interesting 
smce  this  is  the  sensitivity  predicted  for  precise  counting  on  the  basis  of  RSD„  only. 
As  was  illustrated  in  figure  13,  the  counting  precision  at  this  sensitivity  is  much 
greater  than  predicted  by  RSD„  for  small  numbers  of  atoms,  since  the  overlap 
between  distributions  of  different  values  of  N  is  slight. 

Counting  in  the  Presence  of  Noise 

The  presence  of  increasing  levels  of  background  noise  during  T„  will  tend  to 
broaden  the  signal  distributions,  I„  for  fixed  values  of  N^;  hence,  it  is  expected  that 
RSD^  will  increase  as  well  for  a  given  sensitivity  at  the  analyte  concentration  level. 
Table  9  reflects  the  increased  signal  variance  in  increased  levels  of  (i,  \  necessary 
to  achieve  precision  counting  at  the  same  levels  of  ^b  as  shown  in  table  7.  It  is 
interesting  to  note,  through  comparison  of  tables  7  and  9,  that  the  value  of  (I^  )^  does 
not  increase  nearly  as  quickly  as  the  minimum  sensitivity  required  for  SAD  in  the 
"simple  integration"  case. 


Table  9 
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Counting  precision  in  the  presence  of  noise 


^b 

(counts) 

(is)c 
(counts/atom) 

Np 
(atoms) 

RSD,                     1 

Calc'd 

Observed 

1 

35 

3 

0.0995 

0.1001 

2.5 

36 

3 

0.0982 

0.0985 

5 

37 

3 

0.0976 

0.0991 

10 

39 

2 

0.0969 

0.0940 

Tlie  minimum  sensitivity  (counts/atom)  --  both  theoretical  and  observed  values  -- 
needed  to  precisely  count  atoms  at  all  levels  is  given  above.  The  value  of  RSD^,  and 
the  corresponding  N    represent  the  most  imprecise  counting  levels. 
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"Extra"  Variance  in  the  Cylindrical  Probe  Model 

In  the  discussion  thus  far,  it  has  been  assumed  that  0^  and  tj  were  both 
constant  for  all  atoms  which  enter  Vp.  However,  variance  in  either  of  these  two 
parameters  results  in  "extra"  variance  in  the  total  signal  distribution,  \,  as  predicted 
in  eqn.  5.13.  Simulations  based  on  the  Cylindrical  Vp  model  were  specifically 
formulated  to  investigate  this  process  by  making  the  residence  time,  t^,  a  variable 
depending  on  the  path  of  the  atom  through  the  cylindrically-shaped  probe  volume. 
In  such  a  situation,  the  total  variance  of  the  signal,  a,^,  is  given  by 

a]  -I  ^   (A^^<t),)2a2(r..)  ^^-^^ 

A  simulation  based  on  the  Cylindrical  Vp  model  at  the  intrinsic  limit  (simple 
integration  of  signal  over  T^)  was  run  8000  times  for  several  values  of  t^  for  tj  =  t^ 
(CW-LIF)  and  0s  =  1  count  (atom)'^  (time  unit)"^.  Table  10  shows  the  results  of  the 
simulations.  The  observed  values  s,^  and  s^(tj)  are  given  in  the  table;  the  value  of 
S2(tr)  was  used  as  an  estimate  for  a^(Q  in  eqn.  8.1  to  calculate  the  expected  total 
signal  variance.  The  calculated  values  of  a,^  compare  very  well  to  the  observed 
values  for  all  simulations.  There  was  good  agreement  for  simulations  with  varying 
levels  of  Np  and  0^  as  well  (although  the  results  are  not  shown  here). 

The  larger  signal  variance  due  to  variation  in  t,  can  be  clearly  seen  in  figure 
19  for  8000  loops  of  the  simulation  with  values  of  i^  =  12.35  counts/atom  and 
N  =  1  atom.  The  observed  distribution  of  Ij  (intrinsic  case)  and  the  Poisson 
distribution  which  would  be  expected  with  fixed  t,  are  shown  in  the  figure. 


Table  10 
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Comparison  of  observed  and  calculated  signal  variance 
for  the  Cylindrical  V   model 


Signal 
average 
(counts) 

Total  Signal  Variance,  a^ 

Residence  Time,  t^ 

observed 

calculated 

average 

variance 

1.00 

1.05                     1.08 

1.00                     0.08 

2.53 

3.04                     3.05 

2.51                     0.51 

4.98 

6.93                     6.98 

4.99                     2.01 

7.44 

11.70                   11.95 

7.47                     4.54 

10.09 

17.99                   18.20 

10.00                    7.97 

15.02 

33.49                   33.35 

15.01                   18.32 

19.92 

51.55                   52.00 

19.91                   32.06 

50.37 

246.02                 246.36 

50.26                  195.18 

The  value  of  0^  was  held  constant  at  1  count  (atom)"^  (time  unit)"^.    The  units  of 
t^  are  arbitrary  time  units. 


133 


</) 


00  o 

p 

O 

o 


a 

CO 


(Oiqissod    0008)    SdonounooQ 


«irf   a> 

^^ 
a  o 

S  73 

■  ts 


s 
S 


9i 
(U 
J3 


o  -a 

OO   »-" 

II  ja 

a  a 
?  -a 

""I 

•2  '^ 

■ft!    So 
•O    o 

MI'S  ^ 
«  «  « 


o  "  - 
a  '^  o 

■C  -c  ^ 


0\ 


OO 


134 
The  consequences  of  variable  tj  can  be  easily  understood  from  eqn.  8.1  in  light 

of  the  discussion  of  e^  and  RSD^  thus  far.  Obviously,  an  increased  variance  would 

decrease  e ^  and  increase  RSD^,  the  latter  particularly  since  the  total  signal  variance 

is  proportional  to  N^.  A  high  enough  sensitivity  to  result  in  a  true  SAD  technique 

for  constant  tj  would  perhaps  not  be  high  enough  if  the  value  of  a^(tj)  were  large;  in 

addition,  since  a^  is  proportional  to  <i>^,  it  might  be  that  a  rather  large  improvement 

in  sensitivity  would  be  necessary  to  achieve  SAD.    A  more  effective  method  of 

achieving  SAD  might  be  to  decrease  the  "extra"  variance  term  by  decreasing  ajCtj), 

even  at  a  cost  to  S/N.  Finally,  the  value  of  e^  as  a  FOM  for  near-SAD  methods  is 

even  more  clear,  since  S/N  does  not  reflect  this  "extra"  variance.  For  example,  the 

two  distributions  in  fig.  19  each  possess  the  same  S/N  ratio;  if  the  noise  level  were 

such  that  Xj  =  5  counts,  then  the  detection  efficiency  of  the  Poisson-distributed  \ 

would  be  significantly  higher  than  due  to  the  Cylindrical  Vp  model. 

Continuous  Monitoring  of  Atoms  with  CW-LIF 

The  scope  of  an  SAD  method  in  the  presence  of  noise  is  represented  by  the 
value  of  T^  over  which  it  is  possible  to  detect  single  atoms  at  the  required  detection 
efficiency.  Earlier  in  this  chapter,  it  was  seen  that,  when  the  simple  integration 
method  is  used,  the  required  sensitivity  necessary  for  SAD  increased  rapidly  as  0^ 
was  increased  (table  7).  However,  since  individual  atoms  have  residence  times  much 
less  than  T^,  it  is  not  surprising  that  the  sensitivity  required  for  SAD  increases  so 
rapidly  with  relatively  small  increases  in  0^. 
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The  rest  of  this  chapter  will  be  concerned  with  monitoring  atoms  under  the 

following  conditions:  T„  =  8000  time  units,  t^  =  20  units,  and  0^  =  0-05  counts/time 

unit.  This  background  level  results  in  an  average  of  400  counts  over  T„  in  the  simple 

integration  method,  and  a  sensitivity  of  approximately  130  counts/atom  would  be 

necessary  for  SAD  at  the  0.1%  confidence  level.  This  is  obviously  a  very  inefficient 

way  to  detect  single  atoms  which  interact  with  the  laser  during  T„.   In  chapter  5, 

three  alternative  signal  processing  methods  based  on  the  application  of  a  "time 

window",  t^  of  a  size  comparable  to  t^,  were  described.  Two  of  these  methods,  the 

photon  burst  and  the  sliding  sum  signal  processing  methods,  will  be  demonstrated  for 

the  case  of  near-SAD  methods  based  on  CW-LIF  to  detect  single  atoms  during  T^, 

under  these  conditions.^ 

The  Photon  Burst  Method 

False  positives  with  the  photon  burst  method:  choice  of  X,, 

The  theoretical  value  of  a  for  the  photon  burst  method,  for  a  given  level  of 
background  noise,  0^,  can  be  calculated  from  eqns.  5.25-5.28  if  the  form  of  the 
distribution  of  I,,  is  well  characterized.    In  any  event,  many  measurements  of  the 

•  :'■■'.  I 

blank  can  be  used  to  estimate  the  probability  of  type  I  error  at  any  value  of  Xj. 
Table  11  shows  the  comparison  of  the  observed  and  theoretically  calculated  average 
number  of  false  positives  during  T^  for  the  case  of  a  Poisson-distributed  background 


^The  Simple  model  will  be  assumed.  The  effects  of  variation  of  either  <p^  or  tj  would 
be  reduced  values  of  e^,  as  discussed  earher. 
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(i.e.,  shot  noise  limit).  The  agreement  between  the  two  sets  of  values  is  satisfactory 

at  all  values  of  X^.  It  is  important  to  notice  two  significant  facts: 

1.  The  "scope"  of  the  SAD  method  is  still  an  issue;  i.e.,  the  correct  value  of  X^ 
may  change  as  T„  increases.  With  t„  =  t^  the  average  background  count  in 
the  time  window  is  approximately  1  count.  A  value  of  X<j  =  5  counts  would 
be  sufficient  for  only  one  application  of  the  time  window  ~  i.e.,  T„  =  t^: 
guaranteed  detection  of  an  atom  (SAD)  during  its  transit  time  through  the 
laser  beam.  However,  even  with  a  method  based  on  a  time  window  such  that 
t„  =  t^  a  larger  value  of  X^  is  still  necessary  for  larger  values  of  T^,.  If 
Xj  =  5  counts  was  still  used  with  T„  =  8000  units,  then  a  single  atom  could 
never  possibly  be  detected  against  the  average  of  3.6  false  positives  due  to  the 

background  noise. 

.  ■  .y         -■-■'■       ■■  '■ .    r  >  i    ■   i  - 

2.  The  difference  between  the  simple  integration  method  and  the  photon  burst 

method  is  striking.    At  0^  =  0.05,  increasing  X^  to  9  counts  is  all  that  is 
necessary  to  reduce  the  probability  of  false  positives  back  to  approximately 
the  0.1%  level.   Compare  this  value  to  the  much  higher  values  of  X^  given 
above,  or  given  in  table  7  with  much  lower  values  of  0,,. 
Detection  of  single  atoms 

Simulations  of  CW-LIF  according  to  the  general  outline  given  in  chapter  six 
were  performed  to  construct  the  raw  data  file;  the  photon  burst  method  was  then 
used  on  this  array  of  detection  time  of  photoelectrons  during  T„.  The  value  of  Np 
was  fixed  at  1  atom  for  every  measurement  ~  in  other  words,  one  single  atom 


v^^.zs^- 


Table  11 
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Calculated  and  observed  false  positives 
using  the  "photon  burst"  signal  processing  method 


Xd 

Average  number  of  false  positives 

(counts) 

observed 

calculated 

5 

3.628 

3.797                i 

6 

0.693 

0.73195 

7 

0.107 

0.11881 

8 

0.014 

0.01663 

' 

0.001 

0.00203 

Distribution  of  background  counts  is  Poisson.  Conditions  of  measurements:  Tjj^  = 
8000  time  units;  t^  =  21  time  units;  0^  =  0.05  counts^nit  time.  The  observed 
number  of  false  positives  is  the  average  of  1000  measurements. 
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entered  and  exited  Vp  sometime  during  T„.   During  the  course  of  a  measurement, 

a  burst  was  considered  "detected"  if  there  were  at  least  9  counts  in  the  burst;  it  has 
already  been  established  that  it  is  very  unlikely  that  the  background  noise  level  would 
cause  a  burst  to  be  detected  at  this  level  during  T^. 

Table  12  shows  the  expected  and  observed  results  for  1000  measurements  at 
two  different  sensitivity  levels.  The  distribution  of  the  number  of  detected  bursts 
from  the  1000  measurements  are  shown,  along  with  the  expected  values.  The  first 
sensitivity  level  results  in  SAD  if  T^  =  t,  --  i.e.,  the  atom  would  always  be  detected 
during  its  transit  through  Vp.  The  use  of  this  sensitivity  in  the  simulation  clearly 
demonstrates  the  Umited  scope  of  SAD  at  that  sensitivity:  even  though  SAD  at  this 
sensitivity  is  possible  with  T„  =  t^  SAD  is  not  possible  at  this  increased  value  of  T„. 
The  second  sensitivity  corresponds  to  the  value  of  Xg  with  X^  =  9  counts  (B  =  0.1%), 
and  would  result  in  near-certain  detection  of  the  atoms  in  all  1000  measurements  if 
the  burst  values  followed  a  Poisson  distribution. 

In  comparing  the  expected  and  observed  values  it  appears  that  the  value  of 
e^  is  somewhat  less  than  expected  for  both  cases,  and  there  are  a  number  of  cases 
of  "double  counts,"  in  which  there  are  two  detected  photon  bursts  rather  than  the 
expected  one  due  to  Np  alone.  The  expected  number  of  false  positives  is  expected 
to  be  very  low  (1-2  for  1000  measurements),  as  has  already  been  established.  Thus, 
the  double  bursts  must  both  be  due  to  the  presence  of  the  analyte  atom  in  some 
manner. 
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Simulations  of  CW-LIF  with  photon  burst  signal  processing: 
Distribution  of  detected  photon  bursts 


Number  of 

<Ps'  = 

0.665 

</>/  = 

1.000 

detected 
photon  bursts 

Observed 

Expected 

Observed 

Expected 

0 

102 

53 

7 

1 

1 

862 

947 

867 

999 

2 

36 

0 

124 

0 

3 

0 

0 

2 

0 

Conditions:  Simple  Model  (fixed  4>^  and  tj.)  with  Tjjj  =  8000  time  units,  t^  =  20  time 
units,  t^  =  21  time  units,  X^j  =  9  counts.  A  photon  burst  must  have  a  sum  of  at 
least  9  counts  to  be  "detected"  above  the  noise.  One  thousand  measurements  were 
taken. 

^units  of  05  are:  counts  (atom)'^  (time  unit)'^ 
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Figure  20  explains  the  discrepancy  in  both  e^  and  double  peaks.  The  upper 
portion  of  the  figure  was  used  to  first  explain  how  the  photon  burst  method  works 
in  chapter  five.  Obviously,  in  the  particular  situation  shown  in  fig.  20(a),  if  the 
sensitivity  is  high  enough  for  SAD,  then  the  probability  is  high  that  this  particular 
burst  would  be  detected  since  t^  contains  all  the  detected  photons  due  to  the  atom. 
However,  a  single  spurious  background  count  can  result  in  the  splitting  of  the  photon 
burst  due  to  the  analyte  atom;  an  example  is  shown  in  fig.  20(b).  In  the  figure 
shown,  the  signal  due  to  the  atom  is  spread  out  over  two  separate  bursts  due  to  the 
single  noise  count.  In  such  situations,  depending  on  the  number  of  signal  photons 
emitted  and  the  exact  positioning  of  the  time  windows,  there  is  a  possibiUty  that 
either  both  bursts  containing  the  signal  are  detected,  or  that  neither  one  is.  At 
higher  sensitivities  (i.e.,  very  near  to  SAD  conditions),  then  it  is  more  likely  that  the 
divided  signal  will  result  in  a  "double  count"  than  in  no  detected  burst  at  all;  on  the 
other  hand,  at  lower  sensitivities,  there  is  a  greater  chance  that  the  analyte  atom  is 
not  detected  at  all.  The  results  shown  in  table  12  reflect  both  of  these  trends. 

Thus,  the  use  of  the  photon  burst  method  comes  with  these  complications. 
The  effect  of  the  "double  counts"  may  not  be  too  critical  if  the  analyte  concentration 
is  so  low  that  it  is  unlikely  that  two  analyte  atoms  enter  Vp  so  close  together.  In  such 
a  situation,  it  might  be  assumed  that  two  closely-spaced  bursts  such  as  represented 
in  the  figure  were  due  to  single  atoms.  However,  the  problem  becomes  worse  with 
both  increased  noise  and  increasing  analyte  concentration. 
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Sliding  Sum:  Peak  Height  Detection 

False  positives 

The  number  of  false  positives  ~  the  number  of  shding  sum  peaks  with  peak 
heights  of  at  least  X^  counts  ~  are  very  similar  to  the  corresponding  values  in  the 
photon  burst  method.  For  larger  values  of  Xj,  the  Gamma  distribution  can  be  used 
to  calculate  the  theoretical  value  of  a,  but  the  method  is  somewhat  tedious  and  is 
only  an  approximation.  Experimental  determination  of  a  is  the  best  method  to  set 
Xj.  With  Xj  =  5  counts,  an  average  of  4.621  false  positives  were  observed  with  1000 
measurements  at  <f>^,  =  0.05  counts/time  unit  and  t„  =  21  time  units;  this  value  is 
slightly  higher  than  the  corresponding  photon  burst  level.  At  X^  =  9  counts,  the 
calculated  value  of  a  is  0.00372;  there  were  three  false  positives  observed  in  1000 
measurements  when  X^  was  set  at  this  level.  This  value  of  the  signal  detection  limit 
was  used  for  the  sliding  sum  signal  processing.  j  .      i , 

Detection  of  single  atoms  .  ,; 

One  thousand  measurements  with  Np  fixed  at  1  atom  were  run  at  the  same 
two  sensitivity  values  as  was  used  for  the  photon  burst  method.  The  results  and 
expected  values  (same  as  in  photon  burst  case)  are  given  in  table  13.  Unlike  in  the 
case  of  the  photon  burst  method,  there  seems  to  be  no  problem  with  decreased 
values  for  e^,  indicating  that  peak  heights  tend  to  follow  a  Poisson  distribution.  In 
the  case  of  the  higher  sensitivity,  however,  there  are  still  a  few  double  counts.  In 
almost  all  cases,  these  extra  peaks  were  found  to  be  due  to  a  "shoulder"  on  the  peak 
due  to  the  analyte  atom,  in  which  the  value  fell  below  X^  and  then  increased  again 
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Simulation  of  CW-LIF  with  sliding  sum  signal  processing 
Distribution  of  detected  peaks 


Number  of 

"detected" 

sliding  sum  peaks 


0 


0,^  =  0.665 


Observed 


Expected 


58 


53 


942 


947 


0 


0 


0/  =  1.000 


Observed 


Expected 


1 


983 


999 


0 


0 


16 


0 


0 


0 


Conditions:  Simple  Model  (fixed  0^  and  t^)  with  T^^^  =  8000  time  units,  t^  =  20  time 
units,  t^  =  21  time  units,  and  X^  =  9  counts.  Sliding  sum  peak  heights  of  9  counts 
or  higher  are  considered  "detected."   One  thousand  measurements  were  taken. 


^units  of  (f)^  are:  counts  (atom)'^  (time  unit)"^. 
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due  to  a  single  background  count  outside  of  t^.  In  these  cases  it  was  obvious  that  the 

peak  was  due  to  only  a  single  atom  present  in  V  ,  and  the  double  counts  were  more 

an  artifact  of  the  signal  processing  algorithm  used  to  count  "detected"  peaks. 

Unlike  the  photon  burst  method,  the  sliding  sum  peak  maxima  can  be  easily 

used  to  form  a  picture  of  the  signal  distribution  due  to  single  analyte  atoms.  If  the 

analyte  concentration  is  kept  very  low,  then  it  is  unlikely  that  there  are  two  atoms  in 

V  at  any  one  time.  In  such  a  situation,  all  the  sliding  sum  peak  maxima  at  or  above 

Xj  will  be  due  to  single  atoms.   At  the  peak  maximum,  t^  will  be  centered  so  that 

the  window  contains  t^  so  that  the  distribution  of  maxima  at  or  above  X^,  provides 

insight  into  the  distribution  of  I^.     Figure  21  shows  such  a  distribution  for  a 

simulation  of  a  near-SAD  CW-LIF  method:  ij  =  9  counts/atom,  fiy,  =  1.05  counts 

during  t^  so  that  the  distribution  of  peak  maxima  should  have  a  Poisson  distribution 

(the  Simple  model  was  used)  with  a  mean  of  10.05  counts.  Both  the  observed  and 

expected  distributions  are  shown  in  the  figure,  and  match  quite  well.   In  a  normal 

experiment,  the  use  of  curve  fitting  techniques  applied  to  this  distribution  of  heights 

may  make  it  possible  to  obtain  estimates  of  both  ij  and  a^. 

Sliding  Sum:  Peak  Area  Detection 

The  possibility  of  using  peak  areas  of  sliding  sum  peaks  rather  than  peak 
heights  was  discussed  in  chapter  five.  Although  the  use  of  peak  areas  will  not  be 
demonstrated  with  the  computer  simulations,  the  possibility  should  be  included  as  a 
brief  discussion.  It  may  seem  intuitively  more  obvious  to  use  peak  height  detection 
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of  the  raw  data  transformed  by  the  sHding  sum  signal  processing  method,  but  in 

actual  fact  there  is  no  loss  in  information  when  using  the  peak  areas  as  the  signal 

rather  than  peak  heights  (see  eqn.  5.29),  and  there  is  one  important  advantage  to  be 

gained,  which  will  now  be  addressed. 

As  with  the  photon  burst  method,  the  peak  height  method  suffers  when  there 
is  overlap  between  atoms  residing  in  V  .  In  such  a  situation,  it  is  not  easily  to  count 
the  atoms  which  pass  through  Vp.  The  advantage  gained  by  using  the  peak  areas 
rather  than  their  heights  is  that  atoms  with  overlapping  residence  times  within  Vp  can 
be  easily  treated;  the  total  area  due  to  the  peak  encompassing  both  atoms  is  due  to 
the  total  sum  of  the  signal  due  to  the  atoms.  In  such  cases,  the  total  number  of 
atoms  which  contribute  to  a  given  peak  is  calculated  by  equation  5.21  for  Nn,;  all  the 
concepts  of  the  counting  precision,  RSD^,  also  apply  in  this  situation. 

Note  that  as  the  analyte  concentration  increases  further  --  to  the  point  where 
V  very  often  contains  at  least  one  analyte  atom  -  it  might  be  more  practical  to 
revert  to  the  simple  integration  method  once  it  is  certain  that  the  analyte  signal  is 
easily  detectable  above  X^. 

Comparison  of  Signal  Processing  Methods 

There  have  been  a  total  of  four  possible  signal  processing  methods  which  have 
been  mentioned  thus  far  in  this  work  which  can  be  applied  to  SAD  methods:  the 
"simple  integration"  method  (with  T„  >  t^),  and  the  three  time  window  methods:  the 
sequential  time  window  (with  t^  «  ty20),  the  photon  burst  and  the  sliding  sum  (peak 
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or  area)  methods.  It  is  difficult  to  state  unequivocally  which  of  these  is  "better,"  since 
the  performance  of  these  is  highly  dependent  on  the  conditions  and  the  requirements 
of  the  experiment  to  be  performed. 

The  process  of  simply  integrating  the  signal  during  T„  may  be  desirable  for 
small  discrete  sample  or  for  the  application  of  SAD  methods  to  larger  analyte 
concentrations,  as  has  been  mentioned.  The  advantage  of  simplicity  is  obvious;  in 
addition,  the  signal  can  be  processed  in  real  time.  However,  any  desired  temporal 
information  of  the  entry/exit  times  of  individual  atoms  into  Vp  is  restricted  by  the 
duration  of  T^. 

If  the  sensitivity  of  the  SAD  method  is  very  high,  then  the  possibility  of 
sequential  application  of  small  time  windows  over  T^,  should  be  investigated.  If  SAD 
is  possible  during  t^  then  the  atoms  can  be  monitored  in  real  time  as  they  cross 
through  V  .  Obviously,  this  is  the  best  situation  one  can  hope  for  in  an  SAD 
method. 

For  methods  in  which  the  sensitivity  is  low  enough  to  desire  to  increase  t^  to 
about  the  level  of  t^  then  two  possible  methods  have  been  addressed  in  this  chapter, 
and  some  of  their  strengths  and  weaknesses  have  been  demonstrated.  The  two  main 
advantages  of  the  photon  burst  method  are  the  fact  that  it  can  be  a  real-time  method 
[63]  and  that  it  is  fairly  uncomplicated.  However,  the  application  of  this  method  to 
higher  levels  of  analyte  concentration  or  background  noise  is  complicated  by  the 
presence  of  double  counts  due  to  single  atoms.   In  addition,  there  is  some  loss  in 
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detection  efficiency  of  near-SAD  methods  due  to  the  "spUtting"  of  atom's  photon 

burst. 

The  sliding  sum  method  is  a  fairly  powerful  signal  processing  method  to  apply 
to  near-SAD  techniques.  There  is  no  loss  of  detection  efficiency  and  information  can 
be  gained  about  the  distribution  of  the  analyte  signal  due  to  single  atoms.  In 
addition,  the  use  of  peak  areas  allows  the  application  of  this  method  to  atom 
counting  even  at  relatively  high  concentrations  of  analyte.  The  main  disadvantages 
are  that  the  method  is  computationally  more  complex  than  any  of  the 
others  -  particularly  for  long  measurement  times  ~  and  the  method  carmot  be  used 
in  real-time. 

Mention  should  be  made  of  one  other  signal-processing  method  which  has 
been  used  to  detect  photon  bursts  due  to  single  atoms:  autocorrelation  of  the  raw 
data  collected  over  T^  [66,  68].  Use  of  an  autocorrelator  is  relatively  simple  and  yet 
supplies  good  evidence  of  SAD  if  the  analyte  concentration  is  low  enough  that  there 
are  no  more  than  one  atom  in  Vp  at  a  given  time.  Another  advantage  of  this  method 
is  that  the  autocorrelogram  can  supply  a  good  estimate  of  Tj  if  T^  is  long  enough  and 
the  noise  is  low  enough.  However,  there  are  problems  in  attempting  to  use  the 
autocorrelation  function  in  attempting  to  measure  the  total  number  of  atoms  which 
have  passed  through  Vp  during  T„  [68]. 

There  was  no  intent  in  this  work  to  attempt  to  cover  all  possible 
signal-processing  methods  which  can  be  applied  to  the  raw  data  from  an  SAD 
experiment  using  CW-LIF.    There  are  many  more  complicated  signal  processing 
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schemes  or  moving  filters  that  might  be  useful.  The  intern  in  this  chapter  was  merely 
to  demonstrate  that  for  all  of  these  the  detection  theory  presented  in  chapter  five  for 
the  general  SAD  model  is  still  applicable  in  these  instances. 

Conclusions 

The  main  goal  of  this  chapter  was  to  use  simple  Monte  Carlo  simulations  to 
illustrate  the  apphcation  of  the  concepts  introduced  in  chapter  five  to  specific 
numerical  examples  of  SAD  and  near-SAD  methods.  What  has  become  evident 
during  the  course  of  these  experiments,  however,  is  that  such  simulations  can  be  very 
valuable  in  clarifying  some  of  the  subtleties  of  SAD  techniques.  Some  of  the 
examples  seen  in  this  chapter  are  the  presence  of  the  "extra"  variance  with  the 
CyUndrical  V  model,  the  scopes  of  SAD  methods,  and  the  behavior  of  the  two  signal 
processing  methods  chosen. 

However,  the  SAD  theory  presented,  and  the  computer  simulations  used,  are 
only  as  good  as  the  model  upon  which  they  are  based.  A  brief  look  at  some  of  the 
assumptions  inherent  in  the  general  SAD  model  as  presented  is  warranted.  The 
assumptions  of  Poisson  processes  for  both  the  background  noise  and  the  detection 
process  are  not  serious  obstacles.  It  is  quite  likely  that  in  many  situafions  the  first 
assumption  is  not  vaUd;  however,  some  assumption  of  the  form  of  lb  had  to  be  made, 
and  a  distribution  other  than  a  Poisson  merely  signifies  that  a  different  value  of  Xj 
must  be  chosen. 
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The  assumption  of  a  Poisson  detection  process  is  usually  valid  in  the  time 

frame  of  most  experiments.  However,  the  characteristics  of  Poisson  processes  must 

be  kept  firmly  in  mind.  The  presence  of  metastable  traps,  detector  "dead"  time,  or 

space-charge  effects  can  all  alter  the  distribution  of  the  detection  process.    In 

addition,  clustering  effects  or  other  interatomic  or  intermolecular  forces  can  alter  the 

distribution  of  Np. 

Perhaps  the  biggest  assumption  in  the  general  SAD  model  is  that  the  probe 
volume  is  well-defined.  In  most  analytical  experiments,  it  is  usually  the  case  that  the 
boundaries  of  the  probe  volume  are  rather  vague  due  to  either  the 
collection/detection  efficiency  of  laser-induced  events  in  the  atoms,  or  the  spatial 
characteristics  of  the  laser-beam  itself.  For  example,  a  common  method  of  describing 
the  radius  of  a  laser  beam  involves  the  use  of  the  1/e  distance  ~  i.e.,  the  distance  at 
which  the  intensity  decays  to  1/e  of  its  peak  value.  However,  the  laser  beam 
obviously  still  has  some  intensity  at  values  beyond  this  distance,  and  atoms  "outside" 
the  laser  might  still  be  detected.  Likewise,  in  LIF  experiments,  the  collection 
efficiency  of  most  optical  systems  is  a  function  of  the  distance  of  the  emitter  from  the 
central  axis  of  the  system,  and  any  abrupt  "line"  to  define  Vp  is  a  gross  simplification. 

The  concept  of  the  detection  efficiency  itself  is  obviously  only  as  well  defined 
as  the  probe  volume,  despite  the  fact  that  the  requirement  for  a  method  to  be 
considered  SAD  has  been  given  in  terms  of  this  FOM.  It  may  often  be  more  sensible 
to  speak  of  an  "effective"  probe  volume  within  which  almost  all  detected  events  are 
likely  to  occur.   In  such  cases,  it  may  be  impossible  to  achieve  very  high  detection 
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efficiencies,  since  it  will  always  be  possible  for  an  atom  to  skirt  the  edges  of  the  ill- 
defined  "effective"  Vp. 

The  assumptions  of  the  Simple  model  of  course  depend  upon  a  very  sharply- 
defined  probe  volume,  within  which  0^  is  almost  constant,  and  outside  of  which  0^  is 
effectively  zero.  Although  difficult  to  achieve  in  practice,  such  a  situation  might 
often  be  desirable  since  the  "extra"  variance  of  the  signal  will  be  reduced,  with 
consequent  improvement  in  e^  and  RSD^.  Some  of  the  possible  methods  to  reduce 
variance  include  laser  beam  expansion,  the  use  of  spatial  filters  and  apertures,  and 
good  collection  optics  (reduction  of  optical  aberration);  as  mentioned  earher, 
decrease  of  the  variance  of  the  signal  may  more  than  compensate  for  any  loss  in  S/N. 

Despite  the  rather  idealized  nature  of  the  general  SAD  model,  many  of  the 
concepts  of  the  theory  of  SAD  based  upon  this  model  are  of  value  to  the  laser 
spectroscopist.  Even  the  analytical  chemist  interested  only  in  ultra-trace  analysis  by 
laser-based  techniques  should  be  aware  of  many  of  the  concepts  presented  for  SAD, 
even  if  there  is  no  desire  to  count  individual  atoms,  since:  (a)  there  may  be  few 
atoms  or  molecules  of  analyte  being  actively  probed  by  the  laser  at  the  signal 
detection  limit;  and  (b)  the  background  noise  is  frequently  reduced  to  very  low  levels, 
near  the  intrinsic  limit.  Of  course,  the  analytical  or  physical  chemists  who  are 
actually  concerned  with  measuring  signals  or  properties  dealing  with  single  atoms 
struggle  with  concerns  which  do  not  arise  in  more  conventional  experiments;  the 
formulation  of  the  general  SAD  model  and  the  related  concepts  have  illustrated  this 
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fact  most  clearly,  and  might  be  used  to  advantage  in  the  evaluation  or  optimization 

of  any  (near-)SAD  technique. 
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APPENDIX 
RANDOM  NUMBER  GENERATORS 


The  following  two  functions  were  used  to  generate  random  numbers  in 

Microsoft  QuickBASIC.  Random  numbers  from  a  normal  distribution  were  obtained 

from  the  following  function,  which  accepted  the  desired  mean  and  standard  deviation 

as  the  mput:  ■  .  ;      . 

FUNCTION  RndNorm!  (mean!,  StanDev!) 
DO 
RandomA!  =  2!  •  RND  -  1!  ,.■     ' 

RandomB!  =  2!  *  RND  -  1! 

Radius2!  =  RandomA!  ^  2  +  RandomB!  ^2  -     • 

LOOP  UNTIL  Radius2!  <  1!  'Mod.  #2 

Deviate!  =  RandomA!  *  SQR((-2!  *  LOG(Radius2!))  /  Radius2!) 
RndNorm!  =  mean!  +  Deviate!  *  StanDev! 
END  FUNCTION 

This  next  function  returned  integers  from  a  Poisson  distribution: 

FUNCTION  Poisson%  (mean!) 

p!  =  EXP(-mean!) 

N%  =  0 

q!  =  RND 

DO  UNTIL  q!  <  p! 
q!  =  q!  *  RND 
N%  =  N%  +  1 

LOOP 

Poisson%  =  N% 
END  FUNCTION 
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